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ANNOTATION 


This  symposium  contains  materials  on  improved  verifi- 

k 

cation  te  cliques  for  individual  meteorological  instruments, 
more  precise  definition  of  their  parameters,  and  also  on 
methods  of  increasing  the  accuracy  and  precision  of  certain 
meteorological  investigations. 

The  symposium  is  intended  for  use  by  specialists  engaged 
in  precision  tool  engineering  and  meteorological  research. 


EXPERIMENTAL  APPLICATION  OF  THE  METHOD  OF  MATHEMATICAL 
STATISTICS  TO  MICRO STRUCTURAL  FOG  AND  CLOUD  RESEARCH. 


by  P.V.  Dyachenko. 

This  paper  describes  the  results  of  an  inquiry  into  the  problem 
of  fog  and  cloud  particle  size  determination  from  the  viewpoint 
of  statistical  probability. 

The  materials  contained  in  the  paper  make  It  possible  to  improve 
considerably  the  accuracy  and  precision  of  Investigations  relat¬ 
ing  to  the  micro  structure  of  fog  and  clouds. 

Folydi sper sity  is  a  general  and  integral  property  of  natural  fogs 

and  clouds.  So  far  this  fact  has  never  been  questioned  in  any  modem 

work  devoted  to  the  microphysics  of  fog  and  clouds.  Only  the  positions 

taken  by  I.P.  Smirnov  in  his  article  reviewed  in  detail  by  this  author 

on  an  earlier  occasion  present  an  isolated  exception. 

The  unquestionable  polydispersity  of  hatural  fog  and  clouds  entails 

the  requirement  for  a  consistent  approach  to  the  problem  of  slZ6  deter- 
the 

mination  of/particles  which  form  them  from  the  positions  of  statistical 
probability. 

The  statistical  inconsistency  manifested  by  some  authors  engaged 
in  microphysical  studies  of  clouds  is  responsible  for  a  considerable 
number  of  inaccuracies  and  errors. 

Houghton  /31/,  for  example,  has  Introduced  a  new  method  for  plot¬ 
ting  experimental  distribution  curves,  different  from  the  method  which 
is  considered  classical  and  is  employed  in  Europe  and  the  USSR.  This 
has  caused  not  only  Houghton  himself,  but  also  a  number  of  other  in¬ 
vestigators  who  have  accepted  his  conclusions  without  a  critical  analys¬ 
is,  to  develop  erroneous  theories  / 2,  16,  27/. 

A  physically  well  substantiated  selection  of  one  or  another  mean 

size  suited  for  a  given  problem  is  essential  for  micro  phy  sic  al  investi- 

the  practiced 

gations  of  fog  and  clouds.  Analysis  of  VHlBWHiV  application  of  t>he 
mean  sizes  in  fog  and  cloud  microphysics  shows  that  with  a  rare  exception 
/lo/  the  mean  sizes  are  improperly  used. 

FTJ)-TT-  C,  ~i.-S-a0//t  V—  / 
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Of  no  lesser  Importance  Is  the  matter  of  determining  values 

c/as<$ 

for  the  width  of  the  Interval  and  the  volume  of  the  sample  of 

the  investigated  set  of  particles  which  Could  assure  a  given  statistical 
accuracy  of  the  measurement  results.  The  absence  of  this  type  of  criteria 
has  lead  in  several  cases  to  the  appearance  of  statistically  inadequate 
data  (27),  (28). 

Special  attention  was  devoted  by  us  to  H.  KiBhler's  ideas  on  the 

relationship  of  masses  (33).  It  i  s  well  known  that  on  the  basis  of  a 

statistical  analysis  of  the  results  of  his  observations  Ktihler  has 

about  f 

formulated  a  hypothesise^  the  predominant  coalescence  of  equal-dze  drops. 

From  the  statistical  viewpoint  this  hypothesis  did  not  seem  to  give  rise 

to  any  objections.  It  is  known  that  in  his  work  Kohler  has  thanked  the 

famous  Ebglish  statistician,  K. Pearson,  for  his  assistance  in  dealing  with 

problems  of  mathematical  statistics  arising  in  the  development  of  Kohler's 

ideas.  Kohler's  theory  has  found  a.  large  number  of  supporters;  it  was 

presented  in  many  textbooks,  and  was  given  further  substantiation  in 

certain  scientific  papers  (16 ) ,  (22)  in  the  form  of  "the  law  of  multiple 

masses"  . 

In  our  work  we  were  oriented  by  the  guiding  ideas  of  Prof.  P.N. 
Tverskoy,  which  are,  possibly,  unique  in  our  own  and  foreign  literature^ 
(19).  He  doubted  the  correctness  of  Kiihler's  hypothesis  and  wrote  thus 
about  the  multiple  relationship  of  masses!  "Hence,  it  is  necessary  to 
look  for  some  other  po  s  sible  explanation  of  the  Indicated  observed  fact, 
or  —  In  questioning  its  realism  —  to  expose  it  to  verification  by  means 
of  new  observations  with  the  aid  of  more  accurate  methods  than  those 
employed  until  now"  (page  29). 

In  the  light  of  the  results  obtained  by  us,  it  became  clear  that 
Kohler  has  made  his  basic  errors  precisely  in  dealing  with  the  statistic¬ 
al  problem  where  his  views  remained  unchallenged. 
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An  important  connection  between  the  character  of  experimental  dis- 

tribution  and  the  form  of  the  investigated  clodiness  was  repeatedly 

established  in  the  papers  published  by  several  authors  (1,7,15,16 ,22, 

28,  and  others).  It  is  obvious  that  more  precise  data  in  this  respect 

will  be  forthcoming  only  after  the  accuracy  of  measurements  has  been 

subsequently  improved.  This,  on  the  other  hand,  calls  for  higher  standards 

in  the  technique  of  observation  and  the  treatment  of  the  results.  In 

fact,  the  objective  of  this  work  is  t>o  provide  the  technique  for  experi- 
parti^lt  size/  _ 

men  tal/de  termination  in  a  polydisperse  system  with  such 

statistical  substantiation  as  would  contribute  to  further  development  of 
sampling  methods  for  observation  of  natural  fo^  and  cloud  particles. 

Along  with  the  solution  of  this  principal  problem,  it  became 
possible  for  us  to  subject  to  critical  analysis  a  number  of  ideas, 
familiar  in  ?o£  and  cloud  micro  sphy  sic  s,  which  are  actually  wrong. 

In  addition  to  the  generally  known  lositions  of  mathematical  statistics,  we  have 
utilized  in  our  work  and  developed  in  application  to  the  requirements  of  the  problem 
at  hand  the  ideas  of  A.K.  Mitropol’skiy  (12)  on  the  envelopes  of  errors,  on  the 
character  of  experimental  distributions,  the  requirements  they  are  supposed  to  meet, 
and  so  on,  and  also  the  ideas  formulated  by  V.I.  Romanovskiy  (13)  on  the  mean  values  and 
their  meaning. 


The  dimensionless  method  was  used  in  presenting  all  the  materials  in  this  paper. 

As  a  result  of  this  we  have  succeeded  without  too  much  of  a  difficulty  in  introducing  a 
few  considerations,  a  move  precise  definiting  of  which  proves  to  be  rather  cumbersome  in  te 


terms  of’  conventional  treatment. 


The  present  article  reports  the  results  of  a  thesis 


study  l$l 


prepared  under  the  ' 


guidance  of  Professor  P.N.  Tverskoy  and  defended  in  1950. 

l.Mathod  of  Experimental  Distribution  Curve  Construction.  ' 

Whereas  in  the  USSR  and  in  Europe  it  is  customary  to  construct  size  distribution 
curves  for  fog  and  cloud  particles,  the  standard  practice  in  the  U.S.  is  to  deal 
in  terma  of  the,  so  called,  "volume 


© 
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distribution”.  The  latter  method  is  developed  mainly  in  Houghton’s 

works  (31).  Initially  (in  1932)  he  also  used  to  plot  size  distribution 
is  slfflcient  substantiation 

curves,  and  it/on}.y  later  (in  1938 );  and  without  ■■■■■■■■■■■■■II 
of  his  reasons  over 

■■■■^that  he  has  switched  to  the  new  method  of  curve  construction. 
This  fact  was  duly  registered  by  Ye.S.  Selezneva  (16)  and  B.V.  Kiryukhin 
(9)  in  their  papers.  However,  erroneous  ideas  still  continue  to  appear 
in  the  studies  of  some  European  authors.  These  ideas  are  similar  to 
those  revealed  by  Ye.S.  Selezneva  (16)  in  M.  Diem’s  writings  (27)  wherein 
the  author  assumed  that  the  American  fogs  consist  of  larger  particles. 

In  G-ertner's  book  (2)  one  finds  the  following  lines:  ”lt  is  difficult 
to  understand  from  Houghton’s  works  why  is  it  that  in  his  more  recent 
fog  investigations  he  has  found  the  predominant  size  of  drops  to  be 
about  ^  -  25  microns,  while  in  his  first  paper  the  highest  maxima  were 
registered  at  y  -  3*1,  6.2,  9*3  microns.  Whereas  his  earlier  measure¬ 
ments  are  in  agreement  with  the  measurement  results  of  other  investigat¬ 
ors,  no  such  correspondence  is  to  be  recorded  with  respect  to  his  new 
measurements”  (page  63). 

A  comparison  of  the  two  mentioned  methods  to  determine  their 
characteristics,  merits,  and  disadvantages,  no  doubt,  appears  to  be 
interesting.  Before  considering  the  differences  between  the  two 
methods  of  experimental  curve  construction,  we  shall  dwell  for  awhile 
on  certain  possible  methods  of  distribution  curve  plotting  for  a 
certain  set  of  real  particles. 

For  this  purpose  we  shall  first  analyse  a  unimodal  di  attribution 


of  particles  described  by  a  given  generalized  variate  X. . 
descrete  _ 

For  a  (■■  variate 'X,  which  is  numerically  equal  to 

(1) 

continuous 

where  t  -■■■■■■■  variate,  and  t  -  its  most  probable  dimensional. 


value,  a.  distribution  takes  place  which  for  a  set  of  N  particles  may 


f 


in  a  general  case  be  written  as 

dN  .  .  ,  .  . 

_=*«(*)<*:,  (2) 

occo&fzeuc*- 

where  —  is  the  ■■■■■■■■  probability  of  a  particle  with  the  variate 
N 

value  ranging  from  X  to  +  it}  n(X)  being  an  arbitrary  function  satis¬ 
fying  the  following  conditions: 

co 

a)  a)  /*«(t)dt=l;  ( 2a ) 

the  5 

b) / single  solution  of  equation  y1  >  -  0  is  root  X  -  1  correspond- 

d(r) 

ing  to  the  most  probable  value  of  t0 ; 

c)  function  n(X  )  is  continuous,  single-valued,  and  positive  over 
the  interval  0  <X  <co. 

Ratio  (2)  is  a  differential  expression  of  the  law  governing  the 
distribution  of  the  analyzed  variate  in  the  investigated  set  of  parti¬ 
cles.  On  drawing  a  random  sample  of  finite  volume  comprising  parti¬ 
cles  from  such  a  set,  we  shall  obtain  N^.  values  of  the  variate  which 
form  an  experimental  distribution.  From  the  physical  point  of  view 
this  distribution  may  be  considered  as  continuous  since  with  increasing 
volume  of  the  sample  it  may  be  as  near  to  such  a  distribution  as  desired. 
In  dividing  this  distribution  into  intervals  of  finite  width,  it  is 
possible  for  the  distribution  of  the  variate  to  obtain  frequencies  n^ 
or  relative  frequencies  ,  which,  correspond  to  the  mid-points  of  the 

Zni 

intervals.  Mi  experimental  distribution  obtained  in  this  manner  is 
usually  represented  graphically  in  the  form  of  a  stepped  curve  (or 
frequency  polygon).  Plotted  along  the  X-axis  are  the  variate  values, 
and  along  the  Y-axis  -  the  frequencies  and  relative  frequencies.  The 
distributions  which  satisfy  (2),  if  graphically  represented,  appear  to 
be  centered,  l.e.  their  modes  are  being  reduced  to  the  unit  value  of 

ef/j  c  ii  d  t&- 

the  variate.  .411  the  means  and  the  sizes  of  the  intervals 

may  be  measured  in  fractions  of  the  most  probable  variate  value  for  a 
given  population. 


S 


A  dimensionless  presentation  ,  therefore,  considerably  simplifies 
the  possible  comparison  of  Individual  distributions  regardless  of  the 
method  of  their  assignment.  Henceforward  we  shall  use  the  method  of 
dimensionless  presentation  and  resort  to  dimensional  statements  only 
when  necessary,  mainly  with  respect  to  the  final  results. 

The  basic  mean  values  which  interest  us  may  be  calculated  from  the 
following  formulae: 

the  arithmetic  mean  value  of  the  variate 

=  ,  (3) 

o 

the  mean  square  value  _ 

(4) 

the  mean  cubic  value 

X3  — 

In  restricting  ourselves  to  the  above  mean  values,  we  shall  pass 
on  directly  to  an  analysis  of  the  main  possible  methods  of  construct¬ 
ing  experimental  curves  which  may  be  plotted  in  investigating  a  set 
of  real  particles. 

Distribution  of  Sizes. 

If  in  formula  (1 )  we  assume 

t  =  p  =  i!  (6) 

where  r  is  the  radius  of  the  particles  In  a  set,  rQ  -  the  most  pro¬ 
bable  radius,  we  obtain  a  size  distribution 

3v~  =  «(p)rfp.  (7) 

Function  njp)  in  (7)  will  continue  to  satisfy  the  conditions  (2). 
The  method  of  experimental  curve  construction  adopted  in  the  Soviet 
Union  and  in  Europe  corresponds  to  such  a  distribution.  Here  the  pure 
variate  values  —  the  particle  radii  r  (or  diameters)  —  are  ploted 
against  the  X-axis,  while  the  frequencies  ©r  relative  frequencies  for 
the  experimental  curves,  and  the  probabilities  for  the  distributions 
given  in  the  differential  form,  are  plotted  against  the  Y-axis. 


V  fn(x)z3d'--kV^r- 


ip 


Expressions  for  the  size  distribution  means  may  be  obtained  from 

by 

(3),  (4),  and  (5),  by  replacing  in  them  't  p .  They  assume  the  form: 


the  arithmetic  mean  slate  distribution  radius 


UU 

P •'=/  "(P)PrfP  =  ^. 


the  mean  square  radius 


P2 


the  mean  cubic  radius 


: }/"/ w  <p>  p2  rfp = v  |/"5r  > 

•  n  •  • 


(8) 

(9) 


rt(rj)  rfidp  = 


,  i  ,3/sy? 
•  ''o  K 


do) 


More  detailed  information  on 


size  distribution  based  on  the 


method  of  experimental  curve  con  struct ion  re  suiting • from  this  distri- 
but  ion ,  as  well  as  the  data  on'  the  of  distribution  (7) 


will  be  discussed  below  in  the  corresponding  sections  devoted  to  this 
problem.  Here  we  shall  only  state  that,  the  mean' values  (8),  (9),  and 
(10)  have  an  entirely'  definite  physical  meaning  and  constitute  import¬ 
ant  physical  characteristics  of  distribution  (7). 

Distribution  of  Surface  Areas. 


a'trj.c  rc  Cel 


In  passing  over  to  a  new  variate  s,  numerically 

equal  to  the  ratio  of  the  particle,  surface  (sectional  area)  to  the 


mo st  probable  value  of  this  surface  sQ ,  we  shall  obtain  a  distribution 


of  surface  areas: 


an  ,  , 

-yy-  =  n(s)ds. 


(11) 


■  The  experimental  curve  for  this  distribution  may  also  be  represent¬ 
ed  in  the  form  of  a  stepped  curve  or  frequency  polygon.  As  in  the 
previous-  case  the  frequencies  and  relative  frequencies  are  plotted 
against  the  Y-axis.  Distribution  ,  (11 )  has  its  own  system  of  mean  values 


The  arithmetic  mean  particle  surface  is' equal  to 


Sj  =  j"  n  (s)  sds  ss 


4n  E  n{rf 


SoZrii 


(na) 


and  coincides  numerically,  as  it  is  easy  to  see,  with  the  surface  area 


of  a  particle  having  mean  quadratic  size  distribution  radius  (9). 

CL 


The  mean  square  surface  are  of  the  particle 
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(llb) 


and  the  mean  cubic  surface  area 


s2=l '/  f n  (s>  s'ds  =^Y 

o 

re  a  _ 

-  ,3/  r  t  \  ^  ~4n  i3/  s"^ 

S3=|/  ./  «(S)SS^^~]/ 


(11°) 


have  no  simple  and  obvious  physical  meaning,  so  that  their  use  as 

physical  characteristics  presents  no  interest.  The  mode  of  surface 

area  distribution  and  the  only  mean  s,  , which  has  an  obvious  physical 
for  skewed  dlstlbutlons 

meaning  do  coincide. 

The  method  of  experimental  curve  construction  based  on  the  utilization 
of  surface  area  distributions  has  found  no  application,  mainly  due  to 
the  absence  of  an  obvious  physical  meaning  of  the  means. 

Distribution  of  Volumes* 

In  a  similar  way,  through  replacement  of  v  for'C  (dimensionless 
volume),  it  is  possible  to  obtain  a  volume  distribution  for  a  given 


set  of  particles  from  ratio  (2) 


dN  /  V  . 


(12) 


with  the  same  conditions  for  n(v),  as  in  (2).  By  analogy  with  the 

preceding  cases  the  dimensionless  volume  vis  equal  to  the  ratio  of 

..  .  4*r3 . 


the  volume  of  particle 


to  its  most  probable  value  vQ  • 


The  volume  distribution  means  of  (12)  turn  out  to  be  equal  to 


—  C  47t2nlr? 

n(v)vdv^- 

0 

for  the  arithmetic  mean  of  the  particle  volume; 

v2=YTn^v2dv=  ikY  -^r 

0 

for  the  mean  square  volume,  and 

v3=y  J  n(v)v*dv  ^3 -y  -gsr 
for  the  mean  cubic  volume  of  particle. 


(13) 


(13a) 


(13b) 


The  arithmetic  mean  volume  v-j_  is  numerically  equal  to  the  volume 
of  a  particle  having  the  mean  cubic  size  distribution  radius  (10  )•  The 
remaining  means,  as  in  the  preceding  case|  ,  have  no  obvious  physical 
meaning  and  are,  therefore,  of  no  practical  interest.  The  meaning¬ 

ful  mean  volume  v-j_  may  be  obtained  only  through  computations  from 

8 


formulae  (10)  and  (13)  and  In  a  general  case  does  not  coincide  with 

mode  v  of  volume  distribution.  Plotting  of  experimental  curves 
o 

according  to  this,  and  the  preceding^  methods  is,  for  obvious  reasons, 
of  no  practical  research  interest,  since  there  is  no  advantage  in  pro¬ 
ducing  a  distribution  with  only  one  mean  value  which  is  present  in 
size  distribution. 

In  conclusion  we  shall  revert  to  the  method  of  experimental  curve 
construction  adopted  in  Houghton's  papers. 

Houghton 1  s  experimental  curves  differ  from  the  volume  distribution 
curves  analyzed  above  mainly  in  that  he  plots  the  linear  particle^ 
.dimensions  (diameters)  along  the  X-axis,  and  not  the  volumes.  This  fact 
does  not  essentially  affect  the  matter  since  it  is  the  volume  that 
functions  as  the  variate  in  Houghton's  distribution.  An  experimental 

<?/W  pc/'fcti./ 

curve  constructed  by  Houghton's  method  represents  a  distribution 

in  which  the  values  of  the  volume  as  the  variate  are  divided  into  inter- 
vals  whose  width  (Cn  terms  of  volume)  is  not  constant.  This  can  be 
readily  seen  if  along  with  a  linear  size  scale, one  should  plot  against 
the  X-axis  of  Houghton's  curve  also  a  scale  of  volumes  which  will  prove 
to  be  non-uniform. 

According  to  statistical  requirements  the  width  of  the  X-axis 

gwe-ru 

intervals  (classes)  must  be  the  same  in  any  distribution.  The  values 

A 

oftthe  distribution  variate  plotted  along  this  axis  must  be  of  equal 

magnitude.  In  his  method  Houghton  has  violated  this  important  condition. 

variaifc 

Quantities  which  have  the  dimensionality  of  the  volume) 

may  be  plotted  against  the  Y-axis  of  Houghton's  curves,  this  is  wrong. 

In  speaking  of  statistical  distributions  A*K.  Mitropol '  skly  (12)  makes 
this  quite  definite  statement:  "The  values  of  a  statistical  quantity 
represented  along  the  X-axis  are  as  a  general  rule  concrete  numbers, 
for  example,  centimeters  of  growth,  kilograms  of  weight,  and  so  on.  And 
by  contrast,  the  frequencies  represented  along  the  Y-axis  are  alway s 


abstract  numbers  indicating  how  many  times  that  or  another  value  of 


of  the  statistical  quantity  was  encountered".  In  Houghton's  curve s. 
whenever  one  turns  to  .the  frequencie  s  ( the  Y-axis  becomes  dimensional 
and  this  cannot  be  considered  correct. 

Therefore,  strictly  spoken,  Houghton's  curves  cannot  be  regarded 
as  experimental  from  a  statistical  point  of  view.  Furthermore,  the 
system  of  the  means,  as  it  logically  follows  from  the  above  discussion, 
is  absent  in  Houghton's  method. 

Instead  of  mean  volume  v^,or  the  corresponding  radius  r^,  Houghton 
utilized  a  mean  size  which,  as  it  will  be  shown  later  on,  is  the  mean 
size  of  distribution  (7)  having  no  direct  bearing  on  the  distribution 
of  volumes. 

In  the  light  of  the  aforesaid  we  arrive  at  the  obvious  conclusion 
that  the  method  of  experimental  curve  construction  in  the  form  in  which 
it  is  presented  in  Houghton's  papers  /3l/  is  in  its  very  essence  statist¬ 
ically  inconsistent.  At  the  same  time,  the  method  adopted  in  the  USSR 
and  Europe,  and  known  to  be  classical,  does  not  contradict  the  basic 
principles  of  mathematical  statistics.  This  is  the  reason  why  only 
problems  related  to  research  in  size  distribution  will  be  discussed 
further  in  this  text. 

2 .  Peakedness  and  Height  of  the  Distribution  Curves. 

The  character  of  experimental  distributions  produced  as  a  result  of 
natural  fog  and  cloud  particle  size  measurements  is,  as  a  rule,  described 
(also  in  discussing  the  peakedness  of  curves  in  which  we  are  interested) 
mainly  from  the  qualitative  point  of  view  /l 6/ .  We  know  only  of  one 
attempt  at  quantitative  estimate  made  by  Houghton  /31/.  He  was  determin¬ 
ing  the  proportion  of  the  total  amount  of  liquid  water  represented  by 
drops  contained  in  a  10-micron  band  centered  about  the  peafe>(  of  a 
relative-volume  distribution  curve.  The  width  of  the  band  consisting 


of  two  adjacent  intervals  —  one  of  which  is  situated  to  the  left,  the 
other  to  the  right  of  the  curve  maximum  --  was  made  arbitrarily  without 
any  substantiation  whatever.  Incest  of  the  measurements  were  carried  out 


by  Houghton  in  subdividing  the  volume  distribution  into  5-niicron  inter- 
l  -V\  4-ef'yn  & 

vals  di ameter )  .  Moreover^  he  has  found  that  the  "sharpness” 

estimated  in  this  manner  in  terms  of  percentages  varies  from  1 6  to  47, 
l.e.  increasing  3  times,  with  the  most  frequently  encountered  value  being 
30  %. 

One  might  in  a  similar  way  estimate  the  sharpness  of  tit  a  distri¬ 
bution  curve  constructed  by  the  method  of  plotting  experimental  curves 
as  adopted  in  the  USSR*  It  would  seem  sufficient  for  this  purpose  to 
determine  the  proportion  of  the  particles  comprized  in  the  band  of  fixed 
size  • 

Our  investigations  show  that  Houghton’s  test  has  no  direct  bearing 
on  the  evaluation  of  the  distribution  curve  peakedness.  The  discussions 
in  this  section  actually  represent  an  attempt  at  developing  a  new  peak¬ 
edness  test  suited  both  for  distributions  given  in  a  differential  form, 
and  for  experimental  distributions  expressed  in  the  form  of  finite 
differences. 

In  accordance  with  the  method  adopted  in  the  Soviet  Union  and  Europe 
for  constructing  emperlc  distributions,  we  shall  take  the  most  probable 
radius  rQ  of  size  distribution  (7)  as  the  basic  initial  size.  .411  other 
sizes,  including  thereby  any  means,  we  shall  measure  in  units  of  the 
most  probable  radius. 

Expression  (7)  will  from  now  on  be  considered  as  the^Mt  distribut¬ 
ion  law  of  a  population  consisting  of  an  infinite  number  of  particles 
from  which  a  random  sample  of  limited  volume  is  being  taken.  With  the 
exception  of  specially  qualified  cases,  distribution  (7)  will  be  assumed 
as  unimodal*  Consequently  function  n(^)  will  have  one  maximum  at  ^  -  1  • 

Let  it  be  mentioned  in  passing  that  function  n(^),and  only  it,  may 

in  full  determine  the  character  of  distribution  (7)«  The  distributions 

with  different  r0Ade scribed  by  one  and  the  same  function  n(^)  must  have 

a  common  peakedness.  There  can  be  no  inverse  effect,  i.e.  distributions 
pwwefj  _ 

with^tfqual  height  may  be  described  by  different  functions.  This 

II 


is  the  reason  for  the  bounded  character  of  the  peakedness  characteristic 


applicable  only  for  purposes  of  distribution  curve  comparison  in  terms 
altitude 

of  the  curve  at  the  maximum. 

dN 

It  is  easy  to  see.  that  in  (7)  ratio  ■jj-  represents  the  probability, 
dN 

while  is  the  probability  density  at  a  given  point  in  the  distribut¬ 

ion  curve.  The  probability  density  at  a  given  point  is  understood  to 
mean  the  probability  value  related  to  the  unit  value  of  variate  p  In 
our  case,  it  turns  out  to  be  equal  to  n(^)  for  distribution  (7). 

Fbr  unimodal  distributions  which  interest  us,  probability  density 
n  (f)  in  accordance  with  (6)  has  it$  own  single  maximum  at  point  p  «  1 
corresponding  to  the  peak  of  distribution  (7)»  For  a  characteristic 
providing  a  single-valued  determination  of  distribution  (7)  peakedness 

we  shall  use  the  numerical  value  of  the  probability  density  maximum  at 

altitude 

the  maximum  point  p  -  1.  This  value  then  will  be  the  of  the 

distribution  curve  equal  to 


!?■=> 


dNm 

Ndf 

dr 
'o 


(14) 


or,  respectively,  in  dimensional  form  when  dp  -  ^  ,  in  accordance 
with  (6) 

u  dNmr0 
Ndr  ’ 


(15) 


where  dNm  is  the  number  of  particles  with  sizes  ranging  from  p  -  1  to 
1  +  dp,  or  (which  is  the  same)  from  r  -  rQ  to  r  -  rQ  +  dr. 

Expressions  (14 )  and  (15)  hold  only  for  the  case  when  the  dis¬ 
tribution  is  given  in  a  differential  form.  If  we  have  to  deal  with  an 
experimental  distribution  in  which  the  subdivision  Interval  is  not 
infinitesimal,  but  finite,  then  formula  (15)  may  be  written  as 


H  m 

Nk\?  ’ 


(16) 


where  H  is  the  value  of  the  curve  altitude  averaged  for  an  interval  of 
finite  width  AP  and  is  numerically  equal  to  the  relative  ferequency 
density  maximum  at  this  interval;  nm  -  number  of  particles  measured 
in  the  region  of  Interval  bf>  maximum;  -  tbtal  number  of  particles 


/X 


(I6a) 


in  k  intervals,  or  the  sample  volume  numerically  equal  to 

i  —  I 

In  passing  over  to  dimensional  presentation,  we  obtain  a  final 

expression  suitable  for  empiric  distribution  curve  altitude  determinat- 

over 

ion,  when  this  distribution  has  been  measured  finite  intervals  A  r  - 


-  *i  -  r^. 


nmro 
Nkbr  ‘ 


(17) 


/ntp)/p*f 


Further  on  we  shall  analyze  the  errors  arising 
in  cdn  sequence  of  the  finiteness  of  and 
In  the  meantime  we  shall  note  only  that  if  the 


0  1  1  !■ 

I  2 

Fig.  1  Curve  altitude 


sample  volumes  are  sufficiently  large  and  the 


intervals  small,  formula  (17)  in  comparison  with 

C2.yp  j  i  c  j  t 

the  relationship  (14)  yields  such  results 


the  inaccuracies  of  which  may  be  conveniently  disregarded. 

Thus  the  numerical  value  of  the  curve  altitude  may  be  determined 
from  formula  (14 )  for  distributions  given  in  differential  form,  and 
from  (17)  for  experimental  distributions. 

In  the  first  case  we  calculate  the  probability  density  maximum. 


in  the  second  that  of  the  relative  frequency  density; 

ApEn;  ' 


(17a) 


If  we  construct  a  rectangle  the  height  of  which  is  equal  to  the 

distribution  curve  maximum  (Fig.l),  and  whose  width  is  assumed  to  be 

ratio  of  the  area 

equal  to  the  most  probable  radius  rQ ,  then  the  ^e 

constructed  rectangle  to  that  of  the  curve  (the  area  enclosed  between 

be 

the  curve  and  the  X-axis  will  for.  brevity  sake/ref  erred  to  by  us  as  the 
curve  area)  actually  characterizes  the  pointedness  of  the  distribution 
from  the  geometrical  point  of  view.  This  is  understandable,  when  one 
considers  the  fact  that  the  probability  density  in  point  of  the 

distribution  curve  is  numerically  equal  to  the  area  with  unit  base  and 

an  altidtude  equal  to  the  curve  ordinate  at  a  given  point,  provided  the 
area  of  the  entire  curve  is  regards  as  unit  area. 

13 


One  may  define  H  -  oo  for  distributions  with  the  gaximum  possible 
peakedness  (monodi  sperse  fog)  and  H  -  0  for  distributions  with  minimum 
peakedness  as  the  limit  values  for  curve  altitude  to  serve  as  parameters 
of  peakedness.  In  somewhat  forestalling  the  discussion,  we  would  like 
to  note  that ,  according  to  the  materials  covering  a  large  number  of 
experimental  distributions  relating  to  fogs  and  clouds,  the  most  probable 
curve  altitude  value  turned  out  to  be  equal  to  unity,  with  deviations 
rarely  exceeding  in  one  ore  another  direction. 

A  comparison  of  our  criterion  with  that  of  Houghton's  shows  that 

the  method  of  peakedness  evaluation  based  on  the  measurement  of  the 

proportion  of  particles, or  the  percentage  of  water  content,  available 

within  a  band  of  fixed  width  is  wrong.  The  resulting  quantity  in  this 

case  turns  out  to  be  dependent  on  the  most  probable  size,  and  the  size,  10 

11 

as  it  is  well  known,  is  subject  to  much  greater  changes  than  the  peaked¬ 
ness  which  interested  Houghton.  The  dependence  of  Houghton' b  character¬ 
istic  on  the  drop  sizes  is  obvious.  Indeed,  the  area  of  the  entire 
distribution  curve  numerically  expresses  a  quantity  which  is  dimension¬ 
less.  Portions  of  this  area  may  be  estimated  only  by  means  of  integrals 
of  the  folio  wing  form 

f  n(o)dp,  (I7b) 

where  ^q  and  ^2  are  the  finite  integration  limits. 

So  long  as  the  integration  limits  ^q  and  ^2  are  dimensionless,  the 
portions  of  the  distribution  curve  area  will  also  be  dimen sionelss,  but 
as  soon  as  pq  and  are  replaced  by  dimensional  limits,  the  curve  area 
portions  also  assume  a  dimensionality. 

If  we  should  use  a  dimensional  ratio 

H  nm 

-F^-N^r  (17°) 

to  estimate  the  distribution  curve  peakedness,  we  would  obtain  a  criter¬ 
ion  possessing  all  the  defects  of  the  ratio  used  by  Houghton.  The  latter 

ratio  also  appears  to  be  varying  over  a  larger  range  than  H,  in  view  of 
the  changing  most  probable  size. 


1 4- 


I 


Let  us  now  pass  to  the  relationship  between 
the  width  of  the  curve  and  the  height  as  introduced 
by  us. 

^■■■■iFor  the  characteristic  of  a  single 
peak  curve  it  is  customary  in  mathematics  to  use 
the  curve's  width  numerically  equal  to  the  differ¬ 
ence  of  abscissae  -p  1  in  which  function  n(^) 
succeeds  in  diminishing  e  times  as  compared  to  the 

Fig.  2. Width  of  curve 

maximum  value  (Fig.2).  A  comparison  of  parameter 

H  introduced  by  us  and  the  width  of  curve  j>"-  presents  an  obvious 
Interest.  Fbr  this  purpose  we  shall  revert  to  the  class  of  distribut¬ 
ions  originally  proposed  and  used  by  us  (5)  in  1950*  Its  form  is 

=  n  (P)  dp  =  Crjme  "  ‘  dp,  (18) 

where  C  is  the  normalizing  distribution  constant,  and  m  and  n  ar  the 
parameters  of  distribution. 

It  is  easy  to  show  that  the  normalizing  constant  C  is  Itself  a 
function  only  of  two  parameters  m  and  n. 


ffl  +  1 


dN  _  n  ( m  \  n  m  n  •  , 

N  r  +  lj  in)  P  1  (19) 


In  accordance  with  the  definition  for  the 

7'h&  Co  five  ft*-  T/Tootr 

class  of  distributions  (l9)/will  appear  to  bem equal  to 


tw  ne  hl  \  n  ,  , 

H  ~  r  fm  +  1\  in)  '  (20  ) 

\  n  )  curve 

width  of  class  (19)/we  shall  write  a 


In  order  to  find  the 
ratio  of  the  probability  density  maximum  to  that  value  of  probability 
density  which  is  e  times  smaller  than  the  maximum.  On  the  basis  of 
(19)  and  (2o)  this  ratio  is  equal  to 

(21) 


e. 


P  me 


— -p  n 
n  * 


In  transforming  (21),  we  arrive  at  311  equation  relative  top 

1  /  ,n  w+n\ 

p  =  en  l  m  ) .  (22 ) 

Equation  (22)  is  transcendental  with  respect  to  the  Independent 


variable  p  and  it  becomes,  therefore,  impossible  to  establish  an  accurate 
relation  between  H  and  p M  -p  1  in  a  general  form.  Important  for  us  is 
the  fact  that  for  class  (19)  curves  with  curve  altitude  close  to  unity 
(and  they  are  the  object  of  our  greatest  interest)  the  product  )H 

turns  out  to  be  a  practically  invariable  quantity.  Graphically  this 
fact  is  demonstrated  in  Table  1  data  calculated  for  distributions  with 
varying  m  and  n  from  formulae  (14 )  and  (22). 


Table  1 


n/n 

n  (p) 

?' 

P" 

P»  -  P' 

H 

(?"  -  P')  H 

1 

-4 

0,398 

1,775 

1,377 

0,83 

1,14 

2 

p2e  3 

0,45 

1,566 

1,116 

1,025 

1,15 

3 

_JL  • 

0,551 

! 

1,537 

0,986 

1,16 

1,14 

4 

p5e  2  P 

0,594 

1,481 

0,887 

1,28 

1,14 

5 

P6e-6p 

0,53 

1,69 

1,16 

0,972 

1,13 

6 

p3e— 3? 

0,39 

2,05 

1,66 

0,675 

1,12 

7 

p2e— 

0,30 

2,36 

2,06 

|  0,541 

i 

1,11 

/  4 


It  follows  from  Table  1  data,  that  the  following  relationship 

(p"  —  r/)Hss  1,13,  (23) 

takes  place  for  class  (19)  distributions  if  the  peakedness  is  moderate 

and  varies  around  unity.  The  meaning  of  this  relationship,  as  it  will 
be  easily  perceived,  consists  in  that  the  unimodal  distribution  curve 
area  is  proportional  to  the  product  of  the  curve's  width  by  its  height 
in  the  domain  of  the  maximum.  We  shall  note  in  passing  that  for  a 
normal  Russian  distribution  this  product  does  not  differ  from 
that  determined  by  us  and  is  equal  to  1.13* 

The  curve  altitude  as  a  parameter  of  distribution  curves,  including 
also  the  experimental  ones,  is  more  convenient  for  natural  fog  and  cloud 
particle  size  investigation  than  the  curve  width.  This  is  so,  because 
in  the  case  of  experimental  distributions  the  height  of  the  curve  may 
be  measured  with  a  considerably  lesser  number  of  errors  owing  to  the 

concentration  of  particles  about  the 


greater 
maximum. 

For  the  analitlcally  defined  distributions 


the  altitude  of  the  curve  may  also  be  determined  2>° 
in  a  simpler  way  and  more  precisely  than  its 
width • 

Should  the  investigator  be  interested  in  the 
width  of  the  curve,  then  it  would  appear  to  us 


f,S 


1,0 


more  effective  for  him  to  determine  it  by  method  os 

the 

of  calculation  from  approximate  formula* 


~_*2S  '22  • 

H  .  •.  H 


*  1637. 
:  1837  • 

’•  ’.’2455' 
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(23a) 


500  1000 

Fig. 3*  Peakedness  of 

experimental  curves 


ensuing  from  (17)  and  (23).  In  general,  however,  it  is  far  better  to 
deal  directly  with  the  altitude  of  the  curve. 

In  conclusion  we  shall  estimate  the  peakedness  of  experimental  curves. 
We  have  utilized  the  distributions  produced  by  a»M.  Borovikov  (1)  to 


serve  as  experimental  material.  These  curves  were  plotted  by  him  on  the 
basis  of  observations  carried  out  in  different  forms  of  clouds.  Shewn  in 

n 


Fig*  3  is  a  graph  Illustrating  the  effect  of  the  sample  volume  on  the 
magnitude  of  the  factor  calculated  from  formula  (17  )•  From  the  graph's 
data  it  appears  that  the  most  probable  value  is  that  of  H  -  1  with 
deviations  practically  not  exceeding  -  30 %  for  distributions  with  a  large 
sample  volume.  No  substantial  effect  of  the  cloud  particle  size  was  to 
be  registered. 

Thus,  the  height  of  the  curve  is  a  characteristic  which  is  constant 
for  natural  fog  and  cloud  particle  size  distributions,  provided  the 
volume  of  the  experimental  distribution  sample  is  sufficiently  large. 

This  characteristic  is  auxiliary  and  will  be  repeatedly  utilized  in 
the  further  development  of  the  discussion. 


3 •  Effect  of  the  Finltness  of  the  Interval  Width. 

It  is  well  known  that  any  large  population  consisting  of  an  infinite 
number  of  particle  s,who  se  sizes  may  interest  us, may  be  characterized  by, 

flkft~£cT/rtG,  BoTic.nl 

an  experimental  curve  (stepped  curve  or  frequency  polygon  )/ with  any  pre- 
assigned  accuracy,  If  the  size  measurements 

effected  in  the  form  of  random  sampling  are  carried  out  for  a  sufficiently 
large  number  of  particles,  with  the  width  of  the  interval  being  sufficient¬ 
ly  small.  In  spite  of  a  considerable  progress  achieved  in  the  practice  of 
natural  fog  and  cloud  particle  size  measurement,  both  these  requirements 
are  not  always  satisfied  to  a  sufficient  degree.  There  occur  measurements 
in  which  the  total  number  of  particles  -Ecii  is  too  small  (27),  (28), 
and  the  width  of  the  interval  is  unnecessarily  large  (29).  This  leads 
to  appreciable  errors  which  were  not  studied  with  respect  to  cloud  and 


fog  particle  distributions. 

p roblem  of  how  the/  affect s 

An  investigation  of  the/finltene ss  of  the  Interval  width  WKKKM  the 


construction  of  an  experimental  distribution  curve/ 


does  not  appear  to  us  to  be  fundament  ally  possible.  Hence  we  shall 


restrict  ourselves  to  an  analysis  of  a  distribution  curve  peak  statistic- 

ally  so  that  the  limitedness  of  the  sample  volume  may  be  totally 

,<i  disregarded. 


We  shall  consider  the  curve  altitude  (the  concept  of  which  was  defined 
In  the  preceding  section)  as  the  investigated  characteristic  whose  magni- 
tude  depends  substantially  on  the  width  of  the  interval*  \ 

If  we  denote  by  the  hight  of  distribution  (?)  Investigated  on 

the  basis  of  the  data  resulting  from  an  infinite  sample  by  way  of  divid¬ 


ing  it  into  Intervals  of  finite  width  then,  obviously. 


(24) 


where  H  is  the  height  of  the  curve  of  the  same  distribution  calculated 
precision 

from  the  formula  (14 )  in  the  assumption  that  the  width  of  the 

Interval  is  Infinite simal . 


Inequality  (24)  turns  into  equation 


(25) 


only  with  vanishing  A  p. 


This  means  that  in  particle  size  measurement  involving  any  unimodal 
distribution  (7)  by  subdivlndlng  it  into  Intervals  of  finite  width  Ap , 
no  matter  how  large  the  sample  volume,  there  occurs  a  distortion  or 
deformation  of  the  experimental  curve  as  compared  to  the  curve  of  a 
distribution  given  in  differential  form  and  corresponding  to  the  populat¬ 
ion  distribution.  The  immediate  objective  of  this  section  actually  is 
to  determine  the  magnitude  of  this  deformation  which  in  the  main  boils 
down  to  the  lowering  of  the  maximum  depending  on  the  width  of  the  interval. 
This  problem  is  equivalent  to  an  evaluation  of  the  error  arising  from 
the  utilization  of  equation  (25)  instead  of  inequality  (24). 

Remembering  that  nm  is  the  number  of  particles  measured  in  the 
vicinity  of  the  maximum  in  interval  Ap,  we  shall  find  the  expression  for 


nm  for  the  case  of  distribution  (7)*  Pursuant  to  the  definition 

l  +  4p 

nm=^N  f  «(P)4. 


(26) 


Then  formula  (7),  with  consideration  of  (26),  may  be  rewritten  in 


the  form  of 


i-Mp 

j  "(p)rfp 

l 


(27) 


where  HJSp  1  s  the  ctirve  altitude  of  a  distribution  distorted  owing  to 
the  finitness  of  the  interval  width. 

In  utilizing  (14)  for  H  and  (27)  for  H^p  we  shall  compose  an  express¬ 
ion  for  the  maximum  possible  relative  error  in  determining  the  curve 


altitude  : 


A  H  H  —  Hl  p 

H  ~~  H 


1+Ap 

J  «(p)rfP 

1 _ _ 

iPl«(p)lp  =  l  ‘ 


(28) 


AH  $n  the  basis  of  quite  obvious  considerat- 

/  T/-/e  /vitfG/v  n jU-p  e  c<a-  £/<rio,z.  ^  ~ 

ion  s. /depends  on  the  position  occupied  by  the  ordinate  of  the  distribut¬ 
ion  (7)  maximum  (the  distribution  mode)  in  the  arbitrarily  selected,  and( 
therefore^  also  arbitrarily  superimposed,  interval  of  width.  The  inte¬ 
gration  limits  in  (26)  are  so  selected  as  to  assure  that  formula  (28) 
should  correspond  to  the  case  of  maximum  possible  error  which  takes 
place  when  the  mode  of  distribution  (7)  coincides  with  the  boundary  of 

it 

two  ajacent  intervals. 

A 

Of  the  two  possible  intervals  xire  select  the  right-hand  one,  since 
measurements  of  the  fine  particles  in  the  left-  hand  interval  are  not 
always  reliable.  Moreover,  it  Is  interesting  to  investigate  the  case 
involving  intervals  the  width  of  which  is  comparable  <fco  rQ  and  larger 
thin  it.  This,  however,  is  possible  only  for  intervals  situated  to  the 
right  from  the  mode. 

Witja  reference  to  the  character  of  experimental  distributions  Ye.S. 

Selezneva  (16)  writes:  "Most  cloud  forms  are  characterized  by  distribut¬ 
or 

ion  curves  with  single  distinct  maximum.  The  structural  distinctions 

^  SK’e^ASe  is 

are  manifested  in  the  sharpness  of  the  maximum  and  the  nature  of  fgggKttKM 

sometimes  these  curves  are  narrow,  almost  symelsrical;  In  other  cases 

sice>AJ'ne  ss 

the  maximum  is  smoothed  out  and  the  curve  shows  a  sharp  in  the 

direction  of  large  drops"  (page  36).  To  obtain  a  dependence  1  (Af ) 

it  is  necessary  to  apply  (|§jp  concrete  distribution  law  to  (28)  and  then 


such  one  of  them  which  corresponds  to  the  most  symmetrical  experimental 
distributions.  In  the  capacity  of  such  a  law  we  have  utilized  the  size 


distribution 


=  2  p2e  1 


(29) 


which  constitutes  a  particular  case  of  class  (19)  for  m  -  2  and  n  ■  3« 
Ita  has  already  been  used  for  some  theoretical  calculations  by  N.S. 
Shishkin  (23)  in  the  form  of  volume  distribution  and  compared  with 
experimental  data.  This  distribution  in  terms  of  peakedness  (H  -  1*025) 

<>*re  00 

and  the  absence  of  distinctly  expressed  is  quite  suitable  for 

our  purposes. 

In  substituting  the  expression  for  porbabillty  density 


from  distribution  (29)  into  (28) 


_JL  * 

n  (0)  =  2p2e  3  P 


(29a) 


(30) 


as  well  as  the  maximum  probability  density  value 

l«(p)U,  =  2e_^.  (29b) 

After  integration  and  transformation  we  obtain 

-Si?  (i  +  ap  +  Ahh 

\H  2Ap_  l+e  V  3  !  (30  ) 

H  2\f 

Formula  (30)  gives  the  expression  for  the  maximum  possible  relative 

error  in  the  determination  of  the  curve  altitude,  which  arises  as  a 

result  of  the  finitness  of  the  interval  width.  This  formula  holds 

only  for  distributions  similar  to  the  model  distribution  (29)  in 

terms  of  peakedness.  G-lven  below  In  Table  2  is  the  dependence  of  the 

experimental  distribution  curve  compression  on  the  magnitude  of  interval 
cfhS/Ali 

where]  -jj-  values  are  calculated  from  formula  (30). 

The  results  of  the  first  half  of  the  table  for  the  case  when  0.1^ 

•  5  are  represented  In  Fig .4  In  the  form  of  a  curve  of  relative 
errors  accompanied  by  experimental  data  (points)  obtained  from  the 
e>qperiment al  curves  produced  by  A»M.  Borovikov  (l). 

As  it  may  be  seen  from  the  graph,  measurement  of  fog  and  cloud 
particle  sizes  with  utilization  of  intervals  comparable  in  width  with 
the  magnitude  of  the  most  probable  radius  may  lead  in  individual  cases 
to  considerable  errors  in  the  region  of  the  maximum.  On  the  contrary, 
if^^ff^0.2  -  0*25  the  errors  do  not  exceed  4  -  7 %•  Findelsen  (29) 


has  plotted  experimental  curve s  in  which  the  interval  width  exceeds 

the  magnitude  of  the  most  probable  size  rQ  •  For  some  of  his  curves 

2,  and  they,  therefore,  have  no  peaks  at  all. 

An  Ar  one 

In  selecting  the  width  of  the  interval  should  be  .guided 

ro 

by  the  data  from  Table  2  or  the  graph  in  Fig  .4  to  make  sure  that  the 

\% 

chosen  interval  not  be  toolarge.  But  there  exists  a  circum¬ 

stance  which  should  induce  the  observer  to  resist  being  carried  away 
in  selecting  ■■■■■I  overly  small  intervals.  Indeed,  if  the  experiment¬ 
al  distribution  is  represented,  as  this  is  usually  the  case,  in  the  form 
of  a  stepped  curve  or  frequency  polygon,  then  the  smaller  the  interval, 
the  more  preaise  is  the  position  of  the  maximum  if  the  sample  volume 

is  sufficiently  large  so  that  the  error  resulting  from  the  insuficient 
filling 

statl stical  of  the  interval  becomes  negligibly  small  at  least 

for  the  intervals  near  the  experimental  curve  maximum.  Under  these 

conditions  the  position  of  the  maximum  cannot  be,  in  principle,  defined 

with  an  accuracy  exceeding  one-half  of  the  Interval  width.  However, 

whenever  the  nj_  values  in  the  vicinity  of  the  experimental  curve  maximum 

are  small  —  and  this  is  possible  at  any  arbitrarily  selected  large 

Nk  so  long  as  the  width  of  the  interval  is  sufficiently  small  — 

the  advantage  in  ikK  determining  the  maximum's  position  resulting  from 

the  reduced  interval  width  turns  out  to  be  but  apparent.  This  is  so 

because  with  a  decrease  of  the  interval  width,  the  insufficient  statistlc- 
fllllng  ' 

the  maximum  and  its  environment  begins  to  manifest  itself 

and  reduces  to  nought  the  effect  resulting  from  the  reduction  of  the 

Interval  width.  In  the  next  section  we  shall  consider  the  problem 

due  to 

concerning  the  error  in  n<  determination  the  insufficient 

filling 

statistical  or  (which  is  the  same)  to  the  finiteness  of 

the  sample  volume. 

Consequently,  the  width  of  the  interval  cannot  be  selected  arbitrar¬ 
ily,  and  must  be  comensurate  to  the  value  of  the  most  probable  radius  rQ . 


v.  rtco  W'wa 

TheVtiest  values  for  the  width  of  the  interval  in  size  measurements  of 
natural 


fog  and' cloud  particles  are  -  0.2  -  0.25,  which  still  do  not 
cause  any  noticeable  deformation  of  the  distribution  curve ,  and  at  the 

■k  maximum  statistical 


pame  time  assure  the 
concentration 

^■■■ftin  each  of  its  Intervals  taken  sepa¬ 
rately  . 

The  problem  of  the  interval  width  is  quite 
evidently  connected  with  the  question  relative 
to  the  number  of  intervals  in  an  experimental 
distribution.  The  necessary  number  of  intervals 
at  Apsi  0.2  -  0*25ro  appears  on  the  average  to  be 
equal  to  12-16  for  experimental  distributions 

a tox¬ 
in  which  the  maximum  variate  value  does  exceed 

3-4  rQ  .  Our  estimates  on  the  number  of  required 


Intervals  for  all  practical  purposes  coincide 


Fig .  4 .  Effect  of  the 
interval  width. 

with  the  results  produced  earlier  by  Prof.  A-K.  Mitropol '  skiy  by  a 


different  method. 


Table  2. 


Ap 

0,1 

0,2 

0,3 

0,4 

0,5 

0,6 

0,7 

0,8 

0,9 

„  AH 

100  -7T 

0,6 

3,8 

8,5 

1 

14,1 

20,6 

27,2 

33,9 

40,0 

45,5 

Ap 

1,0 

1,5 

2 

4  j 

8 

16 

32 

64 

CO 

< 

50,4 

66,7 

75,0 

87,5 

93,8 

1 

97,7 

98,5 

99,2 

100 

79 
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4  •  Influence  of  the  Finltne.,ss  of  a  Random-Sample  Volume . 

The  effect  of  the  volume  finitness  of  a  random  sample  in  particle- 
size  measurments  in  fog  and  clouds  as  polydi spersed  systems  may  be  mani¬ 
fested  in  such  a  possible  distortion  of  measurement  results  that  this 
distortion  could  be  considered  as  an  accidental  error  whose  magnitude 
must  be  determined  with  recourse  to  methods  of  probability  theory  and 
mathematical  statistics. 


Let  n^,  as  before,  stand  for  the  frequencies  characterising  the  i 
of  particles  comprized  in  a  finite  number  k  of  intervals  with  the  total 
number  of  particles  -  ^n^  • 

We  shall  consider  one  such  arbitrarily  selected  interval. 

It  is  obvious  that  the  proportion  of  particles  p  of  size  is  un¬ 
known  to  us.  Therefore,  in  utilizing  the  experimental  data,  we  are 


iHer 


n^ 


as 


compelled  to  regard  the  magnitude  of  relative  frequency  -  P 

an  approximate  value  of  the  proportion  of  particles  p  =»  n(p)Ap\  ’which  is 
unknown  to  us,  i.e.  to  assume  in  compliance  wi th  the  law  of  large  numbers 
that 

P  P1 


To  estimate  this  approximate  value  it  is  necessary  to  select  before¬ 
hand  the  value  of  reliability  0(  so  close  to  unity  as  to  be  able  to  dis¬ 
count  the  event  with  probability  1  -  eC  as  practically  impossibe,  and  to 

determine  such  deviations  of  the  unknown  proportion  of  particles  from 

ni 

the  observed  relative  frequency  p1  -  which  may  be  expected  with 

k 

probability^  in  the  case  of  a  volume  random  sample  taken  from  an 
infinite  general  population. 

It  is  known  from  the  probability  theory  that  the  magnitude  of  these 
deviations  may  be  determined  from  formula 

(31 ) 

where  €>p  is  the  standard  propprtion  p;  ■  <P  (x)  is  the  probability  that 
the  unknown  proportion  falls  between  the  limits 

P'  -  x°P<P<P'  +  xip.  (3!a) 


The  numerical  values  of  parameter  x,  depends^®  on  the  selected 

A 

reliability,  are  determinable  from  the  table  °/ 

0 

which  are  usually  provided  in  the  reference  section  of  textbooks  in 
the  theory  of  probability. 

Expression  (31)  is  hardly  suitable  for  practical  application 

_on_ 

since  the  errors  appears  to  depend  not  only  4^01  sample  volume  N^, 

but  also  on  the  magnitude  of  the  relative  freauency  —  . 

Nfc 

Formula  (31)  yields  an  expression  for  relative  error  Dj_  more 


convenient  for  practical  utilization 


(31°) 


In  substituting  p'  -  — ^  into  the  above  expression,  we  obtain 

Nk 

D^-^Vl-p'.  (32) 

V  n, 

From  (32)  it  follows  that  for  the  curve  sections  of  an  experimental 
or  rxKdora  sampling  distribution,  where  p'  ■  ^  ia  at  its  maximum  and 

Nk 

comparable  to  unity  (this  is  possible  only  for  excessively  large  inter¬ 
vals),  the  relative  error  value  in  measurement  of  depends  both  on 
n.  itself,  and  on  relative  frequency  p*  -  l.e.  in  the  final  count 

also  on  the  number  of  particles  Nk  on  the  measurement  of  which  the 
investigated  experimental  distribution  is  actually  based. 

In  the  remaining  points,  or  more  precisely  in  the  experimental  curve 
intervals  where  p 1  ^  1  or  n i  ^  Nk  (and  in  the  case  of  small  intervals 
and  a  moderate  peakedness  also  over  the  entire  length  of  the  experimental 
curb6 ) ,  the  radio  and  1-p’  hardly  differs  from  unity,  and  the  square  root 


from  this  quantity  is  the  more  so  nearer  to  unity. 


Thus,  for  an  experimental  distribution  with  arbitrary  ( 
shape  of  a  single-peak  .  Ar 


peakedness 


the  condition 


curve,  for  with  intervals 
n,CNk, 


so  small  that 


(33) 


is  secured  everywhere,  including  the  maximum  itself,  formula  (32)  may 
be  rewritten  in  the  following  final  form 


434) 


Dl  Vnt  : 


On  the  basis  of  (16 )  condition  (33)  for  the  case  n-j_  -  nm  may  be 


rewritten  in  the  form  of  inequality 


H&f>  <S  1, 


which  holds  for  arbitrary  unimodal  distributions.  For  a 
unit  peakedness 

(33)  assumes  the  aspect  of 

ApCl. 


distribution  with 


(35) 


From  formula  (34)  it  follows  that  the  relative  error  in  determining 
nj_  for  an  arbitrary  distribution  with  n^  Njj  is  not  the  same  for  differ- 

IAJl<recLA 

ent  experimental  curve  intervals,  as  thisYwas  to  be  expected.  At  a  given 

reliability^  error  becomes  a  function  only  of  nj_  —  the  number  of 

particles  contained  in  the  given  interval  in  which  we  are  interested  — 

on 

and  is  independent  of  Nk  -  the  number  of  particle s  0MI  the  measurement 
of  which  the  investigated  experimental  distribution  was  actually  based. 
Thus^we  come  to  the  conclusion  that  the  relative  error  in  determining nj_ , 
if  the  intervals  are  sufficiently  small,  does  not  essentially  depend  on 
the  distribution  law,  nor  on  the  volume  of  fcandom  sample.  Consequently, 
formula  (34)  is  applicable  to  any  given  experimental  distribution,  pro¬ 
viding  condition  (33)  is  satisfied 

In  analyzing  the  errors  in  an  interval,  the  numbers  n.  —  if  they 
particle  concentration  In/  which  are/ 

denote  the  measure  o f  interval s,/ small  in  comparison 

with  the  most  probable  radius  —  may,  in  approximation,  be  considered 

as  independent,  and  the  polydisperse  fog  regarded  as  consisting  of  a 

finite  number  of  independent  monodisperse  components.  In  passing,  we 

shall  note  that  formula  (34)  gives  higer  relative  error  values  if 

i  5 

condition  (35 )  viol  ated.  In  this  case  the  excess  value  is  greater  for 
the  maximum  zone  and  smaller  for  the  remaining  intervals  of  the  experimental 
di stribution . 


It,  therefore,  becomes  apparent  that  the  errors  calculated  dfh  the 
basis  of  formula  (34)  —  if  the  Intervals  i  are  sufficiently  small  — 


are  applicable  to  any  experimental  distribution  and  to  any  of  its  inter¬ 
vals  with  ail  accuracy  sufficient  for  all  practical  purposes,  if  the  peak¬ 
edness  of  the  investigated  experimental  distribution  differs  littl6  from 


| 


fg  unit  peakedne  ss  • 

particle  concentration  in  an/ 

The  relative  error  in  determining  the/InTervalTJBjHgjjjjj^'who se 
magnitude  i s  numerically  determinable  by  expression  (34),  is  a  function 


of  frequency  ni 


Hence,  it  is  not  identical  for  the  different  intervals 


or  regions  of  one  and  the  same  experimental  curve. 

In  conclusion  we  shall  dwell  briefly  on  the  matter  of  reliability. 
From  a  statistical  point  of  view  there  necessarily  must  be  a  correspond¬ 
ence  between  any  error  and  the  reliability  with  which  it  is  estimated. 

The  very  essence  of  the  reliability  concept  indicates  that  it  should  be 
so  close  to  unity  in  the  solution  of  practical  problems  that  the  events 
with  1-oi  probability  could  be  considered  as  practically  unfeasible. 

V.I.  Romano vskiy  (14)  remarks  that  depending  on  the  practical  requirements 
e<is  being  assumed  equal  to  0*95,  0*99,  or  0  *999  •  For  our  calculations 
we  have  used  the  value  of  <*  -  0»95» 

In  all  experimental  stldies  known  to  us  which  deal  ivith  mathematical 


statistics  in  relation  to  the  structure  of  natural  fog,  clouds,  and  rain 
(32),  (34),  (29),  (10),  as  well  as  those  devoted  to  research  in  disperse 
system  concentrations  (3),  the  ft  Rirftmrtit  A?£/W  error  is  usually  discussed 
in  estimating  random  errors.  Important  is  the  fact  that  reliability  e< 

(a  definition  of  this  concept,  regrettably,  is  lacking  in  most  works)  in 
this  case  turns  out  to  be  congruent  with  its  complement  with  respect  to 
unity.  Thus,  for  a  deviation  distribution  obeying  the  normal  law, 
reliability  0(,  which  characterizes  the  arithmetic  mean  error,  appears  to 
be  equal  to 

o(  -  0.68,  and  1-ot  -  0*32. 

This  means  that  the  probability  of  measurement  with  an  error  greater 
than  arithmetic  mean  error,  turns  out  to  be  considerable  enough  and 


comparable  to  the  reliability  itself.  In  this  case,  on  the  average,  out 

0 

of  three  measurements  two  may  fit  in,  while  the  third  may  exhibit  a 

deviation  exceeding  the  arithmetic  mean  error.  In  our  case  the  reliablli 

ty  is  so  close  to  unity  that  l-©(  -  0«05»  To  utilize  the  arithmetic  mean 

error  without  due  consideration  of  the  reliability  which  characterizes 

done 

it,  as  this  vasflH|by  Ktihler  (32),  is  obviously  a  mistake.  Etihler  and 
Niderdorfer,  without  paying  any  attention  to  the  long  since  known  concept 
of  reliability,  treated  the  arithmetic  mean  error  in  analyzing  the 
amplitudes  of  a  multipeak  curve  as  if  this  error  had  a  unit  reliability. 
The  observation  results  and  the  conclusions  drawn  from  them  by  BUhler  are 
therefore,  trustworthy  only  to  the  extent  to  which  the  events  with  pro¬ 
bability  1-oC-  0*32  may  be  attributed  practical,  unfeasibility. 

The  utilization  by  Kohler  of  arithmetic  mean  errors  characterized 
by  small  reliability  0^  value  actually  constitutes  the  basic  reason  for 

f.  aJLX-O-cjj 

the  mistake  he  made  in  treating  the  experimental  data.  This  thus 

caused  him  to  formulate  the  physically  untenable  hypothesis  and  is  some¬ 
times  mistakenly  referred  to  as  Kohler's  multiplicity  law. 

A  similar  error  was  in  its  time  made  by  Defant  (26),  and  later  by 

Crap'S.  . 

Niderdorfer  (34)^,  with  regard  to  the  size  distribution  of  rain 

5*  Sufficiency  Test  of  a  Random-Sample  Volume. 

In  the  preceding  sections  we  have  found  a  method  of  estimating  the 

accuracy  of  an  arbitrary  experimental  distribution  curve  in  any  of  its 

intervals,  if  the  interval  is  sufficiently  small* 

The  described  method  may  prove  to  be  useful  only  in  analyses  of  the 
filling 

statistical  for  the  already  plotted  experimental  curves.  In 

those  cases  when  the  investigator  approaches  the  problem  of  experimental 

curve  construction,  it  is  desirable  to  obtain  the  estimates  for  such 

total  number  of  particles  in  measuring  which  the  required  statistical 
filling 

could  be  assured  at  the  maximum. 

Until  now  the  total  number  of  particles  N^  -2»i  in  measurements  of 

At 


fog  and  cloud  particle  sizes  is  known  to  have  been  absolutely  random.  Its 
value  varied  over  a  wide  range  with  totally  unsubstantiated  limits. 

Houghton  /3l/  supposed  that  a  "good"  single-peak  experimental  curve  in 
measuring  the  particle  sizes  present  in  natural  fogs  and  clouds  may  be  pro¬ 
duced  practically  always  when  the  number  of, particles  measured  to  yield  the 

the  width/ 

curve  is  no  less  than  500-1000  with/ of  the  interval  corresponding  to  the 
diameter  of  5  microns.  Apart  from  this  statement  --  the  insufficiency  of 
which  will  become  apparent  in  further  discussion  --  we  crould  rind  no  other- 
information  given  by  Houghton  to  this  effect  not  only  in  terms  of  quantitat¬ 
ive,  but  also  in  terms  of  qualitative  characteristics.  Yet,  the  need  for  a 
substantiated  estimate  of  the  required  sample  volume  is  evident.  No  appreci¬ 
able  advantages  from  the  standpoint  of  accuracy  can  be  derived  in  the  case 
of  larger  and  subsequently  growing  sample  volumes,  whereas  the  laboriousness 
of  measurements  increases  unproductively .  We  shall  define  the  sample  volume 
sufficiency  test  as  a  quantitative  ratio  permitting  to  determine  that  total 
number  of  particles  N^.  of  the  investigated  set  which  is  sufficient  for  the 
experimental  distribution  curve  to  be  characterized  at  the  maximum  by  the 
relative  error  (due  to  the  flniteness  of  the  sample  volume)  whose  magnitude 
does  not  exceed  the  preassigned  value  of  Dm»  Here  we  shall  disregard  the 
error  due  to  the  finiteness  of  the  interval  width  (considered  in  Section  2) 
owing  to  its  negligible  value.  This  is  admissible  with  respect  to  distri¬ 
bution  (29)  when  0*2  to  0.25*  To  construct  the  test  —  the  essence  of 

which  transpires  from  the  above  definition  —  we  shall  make  use  of  the 
height  of  curve  (16 )  in  rewriting  the  equation  in  the  form  providing  a 
solution  with  respect  to  nm ; 

nm  =  ™k^P* 

In  substituting  this  expression  for  nm  in  (34),  and  solving  it  with 


k» 


**  = 


xi 

P  ' 


(35a ) 


respect  to  N 


we  obt  ain 


whence,  upon  transition  to  the  dimensional  aspect,  it  follows  that 


JV*=-=^2- 

HDlXr 


(36) 


Jbrmula  ( 36 )  then,  con stitute s  the  criterion  for  sample  volume 

sufficiency  suited  for  any  unimodal  distribution.  In  order  to  estimate 

the  required  number  of  particles  N^.  it  is  necessary  to  preselect  the 

value  of  measurement  reliability^  sufficiently  close  to  unity,  and  to 

at  the  maximum  Dm* 

assign^  beforehand  the  desired  relative  error  value 

The  magnitude  of  this  error  should  be  negligible  as  compared 
to  unity.  Moreover,  it  is  important  to  know  the  tentative  value  of  the 
most  probable  radius  rQ  to  assure  the  most  effective  selection  ( s6e  Table 
2)  of  the  interval  width  Ar  =  ^pr0  .  The  peakedness  coefficient  in  forrnul; 
(36)  may  be  taken  as  equal  to  its  minimum  value  as  encountered  in  pract¬ 
ice  . 


Formula  (36)  for  H  -  0  *5^1 ,  corre  sponding  to  distributions 

with  the  most  obtuse  peak  (m  ■  2,  n  -  1),  with  reliability  =.  0.95 
(x  =  1.96),  may  be  rewritten  as  follows: 

The  application  of  formula  (37)  for  purposes  of  evaluation  may 
assure  a  substantiated  --  and,  consequently,  also  an  effective  -- 
selection  of  the  sample  volume  for  measurements  of  particle  sizes  in 
natural  fog  and  clouds.  Further  on  in  the  text,  in  a  corre  sponding 
section,  the  basis  of  the  sufficiency  test  justification  will  be  evalu¬ 
ated  by  way  of  comaprison  with  the  experimental  data  from  formula  (3^) 


constituting  the  foundation  of  this  test. 

In  this  manner  it  will  be  seen  that  in  contradistinction  to  Houghton's 

i 

ideas,  the  requred  sample  volume  in  our  test  turns  out  to  be  dependent 
A 

through  parameter  x  on  the  measurement  reliability  on  the  relative 
error  at  the  maximum  Dffl,  on  the  shape  of  the  curve  in  terms  of  the 
peakedness  coefficient,  and,  finally,  on  the  interval- width  to  the 
most-probable-radius-value  rQ  ratio.  Houghton's  considerations  referred 


to  by  us  in  the  beginning  of  this  section  cannot  serve  as  a,  sufficiency 

h 

criterion.  Rather,  this^a  corrrect,  but  sufficiently  gross  rule 
derived  directly  from  ample  empirical  data,  in  which  the  dependence  on 
the  width  of  the  interval  and  the  magnitude  of  the  most  probable  size 
remained  unreflected. 

6 .  Statistical  Filling  and  Multimodality. 

Published  data  on  size  measurements  of  particles  in  foj,  and  clouds 
abound  in  experimental  distributions  the  multimodality  of  which  is 
sometimes  attributed  a  special  meaning  ensuing  from  the  widely  known 
hypothetical  constructions  proposed  by  Kohler  /32/,  /33/  relative  to  the 
multiple  relationship  of  cloud  and  fog  particle  masses. 

Let  us  consider  the  problem  concerning  the  authenticity  of  the 
multimodality  of  experimental  distributions  viewed  from  the  positions 
of  statistical  probability,  as  this  question  was  originally  put  forth 
by  us  A/  in  1950. 

The  complexity  of  the  problem  is  aggravated  by  the  fact  that  multi¬ 
modal  distributions  —  as  noted  by  V.I.  Romanovsky  /13/7-  as  a  rule,  are 
not  studied  in  mathematical  statistics. 

In  the  preceding  section  we  have  derived  an  expression  for  the 
relative  error  in  estimating  the  number  of  particles  filling  an  indivi¬ 
dual  interval.  This  error  comes  as  a  result  of  the  boundedness  of  the 
sample  volume.  This  fact  makes  it  possible  to 
provide  the  experimental  curve  with  two  supp¬ 
lementary  curves.  One  of  them  is  them  is  the 
envelope  of  the  maximum  possible  deviations 
of  the  experimental  curve  from  the  general 
population  cure;  the  other  is  the  envelope 
of  the  minimum  possible  deviations  (fig  .5)* 

The  distance  between  the  above  envelopes  along 

Fig .  5  .  The  envelopes  of 

the  vertical  in  any  interval  is  evidently  errors. 


3/ 


(37a) 


equal  to  the  approximate  double  value  of  the  absolute  error,  i  •  e . 

2A  ni  s  2x  Y ni  ■ 

The  numerical  value  of  the  area  enclosed  between  the  error  envelopes 
for  a  given  distribution  essentially  depend  on  the  sample  volume,  the 
measurement  reliability,  and  the  width  of  the  interval,  as  this  follows 
directly  from  (34).  Reliability  0<.  sufficiently  close  to  unity,  as  this 
is  usually  assumed  when  dealing  with  practical  problems  /l4/  ,  represents 
in  the  case  under  consideration  the  probability  with  which  one  may 
expect  the  appearance  of  an  arbitrary  curve  contour  in  the  zone  bounded 
by  the  error  envelopes.  Hence,  this  zone  may  be  considered  as  a  locus 
where  each  point  may  belong  to  any  arbitrarily  selected  distribution 
curve  provided  the  normalization  condition  for  it  is  observed.  Each 
normalized  distribution  curve  traced  through  this  zone  may  be  a  populat¬ 
ion  curve.  The  probability  for  the  contour  of  a  normalized  distribution 
curve  to  pass  beyond  the  error  envelopes  for  each  individual  Interval, 
apparently,  is  equal  to  1  -  ok*  The  probability  for  the  contour  to 
emerge  beyond  the  lower  envelope  is  smaller  even  the  more  so*  In  accord¬ 
ance  with  A*K.  Mitropol 1  sky 1  s  suggestions  /l2/  we  shall  call  the  zone 
enclosed  between  the  error  envelopes  as  the  zone  of  fluctuations. 

The  portion  of  the  curve  area  bounded  by  the  lower  error  envelope 
and  the  X-axis  is  an  authentic  part  of  the  experimental  curve  area, 
since  the  penetration  of  the  contour  of  any  experimental  curve  from  the 
fluctuation  zone  into  this  region  is  so  improbable  that  it  can  be 
considered  as  practically  impossible. 

A  very  useful  characteristic  of  the  autheticity  zone  is  the  fact 
that  with  infinitely  growing  sample  volume,  and  infinitesimal  intervals, 
it  has  the  whole  area  of  the  population  curve  as  its  limit.  The  zone 
between  the  error  envelopes  (fluctuation  zone)  ultimately  vanishes. 

This  would  not  have  been  possible  had  there  been  no  limit  for  each 
authentic  part  of  any  individual  interval,  and  had  the  relative  errors 

3^ 


not  been  tending  to  zero 


So  far  we  were  dealing  mainly  with  experimental  curves  plotted  on 


the  basis  of  measurements  involving  a  finite  and  constant  number  of 
particles . 

There  arise  a  number  of  questions  in  practical  statistical  measure¬ 
ments  the  solution  of  which  is  possible  only  with  due  consideration  of 
the  variability  of  the  sample  volume.  For  example:  What  is  the  minimum 
number  of  particles  required  to  assure  that  a  curv6  plotted  as  a  result 
of  their  measurement  could  be  considered  as  an  experimental  distribution  ? 
Evidently,  one  particle  is  clearly  insufficient  for  this  purpose.  .An 
accounting  of  the  sample  volume  variability  is  equivalent  to  a  test  in 
which  the  volume  of  a  sample  taken  from  a  fixed  set  of  particles  grows 
as  a  result  of  sequential  incorporation  into  it  of  individual  particles. 

In  such  an  experiment  we  shall  be  interested  mainly  in  the  changes  mani¬ 
fested  by  the  experimental  curve  as  a  result  of  statistical  filling  which 
is  understood  to  mean  a  process  connected  only  with  an  increase  in  the 
sample  volume . 

In  passing  now  on  to  an  aritrary  and  separately  taken  interval  and 
assuming  the  maximum  value  of  the  relative  frequency  density  for  it  to 
be  equal  to  the  probability  density  n(^)  of  the  general  population,  we 
shall  obtain 


(38) 

as  the  mean  value  for  the  authentic  portion  of  the  interval  particle 
concentration,  and  correspondingly  2n((j))D^  for  the  components  of  the 
fluctuation  zone. 


Evidently,  the  difference  1-Di  in  (38) 
characterizes  the  authetic  part  of  probability 
density.  In  the  case  of  an  experimental  distrib¬ 
ution  it  characterizes  the  authentic  part  of  the 

ni. 

relative  frequency  density  — -=■  for  a  given 

/SpNk 

Fig .6  Statistical  interval.  Hence,  this  difference  may  be  consider- 

filling  of  an  interval. 

ed  as  a  measure  of  the  statistical  filling  of  the  curve  in  the  invest!- 


gated  interval  •  Figure  6  shows  a  graph|(  of  function 

1  Dt  /(«,.)=!  -  y-  ,  (35a) 

from  which  it  transpire^  that  --  depending  on  the  relationship  between 

nffl  and  Dm"*  there  are  for  a  given  reliability  value  three  basic  stages  of 

K 

statistical  filling  or  experimental  curve  shaping.  We  shall  now  pass  on 
to  a  discussion  of  these  stages* 

The  Initial  (i)  Shaping  Stage. 

The  most  densely  filled  Interval  contains  nm  <  x2  number  of  particles 

and  has  no  authentic  portion  in  its  filling,  i.e.  cases  are  possible 

when  Anm  >  nm-  Fbr  all  the  remaining  intervals  the  inequalities 

will  be  even  more  sharply  expressed.  The  experimental  curve  in  the 

first  stage  of  shaping  is  characterized  by  absolute  Indefiniteness  due 

to  the  absence  of  a  positive  authenticity  zone  portion.  Such  a  curve,  as 

a  rule,  has  nothing  in  common  with  the  di  stribution  of  the 

random 

universe  from  which  the/s  ample  was  drawn.  This  may  be  manifested  in  the 
fact  that  the  universe  distribution  mode  would  fail  to  coincide  with 
the  the  most  den^  sely  filled  Interval.  The  number  of  peaks  of  such  a 
curve  in  isolated  cases  may  be  comparable  with  the  sample  volume.  It 
may  not  resemble  the  initial  distribution  by  form.  Addition  of  two- 
three  particles  may  substantially  change  the  shape  of  such  a  curve, 
the  number  of  maxima,  and  the  position  of  the  most  densely  filled  interval. 

A  sample  with  a  volume  sufficient  for  an  experimental  curve  in  the 
first  shaping  stage  must  be  considered  as  totally  unrepre sent ative  from 
statistical  point  of  view,  provided,  of  course,  that  there  had  been  no 
subdivision  into  excessively  small  intervals. 

The  Main  (II)  Shaping  Stage. 

The  most  densely  filled  interval  contains  a  number  of  particles 

x2 


within  the  limits 


where  Dm  ^  1  •  Here  the  case  corresponds  to  the  beginning  of  the 


m 


appearance  of  the  authentic  part  of  the  interval  filling,  and  n 


m 


signifies  the  filling  with  error  Dm  at  the  maximum,  when  the  authentic 
portion  of  the  interval  filling  grows  to  (1-Dm)nm.  By  the  end  of  the 
second  shaping  stage  the  experimental  curves  in  their  main  central  part 
have  a  shape  corre sponding  gflV  to  the  population  distribution.  The  most 
densely  filled  interval,  if  it  is  not  too  small,  comprizes  the  mode  of 
the  initial  distribution.  The  multimodality  registered  in  the  central 
part  of  the  experimental  curve  at  the  beginning  of  the  second  shaping 
stage  shifts  by  the  end  of  the  sat^g  into  the  region  of  the  curve  tails. 
Further  addition  of  no  matter  how  large  a  number  of  particles  does  not 
substantially  alter  the  empiric  curve  shape  at  the  end  of  the  second  stage. 

2h.  experimental  curve  constructed  up  to  the  end  of  the  second  shap¬ 
ing  stage  graphically  represents  the  sample  which  is  statistically 
representative,  provided  the  error  Dm  with  which  the  peak  was  determined 
satisfies  the  requirements  of  a  given  problem. 

The  Ultimate  (ill)  Shaping  Stage. 


The  curve  maximum  is  filled  to  the  extent  that 

A-2  ,  . 

~02  <  «„  <  CO. 


(A) 


Moreover,  the  inequality  Dm  continues  to  grow  and  ultimately 

Dm  -  0.  At  this  shaping  stage  the  experimental  curve,  without  exhibit¬ 
ing  any  appreciable  changes,  attains  the  state  of  perfect  Identity  with 
the  universe  distribution.  The  numeric  al^  value  of  any  experimental 

curve  parameter  approaches  asymptotically  to  its  maximum  value.  With 

observed 

regard  to  the  curve  altitude  this  fact  may  be  CHW  in  Fig •  3 • 

Unlike  the  first  two  stages  of  experimental  distribution  curve 
shapping,  the  last  stage  is  characterized  by  the  fact  that  the  size  and 
the  shape  of  the  experimental  distribution  curve  authenticity  zone  turn 
out  to  be  practically  independent  from  the  volume  of  the  random  sample. 
In  the  light  of  the  above  ideas,  only  such  experimental  curves 


as  have  been  plotted  up  to  the  end  of  the  second  and  the  beginning  of 

AS 

the  third  shaping  stagey  (1  .e  *Aare  statistically  filled),  appear  to  be 
suitable  for®fcan  analysis  of  possible  multimodality,  its  character, 
and  causes  of  its  origin. 

Consideration  of  the  character  of  statistical  filling  makes  it 
possible  to  develop  a  reality  test  for  the  multimodality  of  emperical 
di stribution  s • 


Authentic  or  real  multimodality  of  an  arbitrary  experimental  dis¬ 


tribution  is  understood  to  mean 


a  multimodality  only  such  as 


acquires  a  perfectly  definite  character  in  terms  of  number,  shape,  and 
relative  position  of  the  peaks  in  the  process  of  further  growth 
of  the  random  sample  volume. 


If,  however,  as  a  result  of  observations  a  multipeak  experimental 
curve  is  produced  which  with  further  growth  of  random  sample  volume 
turns  into  a  single-peak  curve  only  in  consequence  of  statistical  fill¬ 
ing,  then  the  multimodality  of  an  experimental  distribution  is  not 
authentic  and  we  shall  consider  it  as  accidental  multimodality  result¬ 
ing  from  the  boundedness  of  the  sample  volume,  or  —  which  is  the  same 
--  from  insufficient  statistical  filling  of  the  experimental  distribut¬ 
ion  intervals  which  interest  us. 


Thus,  the  best  method  for  the  solution  of  the  question  relative  to 
the  authenticity  of  the  multimodality  of  experimental  distribut¬ 

ions  i  s  to  incre  ase  the  volume  of  the  random  sample  to  assure  the 


required  statistical  filling. 

When  random  sample  volume  cannot  be  increased,  the  character  of 
the  observed  multimodality  may  be  approximately  estimated  in  the  follow¬ 
ing  manner:  the  contour  of  the  investigated  empirical  curve  is  accompan¬ 
ied  by  the  envelopes  of  errprs  calculated  from  formula  (3^)»  If  it 

to  through  the  tone 

proves  to  be  impo  ssibl e/draw  a  single  one-peak  curve 


enclosed  between  the  two  envelopes  without  the  contour  of  the  curve 


penetrating  outside  thtS  zone  (see  Fig. 5),  then  the  multimodality  of 
such  an  experimental  distribution  can  be  considered  as  practically 
authentic.  One  might  proceed  also  in  the  inverse  order:  by  superposing 
a  normalized,  averaged  single-peak  curve  with  error  envelopes  on  the 
investigated  experimental  distribution.  multimodality  the  ampli¬ 

tudes  of  which  do  not  noticeably  pass  beyond  the  error  envelopes  ashatild 
be  considered  as  accidental.  One  should  bear  in  mind  here,  that  the 

very  concept  of  reliability  presuposes  the 
possibility  of  a  transgression  ofQH  accident¬ 
al  multimodality  amplitudes  in  1- (leases,  if 
the  number  of  the  investigated  experimental 
curves  is  sufficiently  large. 

The  necessary  corollary  an  analysis 

experiment al  distribution  multimodality  con¬ 
sists  in  that  not  any  multimodality  is  authen¬ 
tic  or  real.  Consequently,  ifcst  physical  in- 
the  causes  for 

vestigatlon  of  the  multi- 

modality  of  experimental  curves  must  be  pre¬ 
ceded  by  its  verification  from  the  statistical  point  of  view.  The 
method  of  the  envelopes  is  perfectly  convenient  for  this  purpose.  On 
the  basis  of  the  above  ideas  it  became  possible  to  review  the  extensive 
material  on  experimental  distributions  published  in  the  works  by  HaSen>ann 
/30/,  Chestnaya  /22/,  Zaytsev  /l/.  Diem  /2rj/,  /28/,  and  Borovikov  /l/. 

These  papers  distinguish  themselves  advantageously  from  other  works  by 

a 

the  considerable  number  of  observed  case  s, /detailed  description  and 
sufficiently  comprehensive  Information  on  the  investigated  distributions. 
The  papers  of  other  authors  are  not  always  as  interesting.  Flndeisen  /29/, 
for  example,  has  produced  distributions  many  of  which  have  no  maximum. 

A LI  his  distributions  are  characterized  by  a  considerable  shift  of  the 
curve  to  the  left  from  the  values  obtained  by  other  investigators.  This 


ing  of  experimental  dis¬ 
tributions  (319  curves). 

1  -  after  Diem,  2  -  after 
Zaytsev,  3  -  after  Borovi¬ 
kov. 


is,  apparently,  due  to  the  excessive  heating  of  the  drops  in  the  measuring 

volume  of  his  apparatus.  The  results  of  many  years  of  Houghton’s  /3l/ 
research 

are  reduced  to  three  "typical"  curves  and  the  direct 
observation  data  are  not  published.  Bricard  /25/  also  avoided  to  publish 


his  data  on  distribution  curves.  ALII 


experimental  curves  considered 


by  us  were  produced  by  the  commonly  known  method  of  taking  photomicro¬ 
graphs  of  drops  collected  on  a  glass  slide  coated  by  a  thin  layer  of 

mineral  oil.  Information  on  this  method  and  its  variants  was 

the 

repeatedly  published  in  tfMifepapers  by  Tarayan  /l8/,  Katchenlcov  /8/, 

Fuks  /20/,  Hagemann  /30/,  Selezneva  /l5/,  Kucherov  /lo/,  Chestnaya  /22/ 
and  other  authors.  In  a  more  comprehensive  form  these  data  are  presented 
by  Zaytsev  /6/  and  Borovikov  /l/.  From  the  statistical  viewpoint  not 
all  the  experiment al  distributions  referred  to  above  are  equivalent. 

The  sample  volume  in  each  of  them  is  the  standard  which,  in  the  main, 
determirPs  the  quality  of  the  experimenatl  material.  Represented  in 
Fig.  7  are  the  integral  distribution  curves  produced  by  Diem  /27/,/28/, 
Zaitsev  /l/ ,  ana  Borovikov  /l/.  In  Table  3  one  will  find  the  comparative 
data  showing  the  minimum,  average,  and  maximum  sample  volumes,  as  well 
as  the  number  of  experimental  curves  classified  according  to  the 
stage  s. 

A  comparison  of  the  curves  from  Fig. 7  and  the  data  from  Table  3  shows 


that  the  best 


from  the  standpoint  of  statistical  filling  are 


Borovikov’s  experimental  distributions.  The  worst  are  the  data  given 
in  Diem's  two  papers.  The  results  of  Hagemann' s  /30/  and  Chestnaya' s 
/22/  observations  are  hardly,  distinguishable  from  the  distributions 
produced  by  Zaytsev  /l /  which  occupy  an  intermediate  position.  The 
principal  characteristic  of  Diem's  distributions  is  that  the 
of  most  of  his  curves  was  stopped  at  the  firstcstage.  The  majority  of 
his  curves  were  not  brought  up  to  the  required  level  of  statistical 
filling.  A  common  defect  of  all  the  published  distributions  consists 


in  the  fact  that  sample  volume  of  experimental  curves  was  taken  at  random. 

A 

It  goes  without  saying  that  from  the  experimental  point  of  view  it  is 
easier  to  catch  many  small  drops  and  much  moiE  difficult  to  collect  a 
large  number  of  large  drops.  Experimental  distributions  with  small 
sample  volumes  turned  out  to  be  constructed,  as  a  rule,  for  the  cases 
of  large-drop  clouds,  whereas  large  sample  volumes  were  used  in  the  main 
for  the  cases  of  fine-drop  clouds.  Yet,  it  follows  from  the  sufficiency 
test  xhs  for  the  case  when  the  interval  width  is  constant  that  one  should 
take  the  experimentally  more  difficult  path  in  order  to  produce  statistic¬ 
ally  comparable  experimental  data. 

T  a  b  1  e _ 3. 


Author 

Sample  volume 

Number  of  curves 

mini¬ 

mum 

ave¬ 

rage 

maxi¬ 

mum 

Nk  <  50 

Nk  =  50 

-500 

Nk  >  500 

Total 

Borovikov 

/V  , 

Zaytsev  h/ 
Diem  /27/, 
/28/ 

63 

8 

10 

500 

250 

60 

2'155 

1388 

640 

5 

40 

'\ 

44 

85 

92 

44 

20 

2 

<ss 

no 

1.34 

V le  shall  now  pass  to  the  problem  of  the  multimodality  of  experiment¬ 
al  curves  and  its  connection  with  the  well  known  relationship  of  multiple 
masses  which  was  originally  registered  by  Defant  /26/  with  respect  to 
rain  drops,  and  later  extended  by  Kohler  over  the  particles  of  natural 
fogs  and  clouds. 

The  method  of  error  envelopes  is  most  convenient  for  the  investigat¬ 
ion  of  the  amplitudes  of  multimodal  experimental  distributions.  The 
fluctuation  zone  here  should  be  considered  as  the  region  of  possible 
accidental  multimodality. 

with 

Analysis  of  the  published  experimental  distributions  high  reliabil¬ 
ity  D^shows  that  all  of  them  have  a  multimodality  of  accidental  character. 
The  number  of  cases  involving  the  passage  of  experimental  curves  in 


expressed  in  percentage 


individual  intervals  beyond  the  error  envelopes  , 
values,  is  noticeably  lower  than  the  admissible  5%  ( <*  =  95)* 

In  analyzing  the  possible  errors  of  the  photomicrography  method, 

A»M.  Borovikov  /l/  comes  to  the  conclusion:  "in  our  opinion  the  total 
error  in  measuring  the  number  of  drops  does  not  exceed  10 %,  and  then 
this  error  in  the  overwhelming  majority  of  cases  affects  the  fine-drop 
measurements"  (page  19  )•  This  conclusion  may  be  assumed  as  correct  only 
for  the  intervals  statistically  so  filled  that  error  Dj_  in  them  is 
negligible  in  comparison  with  other  statistical  errors.  For  many  experi¬ 
mental  distributions  analyzed  by  us  error  exceeds  considerably  all 
other  errors  put  together. 

A  structural  analysis  of  multimodality  in  experimental  distribut¬ 
ions  to  determine  its  correspondence  to  the  relationship  of  multiplicity 
of  masses  shows  that  this  relationship  does  not  hold  for  the  entire 
length  of  any  experimental  curve.  This  fact  was  already  registered  in 
Ye»S.  Selezneva's  /l6/  and  1. 1.  Che  stn  ay  3^/2  2/  works.  Of  Interest  is 
the  question  regarding  the  conditions  which  should  be  satisfied  by  a 
statistically  unfilled  experimental  distribution  to  cause  its  multi¬ 
modality  to  satisfy  the  multiplicity  ratio.  For  this  purpose  we  shall 
examine  the  maxima  r and  r2  of  an  arbitrary  multipeak  experimental  curve. 
These  maxima  satisfy  the  multplicity  relation  in  the  case  of  the  obvious 


ratio 


• — y~2  • 


(39) 


The  distance  between  maxima  r-^  and  r g  expressed  through  one  of  them 
may  be  found  from  (39)  as 

r*~  ri=(v’2-  l)r,.  (39a) 


The  same  distance  may  also  be  expressed  in  terms  of  the  Interval 


width  ^r  of  a  given  experimental  distribution  in  the  form  of 

r2  —  rx  =  (l+  1)A  r. 


(40) 


where  is  the  number  of  intermediate  intervals  situated  between  maxima 

rl  ^  ^  A, 


3 _  strictly/ 

.It  follows  from  (39)  and  (4o)  that  maxima  r^  and  r 2  “  V  2r^ /satisfy 
the  mass  multiplicity  ratio  if  condition 

(41) 

V‘2  —  \ 

is  observed. 

Observations  indicate  that  the  sequence  with  which  the  maxima 
alternate  is  fortuitlous.  Furthermore,  most  frequent  are  the  cases 
when  the  maxima  alternate  after  one  interval,  which  corresponds  to  the 
case  d-  1.  Recorded  below  by  way  of  an  example  is  a  table  of  the 
recurrence  of  6.  •  The  data  for  this  table  are  taken  from  the  experimental 
material  produced  by  1. 1.  Chestnaya  /22/ 

Table  4. 


8 

9 

10  and 
more 

Total 

Number  of  |  |  |  w  |  "  |  |  ”  |  ”  |  ‘  “  |  |  3^6 

cases 

For  Chestnaya' s  distributions  constructed  with  intervals  of  1*3 
cf- 

microns  (in  di  ameter )  in  accordance  with  (4l)  one  should  expect  that 
the  mass  multiplicity  relation  may  best  be  satisfied  at  (  .  1  in  the 
size  group  10-12.6  microns,  and  for  i  =  2  in  the  size  group  15-19  microns 
respectively. 

of  87  experimental  curves  of  Chestanya  there  was  only  one  for  which 
the  multiplicity  relation  was  met  over  a  large  segment  of  the  curve. 

We  shall  cite  this  example  in  designating  the  position  of  the  six  peaks 
of  the  curve  by  figures  in  parentheses  to  indicate  the  number  of  inter¬ 
mediate  intervals  between  the  maxima:  7*8(1),  10.4(1),  13»o(l),  15.6(2), 
19*5(2),  23.4  microns. 

This  example  graphically  demonstrates  the  mechanism  of  producing 
experimental  distributions  satisfying  the  multiplicity  relationship. 

It  can  be  shown  without  to  much  of  an  effort  that  for  a  limited  sample 
from  any  set  of  real  particles  (including,  of  course,  a  single -peak  one) 


the  compliance  with  the  mass  multiplicity  relation  may  be  demonstrated 
on  individual  experimental  curves  in  any  preselected  region  of  the 
experimental  curve  and  for  any  maxima  known,  beforehand.  For  this 
purpose  —  apart  from  the  statistical  incompleteness  without  which  it  is 
possible  that  there  might  be  no  accidental  multimodal ity  --  it  is  only 
necessary  to  select  the  interval  width  in  conformity  with  condition  (4l). 

It  is  well  known  that  one  of  the  elementary  requirements  in  mathemat¬ 
ical  statistics  is  that  the  samples  should  be  drawn  randomly.  For  an 
experimental  confirmation  of  the  mass  multiplicity  relation  it  turns  out 
to  be  necessary,  in  violation  of  this  rule,  to  resortjto  a  selective 

separation  of  a  certain  portion  of  the  experimental  data,  and  then 

Kohler 1  s 

precisely  that  part  of  it  which  is  useful  for  the  substantiation  of  Plif 
hypothetical  views.  This  is  the  second  serious  error  of  statistical 
character  made  by  KiShler  in  his  analysis  of  the  structural  characterist¬ 
ics  of  the  multimodality  of  experimental  distributions. 

al  so 

Let  it  be  mentioned  in  passing  that  the  method  used  by  KBhler/offers 


an  equal  opportunity  to  "substantiate"  a  counter  hypothesis  concerning 
the  predominant  non-coalescence  of  equal-size  drops.  To  achieve  this 
goal  it  suffices  to  select  those  cases  from  the  observation  results  in 
which  £he  maxima  of  a  statistically  unfilled  distribution  curve  are 
followed  by  minima  appearing  at  definite  distances. 

Thus,  as  we  can  see,  the  "Kdfhler  law", 

"observed"  mass  multiple 


as  the  allegedly  experimentally 


relation  is  sometimes  referred  to)  has  no  direct 
bearing  on  the  microphysics  of  clouds  and  fog  and  is  not  a  physical  law. 

As  a  matter  of  fact,  the  phenomenon  registered  by  Kohler  is  the  inevitable 
result  of  the  manifested  boundedness  of  the  sample  volume. 

7  •  The  Mean  Sizes  in  Poly  disperse  Fog. 

In  distribution  research  in  natural  fog  and  clouds  it  it  is  custom¬ 


ary  to  make  use  of  the  the  concept  of  the,  so-called,  mean  sizes.  The 
following  designations  may  be  found  in  literature:  the  most  probable, 


mean 

the  predominant,  the  arithmetic  mean,  the/quadratic ,  the  mean  cubic, 

© 

the  mean  size  after  Houghton,  the  average  size.  Recently  N.S.  Shishkin 

radius  of  the  drops 

(23)  has  utilized  the  concept  of  the  responsible  for 

the  greatest  contribution  to  the  liquid-water  content-  By  analogy  one 

might  as  well  speak  of  the  radius  of  the  drops  responsible  for  the 

greatest  contribution  to  the  surface  area-  As  mentioned  before,  different 

methods  of  experimental  distribution  curve  construction  are  adopted  in 

the  Soviet  Union  and  Europe,  on  the  one  hand,  and  in  the  USA,  on  the 

other  hand.  As  a  result  of  this  discrepancy  in  methods  there  is  no 

uniformly  accepted  system  of  mean  size  concepts  and  designations,  and 

this  sometime  leads  to  obvious  errors.  Ye-S-  Selezneva  (16 )  has  pointed 

out  in  her  paper  to  one  of  such  misconceptions  of  M-  Diem  who  believed 

that  the  American  fogs  consist  of  larger  drops  than  the  European-  One 

could  mention  another  absurdity  of  this  kind  encountered  in  J.  Bricard's 

work  (23 ),  where  the  arithmetic  mean  size  was  referred  to  as  simply  an 

"average"  and  was  used  for  the  calculation  of  liquid-water  content  on 

the  basis  of  a  certain  number  of  observed  particles  in  a  unit  volume. 

It  is  evident  that  in  this  case  it  was  necessary  to  use  the  mean  cubic 
inaccuracy 

size.  The  in  liquid-water  content  determination  arising  from  this 

type  of  error  depends  on  the  nature  of  the  distribution  and  in  indivi¬ 
dual  cases  may  be  comen  surate  with  the  evaluated  liquid  water.  Moreover, 
it  should  be  borne  in  mind  that  the  concept  "mean  size"  without  a 
reference  to  the  meaning  of  the  mean  is  hollow  and,  consequently,  errone¬ 
ous.  Owing  to  the  above  mentioned  difference  in  the  methods  of  experi¬ 
mental  curve  construction  one  and  the  same  mean  size  is  sometimes 
designated  in  a  different  way,  while  different  sizes  are  often  identical¬ 
ly  designated.  Thus,  for  example,  Shishkin  /23/  refers  to  the  mode  of 
volume  distribution  as  the  radius  of  the  drops  responsible  for  the 
greatest  contribution  to  liquid-water  contents.  In  Houghton's  paper 
this  is  called  the  predominant  size.  The  predominant  size  in  Houghton's 


+3 


interpretation  and  the  predominant  size  in  works  by  Soviet  investigators 

A 

/l,7,l6/  are  not  one  and  the  same  thing  because  of  the  different  method 
of  e^mprical  curve  construction.  Houghton  /3l/  uses  only  one  "average 
size"  which  he  has  calculated  by  formula 

r  "  S nrt  (  A) 

without  distinguishing  it  from  his  predominant  size^  in  assuming  that 
they  are  coincident.  In  the  section  devoted  to  the  method  of  coronal 
rings,  we  read:  "in  usual  application  it  yields  only  an  average  or 
predominant  particle  size"  (page  5)*  The  only  "average  size"  of 
Houghton  does  not  coincide  with  the  only  "average  size"  of  Bricard  /25/, 
who  has  computed  it  as  the  aritmetic  mean  siZ6,  and  utilized  it  as  the 
mean  cubic  size.  Finally,  one  also  finds  such  studies  /9/,  /l 6/  in  which 
the  mean  cubic  size  is  identified  as  the  mode  of  volume  distribution^ 
with  a  reference  to  which  we  have  begun  this  discussion. 

The  a  monodlsperse 

fog  such  as  would  in  all  cases  be  equivalent  to  the  investigated  poly- 
disperse  fog  is  fundamentally  Impossible.  In  view  of  this  fact  we  are 
compelled  instead  of  a  polydisperse  fog  to  construct  for  each  separate 
case  a  model  of  such  a  monodlsperse  fog  that  would  be  equivalent  to 
polydisperse  fog  only  with  respect  to  a  given  property  called  the 
determining  property  in  statistics. 

The  calculation  operation  for  one  or  another  mean  size  does  not 
determine  only  the  particle  size  of  a  certain  equivalent  monodlsperse 
fog,  but  simultaneously  also  establishes  a  corresponding  equivalence 
condition  with  respect  to  the  given  determining  property. 

To  clarify  the  question  oif  the  correspondence  of  the  exposed 
ideas  on  the  mean  particle  sizes  of  a  polydisperse  fog  to  the  basic 
positions  of  mathematical  statistics  in  this  matter  we  shall  turn  to 
a  pertinent  passage  in  V.I.  Komarovsky' s  book  / 13/:  "The  most  general 
and  precise  definition  of  the  concept  of  the  mean  was  formulated  almost 


simultaneously  (in  1929)  by  the  Italian  scientist  Cazini  and  the  Soviet 
scientist  Boyarskiy.  We  shall  now  quote  this  definition  which  will  be 
called  the  Boy arskiy-Cazini  definition. 

Supposing  we  have  a  certain  statistical  collective  S  with  volume  n 
which  is  characterized  by  values 

xlf  x2f  . . . ,  %  (28.1) 

of  a  quant ative  variable  x.  As  the  determining  property  of  collective 
A 

S  we  shall  designate  such  its  property  related  to  values  (28.1)  which 
must  remain  invariable  at  all  possible  variations  of  values  (28.1). 

Then  the  variable  mean  x  with  respect  tp  the  considered  determin¬ 
ing  property  of  collective  S  will  be  called  that  value  of  variable  x, 
equal  for  all  the  collective  members,  which  may  be  attributed  to  them 
without  changing  the  determining  property  of  the  collective",  (page  95)* 
From  the  great  abundance  of  various  means  which  are  known  in 
mathematical  statistics  we  are  Interested  only  in  those  which  belong  to 
the  characteristic  class  of  means  called  regular.  (According  to  Romaonv- 
skiy ,  included  among  these  means  are  such  which  are  monotone,  associat¬ 
ive,  and  single -valued,  while  the  determining  property  is  given  In  the 
form  of  a  continuous  and  symmetric  function). 

About  such  a  me  an  and  its  correct  statistical  interpretation  Roman- 
ovskly  /13/  writes;  "  This  quantity  is  obtained  when  under  the  condition 
of  invariability  of  the  determining  property,  instead  of  the  given 
collective,  and  abstract  levelled  collective  is  constructed  having  one 
and  the  same  criterion  x  for  all  Its  members.  This  quantity  in  the 
sense  of  preserving  the  determining  property  is  thus  equivalent  to  the 
totality  of  x^  values  in  the  investigated  collective  and,  therefore, 
actually  represents  one  of  the  means  of  describing  the  statistical 
collective  which  characterizes  the  collective  with  respect  to  the 
determining  property  as  completely  as  the  totality  of  x^  values  of 
variable  x"  (page  101). 


The  objective  of  this  paragraph  i  s  to  venture  to  reconsider  once 


again  the  concepts  of  the  basic  mean  particle  sizes  of  polydisperse 

fog  together  with  their  de  signatiox®  in  accordance  with  the  ideas  outlined 
above  <> 


The  Most  Probable  Radius. 


As  stated  earlier  we  have  adopted  the  most  probable  radius  as  the 
initial,  the  unit  radius.  It  becomes  clear  from  the  designation  itself 
that  this  size  corresponds  to  the  radius  of  drops  such  as  are  being 
most  frequently  encountered  in  size  distributions. 

For  distribution  (7)  given  in  analytical  form  the  most  probable 
radius  r  =  rQ ,  or  p  “1#  corresponds  to  the  maximum  of  probability 
density.  Therefore,  its  magnitude  may  be  determined  by  solving  equation 


(42) 


with  respect  to  r. 

In  the  case  of  an  emperic  distribution  based  on  size  measurements 
of  a  finite  number  of  particles,  their  sizes  are  usually  plotted  on  the 


X-axis,  On  the  Y-axis  one  plots  the  relative  frequency 


nl 

X* 


who  se 


density,  if  the  number  of  particles  -  £n-j_  and  small  intervals  Ap  is 
sufficient,  will  be  qqual  to 

tli  .  s. 


Ap  Znt 


the 

where  n(p)  is  the  function  of  probability  density  of 

distribution  (7)  which  in  this  case  is  assumed  to  have  a  range  satisfying 
the  experiment. 

From  the  approximate  character  of  equation  (43),  which  is  an 

with 

expression  of  the  law  of  large  numbers,  it  follows  that  df  n((j))  unknown 
to  us  the  relative  frequency  density  maximum  appears  to  coincide  with 
probability  density  only  for  an  infinite  sample  volume  and  infinitesimal 
interval  width.  Therefore,  the  maximum  of  any  experimental  curvey  and 
consequently,  also  the  most  probable  radius  for  any  experimental  distri¬ 
bution  is  always  determined  approximately,  and  the  more  accurately,  the 


greater  the  sampling  scale  and  the  smaller  interval s.  The  remaining 

Q 

mean  sizes  which  will  be  considered  later  will  be  measured  in  fractions 
of  the  most  probable  radius. 

Apprt  for  size  distributions,  other  distributions  are  also  possible, 
as  noted  in  section  1.  Each  of  these  distributions  has  its  mode  disting¬ 
uishable  from  the  most  probable  size.  The  niode  of  surface-area  distribut¬ 
ion  s  ^  s,  expre  ssed  in  linear  measures,  is  the  predominant  size  in  terms 
of  surface  area,  whereas  the  mode  of  volume  distributions  ^>v  represents 
a  predominant  size  in  terms  of  volume,  and  is  also  expressed  linearly. 

In  future  di scussion  we  shall  refer  to  them  as  such. 

Houghton  /3l/  in  utilizing  a  volume  distribution^  erroneously  refer¬ 
red  to  the  mode  <af  this  distribution  as  the  "predominant"  size.  Lacking 
in  this  designation  is  the  dstx  indication  of  the  criterion  in  terms  of 
which  the  predominance  was  registered.  The  size  distribution  mode  or 
the  most  probable  size  sometimes  /l/,  /l / ,  /l6/ai?e  referred  to  as  the 
predominant  size.  Statistically  this  is  correct,  but  if  one  takes  into 
account  the  presence  of  the  predominant  size  in  Houghton's  investigation, 
then  it  becomes  apparent  that  the  use  of  this  term  is  undesirable  and 
should  be  discontinued  mainly  to  avoid  confusion  similar  to  that 
described  above. 

Linear  Equivalence  Radius,  or  Arithmetic  Mean  Radius. 


Let  us  select  linear  particle- size  s  for  an  investigated  set  to  serve 

as  the  determining  property.  For  this  purpose  Itf  suffices  to  equate  the 

total  length  formed  by  polydisperse  fog  drops  /"j/  to  the  same  number 

of  particles  of  a  monodisperse  fog.  The  condition  of  linear  equivalence 

may  be  expressed  as  ^ 

2A !kf  11  (p)  pdf.  =  2jV*p7  ,  (44) 


where  p is  the  particle  radius  in  monodisperse  fogy  equivalent  to 

given  polydisperse  system  with  respect  to  linear  sizes.  From  the 

expression  _  — 

Pi  =  —  \  n  (p)  prfp. 

0 
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(45) 


ensuing  from  (44)  and  valid  only  for  distributions  given  in  a  different¬ 
ial  form,  it  is  easy  to  pass  over  to  a  form  suitable  for  experimental 
distribution  investigations.  For  this  purpose  in  formula  (45)  the 


probability  density  should  be  replaced  by  its  approximate  value  from 
ni 

(43),  viz.  = — —  which  is  the  relative  frequency  density,  and  in 
accordance  with  (6)  to  assume  p  and  dp  equal  to 


d',  ~  A»  =  . 

h)  '  '  r0 


(45a) 


In  'replacing  the  integral  sign  by  the  summation  sign,  we  have 


r\  hi 


i  -e . 


—  _  Zn,r, 
r‘  ~  ■ 


(45b) 


(46) 


Leaving  aside  the  question  of  the  accuracy  of  the  approximate 

c.  o  p 

>  •  ^  t  i.  \  .  ,  j  * _  —  »  j  > 


relationship  (46)  as  a  matter  transcending  the 


•of  this  dis¬ 


cussion,  we  shall  only  state  that 


with  tending  to  infinity  and 


vani  shing  Ar ,  the  approximate  expression  has  as  its  stochastic  limit 

the  exact  expression  (45),  of  course,  if  the  finite-volume  random 

sample  is  drawn  from  distribution  (7 )  •  The  aritmetic  mean  size 

estimated  on  the  basis  of  formula  (46)  is  determined  with  an  error, 

the  greater,  the  smaller  the  volume  of  the  random  sample. 

-  ni 

The  arithmetic  mean  radius  p^  =  jr—  -g'n. e  radius  of  the  drops  of  a 

monodisperse  fog  such  as  for  the  same  number  of  particles  possesses  a 

sum  of  linear  sizes  identical  to  that  of  a  poljidi sperse  fog.  In  operat¬ 
ing  with  the  concept  of  aritmetic  mean  radius,  we^  by  virtue  of  this 
very  f  act  j  are  actually  dealing  with  a  monodisperse  fog  equivalent  to 
the  investigated  polydisperse  fog  in  terms  of  linear  sizes.  The  basic 
usefulness  of  mean  size  consists  in  the  fact  that  it  (after  a  simple 
calculation)  approximates  most  closely  the  distribution  mode  which  for 
skewed  experimental  distributions  is  usually  not  calculated.  The 
designation  "mean  size"  in  application  to  Fj_  used  in  a  number  of  papers 
/l/,  /25/f  and  so  on,  is  undesirable  in  view  of  the  indefiniteness  of 
this  term  of  reference. 


Surface  Equivalence  Radius,  or  Mean  Quadratic  Radius 


If  the  total  surface  (section)  of  polydisperse  fog  particles  be 
considered  as  the  determining  property,  then  in  a  similar  way  it  ihay 
be  possible  to  obtain  an  expression  for  the  mean  quadratic  radius: 

Vi  =~-  ^  =  \/‘  J  n  (p)  p  -d?  .  (^7  ) 

o 

Applicable  for  experimental  distributions  is  a  formula  expressed 
in  sums  (in  view  of  their  obviousness  we  leave  out  the  conversion 


operation s ) : 


V 


r  2ntrj 


(48) 


The  mean  size  r2  is  the  radius  of  particles  of  a  monodisperse  fog 
equivalent  to  a  given  poljidi  sparse  fog  in  terms  of  surface  area  (section). 
In  operating  with  the  concept  of  mean  square  radius,  we,  by  virtue  of 
this  very  fact,  are  actually  dealing  with  a  monodisperse  fog  such  as, 
in  case  of  the  invariability  of  the  particle  number, 
possesses  a  total  surface  area  (section)  Identical  to  that  of  a  poly¬ 
disperse  fog.  This  size  was  used  so  far  in  the  analyses  of  optical 
problems  /9/,  /lo/,  but  may,  in  principle,  be  applied  in  future  to 
Investigations  of  any  phenomena  the  mechanism  of  which  depends  on  the 
surface  area  (section)  of  polydisperse  fog  particles. 

Volume  Equivalence  Radius,  or  Mean  Cubic  Radius. 

Corresponding  to  the  total  volume  of  polydi  sper  se<jL  fog  particles, 
whieh  is  considered  as  the  determining  property,  is  the  mean  cubic  size 
numerically  equivalent  to  ,  , _ 

J  /oo 

•=  y  J  »(r ')?zdo ,  (49 ) 


r3 


where  r,  is  the  dimensional  value  of  this  mean  size. 

For  experimental  distributions  use  is  made  of  the  obvious  approximate 
expression  in  sums 

-  v/  W? 

(50) 


i/M 

3~  V  Zn, 


In  operating  with  the  concept  of  mean  cubic  size,  we  thereby  reduce 
the  problem  to  an  analysis  of  a  monodisperse  fog  such  as,  in  case  of  the 
same  number  of  particles/  as  is  present  in  a  polydisperse  fog,  has  a 


(3  U 


volume  identical  to  that  of  the  latter.  This  size-ls  used  for  the 

solution  of  problems  connected  with  the  study  of  fog  water,  since  the 
numerically 

latter  coincides/pretty  accurately  with  the  volume  of  particles.  The 
size  r-j  has  found  the  widest  and  aon  si  stent  application  in  the  calculat¬ 
ions  of  natural  fog  water  carried  put  by  N.V.  Kucherov  /lo/  and  B.V. 
Kiryukhin  /9/» 

The  Specific-Surface  Equivalence  Radius. 

As  for  any  disperse  system,  the  dispersity  of  fog  may  be  evaluated 
through  specific  surface  ^  which  is  usually  assumed  equal  to  the  surface 
of  a  matter  dispersed  in  the  volume  of  one  cubic  centimeter. 

Fbr  monodisperse  fog,  the  particle  radius  of  which  is.  equal  to  ^  t 
the  specific  surface  ©  will  be  equal  to  the  surface  of  the  drops  which 
after  coalescence  occupy  a  volume  of  1  cm^,  and  may  be  numerically  deter¬ 
mined  from  the  following  relationship 


N't*# 


—  5 

P* 


where  N1  is  the  number  of  monodisperse  fog  drops  with  radius 
forma  1-cm^  volume  after  coalescence. 

In  the  case  of  polydisperse  fog, its  specific  surface 
determined  in  accordance  with  the  definition  from  expression 


(51 ) 

^  which 

will  be 


S 

n 


CO 

4kN'  j  n  (p)  p*dp  3  J  /:  (p)  pMp 
0 


4-*ap  j  «(p)p3^p 
3  0 


CO 

j  « (p)  p3^p 


(52) 


She  two  mentioned  fogs  will  be  equivalent  in  terms  of  specific  surface 
if  one  assumes 


(52a) 


So 


-50a- 


CO 

J 

.  0 


j  n  (p)  p3tfp 


oo 

I 

0 


J  n  (?)  pidp 


(53) 


In  equating  (51)  and  (52)  and  solving  the  thus  produced  equivalence 

equation  with  respect  to  f  ,  we  obtain 

In  making  use  of  (48)  and  (49), it  is  easy  to  produce 


or 


(54) 


(55) 


From  expressions  (jl)  and  (52)  it  follows  that  the  magnitude  of  the  specific 
surface  of  monodisperse  fog  is  determined  only  by  the  drop  radius,  while  the  mag¬ 
nitude  of  the  specific  surface  of  polydisperse  fog  depends  both  on  particle  radius 
as  well  as  on  the  form  of  the  distribution  function. 

The  magnitude  of  the  specific -surface  equivalence  radius  for  empirical  dis¬ 
tributions  is  determined  by  the  self-evident  approximate  formula 


r 


.2  * 


(56) 


In  calculating  the  size  or  r^  by  formulae  (53),  (54),  and  (55) 
if  the  distribution  is  given  in  analytical  form,  and  from  (55)  and  (56) 
for  experimental  distributions,  we  thereby  find  the  size  of  particles 

0J 

contained  in  monodisperse  specific  surface 

identical  to  that  of  the  investigated  polydisperse  fog.  In  all  cases 
when  the  mechanism  of  the  investigated  fog  phenomenon  is  determined  by 
the  magnitude  of  specific  surface,  we  can  —  in  applying  the  concept 


S'/ 


-51- 


V-) 


of  the  specific-surface  equivalence  radius  instead  of  a 

polydisperse  fog  utilize  a  model  of  monodisperse  fog  equivalent  only  in 
this  respect.  Instead  of  the  specific  surface  we  may  be  interested  in 
the  sum  of  cross-sections  per  1  cur  of  fog  liquid  water  (specific  sect¬ 
ion).  In  this  case/  too  ^he  equivalent  monodisperse  fog  will  be  that^ 
the  particle  size  of  which  is  determinable  from  formulae  (53)-(56). 

This  fact  (we  shall  omit  here  the  derivation  of  the  pertinent  formula) 
makes  it  possible  to  considerably  widen  the  scope  of  application  of 
the  specif ic- surf ace  equivalence  radius  .  We  see  that  the  mean  size 
r  corresponds  to  specific  surface g  selected  as  the  determining  property. 
As  on  previous  occasions  this  is  the  radius  of  particles  present  in 
a  monodisperse  fog  having  the  same  specific  surface  as  the  investigated 
polydisperse  fog. 

Houghton  /3l/  was  first  to  use  the  numerical  value  of  size  r  .  He 
has  also  discovered  the  experimental  numerical  closeness  between  this 
size  ana  the  volume  distribution  mode.  With  reference  to  the  volume 
distribution  curve  plot  he  writes  /3l/*" ••••  the  average  diameter  is 
approximately  equal  to  the  abscissa  corresponding  to  the  maximum  of  the 
drop-size  distribution  curve"  (page  71). 

We  shall  see  later  in  the  text  that  the  numerical  closeness  between 

I 

size  r^  and  volume  distribution  mode  rv  could  have^  lndeed(  been  observed 
in  certain  particular  cases,  but  It  Is  not  a  general  statistical  rule. 

Houghton  was  unable  at  that  time  to  reveal  the  meaning  of  this 

mean  size  and  to  substantiate  the  legitimacy  of  its  application  to 
volume  distributions.  This  fact,  then,  became  subsequently  the  direct 
cause  of  the  confusion  involved  in  the  problem 'of  mean  sizes. 

As  it  fallows  from  our  analysis,  the  size  r  with  Its  inherent 
simple  and  obvious  physical  meaning  constitutes  the  mean  of  size  dlstrl- 

So z 


52- 


(32) 

butions  and  does  not  in  an  obvious  way  become  a.  part  of  the  system  of 
volume  distribution  means. 

The  designation  "average"  size  used  by  Houghton  (and  subsequently 
by  other  author s ) , inaccurate  as  it  is,  has  not  only  failed  in  helping 
to  clarify  the  meaning  of  this  mean,  but,  on  the  contrary,  has  lead  in 
a  number  of  cases  /27/f  / 2/ t  /25/,  to  the  formation  of  erroneous  ideas 
about  the  mean  particle  sizes  of  natural  fog  and  clouds. 

We  have  examined  in  detail  the  mean  sizes  r_ ,  r_  ,  r „ ,  v~, ,  and  r  , 

°  1  ’  d  ’  3  ’  & 

which  are  Quite  sufficient  for  the  solution  of  basic 

problems  arising  in  the  process  of  statistical  investigation  of  natural 
clouds  and  fog  as  polydisperse  systems  analyzed  in  the  form  of  size  dis¬ 
tributions.  Apart  from  the  considered  mean  sizes,  use  is  made  in  some 
studies  of  the  volume  distribution  mode  (see  Section  1)  expressed  in 
linear  measures.  This  is  the  size  already  familiar  to  us:  rv  -  the 
predominant  radius  in  terms  of  volume.  The  moot  probable  volume  vQ 
of  distribution  (12)  corresponding  to  the  maximum  of  this  distribution 
may  be  found  as  the  only  root  of  equation 

A.  n  ( v )  ==  0.  ( 56  a ) 


Then  with  the  aid  of 
from  the  obvious  formula 


vQ  the  predominant  volume  radius  may  be  found 


b 

(56  ) 


but  it  does  not  coincide  with  any  of  the  mean  sizes  in  a  general  case. 

The  predominant  volume  radius  is  of  no  interest  for  purposes  of  practical 
utilization  since  it  is  not  a  mean  of  size  distribution.  The  predominant 
surface  size  rs  is  not  directly  discussed  in  literature,  and  will,  there¬ 
fore,  not  be  discussed  here. 


In  conclusion  we  shall  pass  to 


question  concerning  the  relationship  of  the 


£3 


.(33) 


distribution  means  of  particle  sizes  in  natural  fog  and  clouds. 


8  .  The  Relationship  of  Mean  Sizes . 


We  have  already  shown  In  Section  6  that  the  distributions  of  parent 
population  particles  in  clouds  and  fogs  must  necessatily  be  considered 
as  single-peak  distributions  with  a  continuous  variate.  The  form  of  the 


law,  however,  remains  unknown 

fc>S<2 


There  arises  a 


-  to  define  a  1  c^w  of  distribution  of  the 


parent  population,  i.e.  to  find  an  analytical  curve  which  could  in  the 

best  possible  way  fit  an  empiric  distribution  or  a  group  of  them  if  all 

of  them  are  obtained  from  clouds  structurally  similar  in  form.  It  is 

smoothing 

known  from  mathematical  statistics  that  a  curve  capable  of 
an  experimental  distribution  in  the  best  possible  way/ 

i  s  who 

theoretical  and  empirical  moments  coincide.  The  theoretic  moments  in 
our  case  are  the  mean  sizes  p p2>  p^>  811(1  *  -^'rie  sarae  sizes  result¬ 

ing  from  approximate  relationships  (46),  (48),  (50),  and  (56)  may  be 
considered  as  empirical  moments.  In  solving  the  problem  of  finding 

such  a  curve  it  is  necessary  to  have  a  sufficiently  general  class  of 

depending  on  the  least 

single-hump  curves  satisfying  conditions  (2)  and 
possible  number  of  parameters.  In  the  capacity  of  such  a  class  use  was 
made  of  class  (19)  previously  utilized  in  Section  2,  and  exhibiting  a 
generality  sufficient  for  our  purposes. 

We  shall  find  the  mean  sizes  for  this  class  as  functions  of  para¬ 
meters  m  and  n. 


It  is  evident  that 


d  V=\  -rMN  =  izCNpm+3e  «  p"  dp, 

where  V  -  volume  of  distribution  particles  numerically  equal  to 


(56c) 


rC/V  f  pm+3e  n  '  dP. 


In  replacing  the  variables  by  formula 


we  obtain  for  (57) 


4 -~CN  f  Pm+3e  n 
o 


■  4 

do  =  -~~CN 


m  +4  OO 
n  \  „  /* 


f  zm+3e~indz-= 


i  -J-4  p  / m  ~f~  ^ 


(57aQ 


(57b) 


In  substituting  th6  expression  C  from  (19)  in  the  above  equation. 


we  produce 


3  r  (m  +  4  \ 

'=  —  *N  "  —  -  -  l 

3  \m)  r  i^n  -f  1  ^ 


(58) 


On  the  other  hand  it  is  evident  that 


^=i-^33- 


(58a) 


In  equating  the  above  equation  and  (58),  we  shall  obtain  the 
equivalence  condition  in  solving  which  with  respect  to  ^  we 
shall  find  the  expression  for  the  mean  cubic  size  of  class  (18)  curves 


in  a  general  form 


3  / 

n  N  «  n  /  V  n  1 

m)  V  r(*±i)' 


Similarly  in  considering  surface  and  size  distributions  it  is 
possible  to  obtain  an  expression  for  the  arithmetic  mean  and  mean 
quadratic  sizes.  On  the  basis  of  (5^),  when  p  ^  and-  are  known, 
it  is  easy  to  determine  the  specific-surface  equivalence  radius  "j? 

In  omitting  the  intermediate  calculations,  vie  shall  write  the 
generalized  expressions: 

for  the  arithmetic  mean  size 

i  r(m  +  2\  (60) 

o  =  ( nY"" 

11  \m  )  (m  -f  1  \  ’ 


for  the  mean  quadratic 


11  \m  )  r  /m  +  l  \  ’ 
ilc  \  n  ) 

JL  /  r  i m  ^  3~T 

p.i==(lL\"  1/  [ 

' OT '  f  r  (’”* +  i 


(61) 


\m  /  y  p  ,  m  -f  1  \ 

and  for  the  specific- surface  equivalence  radius 

i_  r('w  +  4^ 

-  /  n  \„  \  n  J 


r(^) 


(62) 


The  mode  of  volume  distribution  j?y,  expressed  in  fractions  of  the 

most  probable  radius  for  class  (19)  curves  may  be  determined  in  an 

elementary  way  and  turns  out  to  be  equal  to 

o  1"/r^+~3  (6  3) 

*  v  V  m 

A  numerical  coincidence  of  sizes  and  observed  by  Houghton 

for  his  experimental  curves,  which,  possibly,  has  served  as  a  basis 

for  treating  these  sizes  as  identical,  occurs  for  class  (19)  curves  at 
n  *  1.  Houghton's  conclusion  anent  the  numerical  coincidence  of  sizes 
and  is  not  contradicted  by  the  experimental  data  produced  by  N.V. 
Kucherov  /lo/,  I«I*  Ghestnaya  /22/,  and  V.^.  Zaytsev  /l / •  We  shall, 
therefore,  seek  the  best  fitting  curve  among  the  n=l  family  of  class 
(19)*  For  this  purpose  we  shall  identify  certain  properties  of  the 
family  of  curves  we  are  interested  in.  We  shall  find  the  expressions 
for  the  mean^  sizes  of  the  family  of  curves  n  =  1  of  class  (19). 

In  taking  advantage  of  the  well  known  formula  T(  x+l  )  -  xi ,  from 
which  it  follows  that 


1' (*-M) 
r  (a) 

(63a 

we  shall 

obtain  instead 

of  (60),  (59),  (61),  and  (62) 

the  following 

preci  se 

and  simple  expressions 

rn  1 

(6*) 

-  _  /(m  2)  ( m  +  1) 

‘  2  m  ’ 

3 

•(65) 

7  V  (rn  -f  3)  (rn  +  2)  (m  +  1) 

*'8  '  m 

(66) 

m  +  3 

13  m 

(67) 

We 

shall  introduce 

an  auxiliary  parameter  £  m  2m 

and  assume 

(rn  +  2)  (m  +  1 )  ss  (m  -f  ^ )2 

(6$a 

and 

(rn  -j-  3)  (rn  -j-  2)  (rn  -j-  1)  ss  (m  -f-  2)3. 

(67b 

-5~6 


Then  instead  of  (64),  (6 5),  (66),  and  j(67 )  we  shall  obtain  even 
simpler  expressions  for  the  mean  sizes: 


Here  (68)  and  (71)  are 


(68) 

Pi — 1  +  2s,  r,  =  r0  (1  -f-2e), 

p2=I+36.  r2  =  ro  ( *  -f"  3s), 

(69) 

P3=1+4s,  r3  ss  r0(l  +4s), 

(70) 

P,=  l+6e,  r3  =  r0(  1  H-  6s). 

(71) 

exact 

expressions.  Inspite  of  the  fact 


that  (69)  and  (70)  are  approximate,  the  resulting  error  for  m  =  2  does 

2  -2P 

not  exceed  2%  .  For  example,  for  curve  p  e  used  by  A*Kh.  Khrgiyan  and 
I.P.  Mazin  /2l/,  the  latter  formulae  give  ^  =  1*50,  =  1*75,  ?3  '  1'50' 

=  2.50,  while  the  exact  values  are  but  slightly  differing  and  are 

T  <S 

equal  to  p  -  1.50,  p2  =  1«73,  ^  =1.96,  =  2.50. 

At 'higher  m  values  the  error  becomes  negligibly  small. 

An  important  property  of  the  n  =  1  family  of  class  (19)  curves  follows 

from  (68),  (69),  (70),  and  (71)*  In  passing  from  one  mean  size  to  another 

of  ^  1 

they  change  by  a  definite  value  multiple  the  integer  6=  Urn-  Moreover, 
the  arithmetic  mean  size  for  the  n=l  family  is  always  larger  than  the 
most  probable  size  by  26rQ,  and  the  specific-surface  mean  -  by  6&rQ  .  This 
enables  the  investigator,  if  the  mean  sizes  are  know,  to  produce  a 
definite  answer  to  the  question  of  whether  a  given  experimental  distri¬ 
bution  may  be  considered  as  satisfying  the  n=l  family  of  class  (19). 

It  is  well  known  that  the  most  probable  size  for  arbitrary  experiment¬ 
al  distributions  cannot  be  calculated.  We  shall  attempt  to  produce  an 
analytical  expression  of  r0  for  the  n=  1  family  for  the  case  of  any  m's. 

It  follows  from  (680,  (70  ),  and  79)  that: 


p.  -  p,  =  4c, 
P,  “  Ps  ~  2s, 


Whence  , 


(71a) 

(71b) 

(71°) 


=  'J,  —  Pi  +  2  (p.  ~  P'3 )  +  “2  (p3  --  Pi)  s  2  (p,  -  -  Pi)  ,  (T1  ) 


1  - 


and  if  one  considers  (64)  and  (67),  then  instead  of  an  approximate 


expression  it  becomes  possible  to  write  an  exact  expression  in  the 


form  of 


(71e) 


Hence,  for  the  most  probable  unit  size  value  one  may  write 

3  /—  ■  -  ,  3pi  —  p 


6e  =  o 


l  (pj  ~  Pi) 


■=■  1. 


(71f ) 


In  passing  on  to  dimensional  presentation,  we  produce  an  express¬ 
ion  for  the  most  probable  size 


(72) 

Formula  (72)  is  exact  from  the  analytical  point  of  view.  Stochastic¬ 
ally  its  accuracy  is  not  inferior  to  that  of  the  well  known  formulae 
(46),  (48),  (50),  and  (55)- 

In  making  use  of  formula  (72),  and  utilizing  the  known  values  of 
r,  and  r -  ,  the  most  probable  size  may  be  determined  more  precisely 
than  from  a  distribution  curve. 

Incorporated  in  Table  5  are  the  calculated  values  of  the  mean 
sizes  for  different  m's 


Table  5 


p 

m 

i 

2 

3 

( 

4  i 

5 

6 

7  1 

8  ; 

9  1 

10 

p  I 

2,00 

1,50 

1,33 

1 ,25 

1  ,20 

1,17 

1,14 

1.12 

1,11 

1,10 

p2 

2,22 

.  1,73 

1,49 

1 ,37 

1,30 

1 ,25 

1,21 

1,18 

1,17 

1,15 

?3 

1  2,88 

1 ,96 

1,64 

1 ,49 

1,39 

1,33 

1 ,28 

1 ,24 

1  ,22 

1,20 

P. 

' 

4,00 

i 

2,50 

2,00 

1 ,75 

1 ,60 

1 ,50 

1,43 

1,38 

1,33 

1,30 

In  order  to  find  the  equation  of  a  curve  satisfying  the  experimental 

distribution,  it  is  necessary^  in  using  formula  (72  )^  to  determine  the 

empirical  moments  £*■  .  and  ,  and  to  make  sure  that  they  are 

Td  >  r0  )  To 

compatible  with  (64)  -  (67)  with  respect  to  their  differences.  If  the 
increments  of  the  moments  satisfy  (64)  -  (67)  then  it  can  be  said  that 


the  experimental  distribution  Is  described  by  the  family  of  curves 

n  =  1  * 


S? 


Thereupon  one  looks  for  such  an  Integer  m  in  Table  5  for  which 
the  mean  size  ratio^  complies  best  of  all  with  the  experiment. 

If  there  are  no  doubts  that  the  experimental  distribution  belongs  to 
the  n  s  1  family,  a  simpler  method  of  parameter  m  determination  may  be 
used  and,  consequently,  also  of  determining  the  analytical  expression 
of  the  curve  for  the  experimental  distribution.  The  essence  of  this 
method  consists  in  the  following.  In  utilizing  (64)  and  (67)  a  ratio 


of  moments  or  mean  sizes  is  composed 

9,  h 


!n  =  ^i-=_  .!±±A  (72&) 

7,  r,  «  +  172  ' 

of  which  for  different  m's  are  given 


the  numerical  values 
in  Table  6 . 


Table  6 . 


III 

1 

2 

! 

3  : 

4 

I 

5 

6 

7 

s 

9  1 

10 

m  — |—  3  j 
rn  +  1  ! 

2,00 

1,67 

1 

1  .50 

1 

1  ,40 

1,33  ! 

1  ,29 

1  ,25 

1  ,22 

1  ,20 

1  .  IS 

In  composing  the  ratio  of  sizes 


for  the  analyzed  experimental 


distribution,  which  is  numerically  equal  to  jk,  and  making  use  of  Table 
6,  It  is  easy  to  find  an  integer  m,  and,  consequently,  the  corresponding 
equation . 

Both  above  methods  are  variants  of  the  method  of  moments.  The 
first  method^  is  described  by  us  for  the  sake  of  consistency  and 
ness  of  the  discussion,  the  second  for  purposes  of  practical  application. 
Both  method  are  practically  equivalent  from  the  point  of  view  of  accuracy 
and  yield  identical  results.  This  is  attributable  to  the  fact  that  p^. 
is  a  moment  of  a  higher  order  thafi  and  ^  and, in  addition  to  that  , 
in  accordance  with  (54)  it  constitutes  an  explicit  function  of  and 
■^2  and,  therefore,  also  of  both  r^  and  r g  which  are  not  explicitly 
utilized  in  the  second  method. 

After  the  selection  of  an  integer  m  from  Table  6  the  particle- size 
distribution  function  for  family  n=l  of  class  (19)  may  be  written  in  the 
form  „.(p) 

(73) 


or  in  the  corresponding  dimensional  form,  a  transition  to  which  does  not 


present  any  particular  difficulty. 

_  _  !«£■ 

The  questioniof  the  accuracy  in  determining  ^  ^ ,  p  ,  and  n(|j>) 

from  the  stochastic  standpoint  are  not  considered  in  this  paper  since 
thfj  tran scendj  the  scope  of  the  investigation.  The  more  so,  because  in 
the  form  of  the  goodne  ss-of—  f  it  test  of  the  empiric  and  theoretical 
distributions  they  are  comprehensively  discussed  in  textbooks  on 
mathematical  statistics. 


All  the  published  distribution  functions  for  fog  and  cloud  particles 
/23,  21,  24,  11,  etc./  are  particular  cases  of  the  class  of  curves  (19) 
and  its  family  (73)  proposed  by  us  /5/« 

In  K.S.  ShifrinSand  N.P.  Bogdanovskay  a 1  s  paper  /24/  a  four-parameter 

is 

class  of  curves  fjp  utilized  ,  from  which  a.  curve  of  best  possible  fit 

may  be  selected.  An  increase  in  the  number  of  parameters  in  this  case 

is  related  to  the  fact  that  the  coefficient  at  the  particle  size  in  the 

constant  of 

exponent  (in  our  designations  -  2)  and  the  normalizing  the 

distribution  are  assumed  to  be  _nits  independent  parameters. 


Actually,  as  we  were  able  to  show, both  these  quantities  are  functions 

of  two  independent  distribution  parameters  and  do  not  possess  the 

property  of  Independence.  Hence,  an  increase  of  the  number  of  parameters 

cannot  be  considered  substantiated  in  this  case. 

* 

In  the  method  proposed  by  L.M.  ^evin  /ll/  the  mean  sizes  or  the 
corresponding  moments  of  the  experimental  distribution  are  not  calculat¬ 
ed.  This  simplifies  the  selection  of  a  curve  from  a  one-parameter 
family.  However,  Levin  does  not  offer  a  method  permitting  to  determine 
whether  the  investigated  experimental  distribution  actually  belongs  to 
a  single -parameter  family  of  curves.  In  view  of  this  fact,  and  also 
because  the  geometric  constructions  introduce  an  indefinite  arbitrari¬ 
ness  in  the  matter  of  curve  parameter  selection,  the  potentialities 

a 

o^  Levin  s  method  are  restricted  In  terms  of  accuracy. 


In  conclusion  we  shall  find  the  relationship  which  makes  it  possible 
on  the  basis  of  mean  sizes  to  determine  by  analytical  method  the  value 
of  parameter  m  and,  consequently,  also  to  find  the  equation  of  the  curve 
smoothing  the  experimental  distribution. 

In  utilizing  (68)  and  (71),  we  note  that 

|T  —  1  -f-7  —  p7=  10e.  (73a) 

In  order  to  produce  the  final  expression  this  equation  must  be 

solved  with  respect  to  ^  which  should  be  replaced  by  m  in  accordance  with 

formal  a  <b  =  — —  ,  whereupon  in  dimensional  form  we  obtain 
2  m 

=  - 1.  (74) 

h  -  n 

Formula  (74),  v/hen  the  mean  sizes  of  the  experimental  distribution 
and  r^  are  known,  permits  the  curve  parameter  m  to  be  calculated  analyt¬ 
ically  up  to  n's  y AL.UG&  tN  A  S/*lp Artb  Exact  /,l /) a/ /V £ rZ. ,  In  this 

case  there  is  no  need  to  resort  to  graphic  constructions,  as  this  is 
required  in  Levin's  method  /ll/. 

9.  A  Few  Bxperlmental  Data. 

All  the  deductions  and  corollaries  ensuing  from  the  consideration 
of  the  problem  of  natural  fog  and  cloud  particle  size  determination  by 
sampling  method  which  were  recorded  so  far  need,  of  course^  to  be  compared 
with  experimental  data.,  just  as  the  results  of  any  theoretical  research. 

In  part  this  was  already  done  in  the  process  of  discussion.  For  example, 
in  Section  2  with  regard  the  height-of-curve  evaluation,  in  Section  3 
in  connection  with  the  problem  of  the  finiteness  of  the  interval  width, 
and  in  Section  6  in  analyzing  the  "law  of  multiple  masses".  A  comparison 
of  the  contents  of  other  sections  with  experimental  data  does  not  appear 
possible  owing  to  the  unavailability  of  published  experimental  materials^ 
suitable  for  this  purpose.  The  experimental  data  discussed  in  this 
section  are,  therefore,  designed  to  fill  in  this  gap. 

The  experimental  part  of  the  study  includes: 
a)  verification  of  the  applicability  of  formula  (34)  to  experimental 


distributions  of  fog  and  cloud  particles; 

b)  investigation  of  the  statistical  filling  of  experimental  curve; 

c )  practical  application  of  the  method  of  moments  in  the  analysis 
of  experimental  distributions. 

serial  sampling 

The  first  two  promblems  were  solved  by  means  of  the 
tx  method  well  known  in  statistics.  This  method,  which  so  far  was  not 
applied  in  microphysical  investigation  of  fog  and  clouds,  essentially 
consists  in  the  following.  A  random  sample  of  large  volume  consisting 
of  a  group  or  series  of  small  samples,  also  randomly  taken,  is  drawn  from 
a  homogeneous  set  of  particles.  To  elucidate  the  question  of  the  curve 
of  best  possible  fit,  several  experimental  distributions  were  produced. 

The  mean  sizes  of  these  distributions  made  it  possible  for  the  first  time 
to  utilize  the  method  of  moments. 

All  particle  size  measurements  were  carried  out  by  the  commonly  known 
/l/  method  of  photomicrography  of  particles  suspended  in  a  thin  layer  of 
oil.  The  filming  was  made  through  GOMAZ  microscope  on  a  standard  positive 
cinematographic  film  (80°  KhD)  with  the  aid  of  a  "Sport"  mirror  camera. 
Total  magnification  (magnification  of  filming  and  projection)  was  equal 
to  1000.  Peg  moisture  was  measured  by  means  of  an  aspiration  psychrometer 
by  the  method  proposed  by  O.V.  Levitskiy  /l7 / ,  and  was  found  to  be  equal 
to  100%  in  all  cases. 

The  experimental  data  are  recorded  in  five  appendices  in  the  form  of 
7  distributions.  They  were  produced  on  the  basis  of  63,547  particle 
measurements.  Incorporated  in  appendices  1,  2,  and  3  are  the  data  on 
size  distributions  of  particles  of  advection  fogs  observed  4n  late  fall 
of  19^6  over  l  4l<f  LatOfefl  in  the  region  of  Priozersk.  The  observat- 

fr  O  ny 

ions  were  carried  out  a  sand  bar  protruding  far  into  the  open  lake. 

The  data  illustrating  the  conditions  under  which  these  fog  obser¬ 
vations  were  carried  out  are  furnished  in  Table  7» 


Table  7. 


Append- 
ice  s . 

Air  temperature 
in  degrees 

Surface 

water 

temperature 
in  degrees 

Wind 

Visibi¬ 

lity 

range, 

m . 

at  2-m  al¬ 
titude  . 

at  water 
level 

velocity 

m/sec. 

2 

3,6 

3,9 

5,4 

8-10 

w 

500 

3 

1  ,0 

3,0 

5,0 

2—3 

sww 

500 

1 

—  1,2 

-0,2  i 

4,2 

3-4 

sw 

200 

In  other  appendices  (4  and  5)  are  represented  four  distributions 
obtained  for  radiation  fog  particles  observed  at  the  outskirts  of  Lenin¬ 
grad  (Novaya  Derevnya)  in  summer  and  fall,  1948.  The  radiation  fogs  were 
observed  in  night-time  at  the  fringe  of  a  large  field  without  any  trees 
or  buildings.  As  a  rule,  there  was  no  wind. 

The  mean  sizes  r-^ ,  ?2,  r^  and  r6  were  calculated  for  all  the  plotted 
experimental  distributions. 

Filming  of  30,000  particles  simultaneously 
drawn  from  natural  advection-origln  fog  was  made  experimental  clarifi¬ 
cation  of  the  questions  connected  with  the  application  of  formula  (34). 
This  fog  was  observed  for  several  hours  without  noticeable  variation  of 
visibility  (300  -  500  m)  and  size  of  particle s;  on  November  3,  1946. 
Sampling  for  oil  tests  was  made  <ffn  a  sand  bar  extending  far  into  the  open 
lake, and  the  pictures  were  taken  right  away  on  the  spot. 

The  particle  size  data  were  represented  in  the  form  of  100  separate 
experimental  curves,  each  curve  covering  a  sample  volume  of  300  particles. 
Complete  information  relative  to  these  curves,  as  well  as  the  general 
distribution,  will  be  found  in  Appendix  1.  On  the  basis  of  these  data 
deviation  distributions  were  obtained  for  each  i-th  interval.  For  the 
sake  of  brevity  deviation  will  be  referred  to  as  the  difference 

/nij  - 

wherein  n^  is  the  mean  or  most  probable  frequency  value  of  the 

i-th  interval,  and  j  -  the  ordinal  number  of  the  curve  in  the  series. 

The  experimental  values  of  relative  error  Dj_  with  reliability  ei»  were 


determined  from  the  deviation  distributions.  An  error  with  reliability  << 

43 


(determined  experimentally)  was  eoiisidered  to  be  that  deviation  which 
was  exceeded  for  the  investigated  interval  in  100c<-cases.  Relative  errors 
Dj_  with  two  reliability  values  :  «< ^  =  0*68,  and  °<2  “  0  .96  were  experiment¬ 
ally  obtained  in  this  manner.  To  avoid  the  effect  of  supplementary 
errors  arising  from  measurement  of  excessively  fine  drops,  only  the 
intervals  situated  to  the  right  of  the  experimental  curve  maximum  were 
subjected  to  investigation.  The  graph  duff  Fig. 8  shows  the  experiment  .ally 
evaluated  errors  Di  accompanied  by  the  corresponding  curves  of 

relative  errors  calculated  from  formula  (34).  A  comparison  of  the  ex¬ 
perimental  data  with  the  theoretical  curve  shows  that  formula  (34) 
satisfies  the  conditions  of  the  experiment  over  a  wide  range  of  frequencies 

For  the  small  frequencies  (fractions  of  ten)  the  real  errors  were  found 
notice  ably 

to  be  exceeding  the  calculated  values.  ,An  explanation  for  this 

experimentally  observed  phenomenon  is  provided  by  V.I.  Romanovskiy  /l4/. 

He  has  shown  theoretically  that  in  case  of  a  small  number  of  measurements 

the  estimates  of  effors,  in  accordance  with  the  classical  statistical 

higher 

theories,  lead  to  results  characterized  by  ^001  precision.  In 
refering  the  reader  interested  in  greater  details  to  the  above  source, 
we  shall  for  purposes  of  comparison,  reproduce  here  a 

theoretical  curve  plotted  in  accordance  with  Romanovskiy 1  s 

OWM 

and  our^exper iment al  data  which  confirm  these  views  (Fig.9). 


Fig .  8 .  Statistical  errors  of  experi¬ 
mental  distributions  in  relat¬ 
ion  to  the  concentration  of  particles 
in  the  interval  and  the  reliability  pa¬ 
rameter. 


Fig  .9  Statistical  errors 

for  the  case  of  small 
frequency  ratios. 


We  can  thus  see  that  formula  (34)  is  entirely  suitable  for  a  practic¬ 


al  estimate  of  errors  In  those  intervals  of  ex¬ 
perimental  distributions  whose  frequencies  are 
not  too  small*  For  frequencies  of  the  order  of 
fractions  of  ten  formula  (34)  yields  appreciably 
underrated  results* 

Let  us  turn  now  to  experimental  distribut¬ 
ions  in  which  it  is  possible  to  trace  the  process 
of  statistical  filling.  As  a  typical  example  we 
shall  examine  the  distribution  obtained  for  the 
particles  of  a  radiation  fog  observed  in  the 
night  of  Aujaust  12,  1948  at  the  outskirts  of 
Leningrad  in  the  region  of  Novaya  Derevnya  (Ap- 

Flg  *10 »  Experimental  pendix  4).  This  distribution  with  the  sample 
distribution  (1000 

particles)and  its  16  volume  involving  1000  particles,  and  a  series  of 
components. 

component  curves,  each  with  100  p article  s,  are  shown  in  Fig. 10* 

The  general  distribution,  as  it  may  be  seen  from  Fig  .10 ,  bas  a  single 
peak,  whereas  the  individual  component  distributions,  almost  all }  have 
several  peaks.  Essential ,  however ,  is  the  fact  that  each  of  these  10 
curves  does  not  resemble  the  initial  distribution  of  which  each  of  them 
is  a  random  sample.  Concentrated  about  the  maximum  of  the  full  curve  are 
about  100  particles.  The  maxima  of  the  individual  curves  of  the  series, 
each,  on  the  average ,  contain*  about  10  particles.  The  calculated  relative 
error  in  the  maximum  determination  in  the  case  of  statistical  filling 
diminishes  from  64%  for  10  particles  to  20%  for  100  particles  per  interval. 
The  degree  of  statistical  concentration  in  each  of  the  10  curves  may  be 
considered  as  the  beginning  of  the  main  shaping  stage  (II  stage). 

Experimental  curves  in  the  initial  shaping  stage  (I  stage)  may  be 
obtained  and  similarly  analyzed,  if  one  of  the  10  curves  with  a  sample 
volume  comprizing  100  particles  be  visualized  as  consisting,  in  turn,  of 


a  series  of  experiments!  curves  with  a  still  smaller  sample  volume. 

We  shall  refrain  from  doing  so  because  there  is  no  need  for  this.  Among 
the  published  experimental  distributions  there  is  a  considerable  number 
of  curves  the  formation  of  which  was  discontinued  at  the  initial  shaping 
form.  Out  of  a  total  of  134  experimental  curves  produced  by  M.  Diem  /27/ 
/28/,  40  are  plotted  on  the  basis  of  measurement  of  50  particles  and 
less,  as  this  may  be  seen  from  Table  3» 

The  position  of  the  interval  with  the  greatest  concentration  of 
particles  varies  in  the  process  of  statistical  filling. 

The  interval  containing  the  mode  of  the  distribution  becomes  finally 
evident  by  the  end  of  the  main  shaping  stage.  The  maxima  of  accidental 
multimodality  in  the  process  of  filling  of  the  experimental  cum  (after 
the  basic  maximum  has  been  finally  formed)  shift  into  the  region  of 
intervals  with  lesser  statistical  concentration.  This  may  be  observed 
on  all  experimental  curves  with  a  sufficient  sample  volume.  The  sub¬ 
sequent  displacement  of  accidental  maxima  into  the  region  of  the  curve 

in 

tails  may  be  observed  the  distribution  of  30,000  particles,  which 

was  already  discussed,  and  also  in  the  cases  of  two  curves  with  a  larger 
sample  volume  (Appendices  2  and  3)* 

Thus,  the  results  obtained  in  the  study  of  the  process  of  experiment¬ 
al  curve  formation  by  method  of  serial  samples  (the  results  v/ere  partly 
represented  in  Fig. 10)  fully  confirm  the  character  of  statistical  fill¬ 
ing  discussed  in  Section  6. 

In  conclusion  we  shall  attempt  to  apply  the  method  of  moments  to 
clarify  the  question  concerning  the  shape  of  the  curve  of  the  best  fit 
for  distributions  in  fog. 

Incorporated  in  Table  8  are  the  ratios  of  the  mean  sizes  obtained 
experimentally _ 

by  us  for  advsction  and  radiation  fogs 

(Appendices  1-5),  and  also  the  values  of  parameter  m,  calculated  for 
our  distributions  from  formula  (^4). 


Table  8. 


Moment  or 
mean  size. 


Pi 

P_2 

P3 

P. 


S 


advection  (1946  ) 


m 


calc  . 


(74) 


radiation  (1948) 


1,35 

1 ,35 

1,24 

[ 

1 ,23 

1,17 

I 

1,12 

1,10 

1,55 

1,51 

1,35 

1 ,34 

1,27 

1,19 

1,15 

1,76 

1,68 

1,40 

1,45 

1,34 

1  ,25 

1,17 

2,23 

2,06 

1,73 

1 ,69 

1,50 

1,37 

1 ,30 

2,1 

2,8 

4,1 

4,9 

6,0 

8,0 

9,9 

3 

2 

1 

4 

1 

1 

5  a 

56 

5b 

Appendix  No  * 

Comparison  of  Table  8  experimental  data  with  the  calculated  data 
of  Table  5  shows  a  good  agreement  of  the  mean  size  ratios. 

The  data  of  Table  8  reveal  a  substantial  difference  in  the  relation¬ 
ship  between  the  mean  sizes  for  radi action  and  advection  fogs.  Two  curves 
are  plotted  in  Pig.  11,  one  of  which  £  m  =  6  )  corresponds  to  the  moderate¬ 
ly  skewed  distributions  (radiation  fog)}  the  other  (  m  =  3  )  complies 

manifesting  a 

with  the  experimental  di stributions  ^(((fconsiderable  skew^ness  (advect¬ 
ion  fog ) . 

The  experimental  curves  for  fogs  of  advection  nature  are  constructed 

in  the  form  of  volume  distributions  and  have  a  mode  displaced  to  the 

right  with  respect  to  ”=  2>  while  for  radiation  fogs  this  ratio  is 

ro 

smaller  averaging  about  1.5»  In  utilizing  these  relationships,  and 
bearing  in  mind  the  fact  that  Houghton  /3l/  plotted  his  volume  distributions 
precisely  for  fogs  of  advection  origin,  we  may  estimate  the  most  probable 
particle  size  of  American  fogs  at  10  -  12*5  microns  (in  radius)  instead 
of  4o  -  50  microns  as  indicated  by  Houghton  /3l/  for  the  predominant 
diameter  in  terms  of  volume.  Houghton's  data,  therefore,  indirectly  con¬ 
firm  the  presence  of  a  considerable  skew  in  the  size  distributions  obser¬ 
ved  by  him  in  advection  fogs. 

The  results  of  our  evaluation  agree  with  the  results  of  a  recalculat¬ 
ion  of  the  American  curves  undertaken  by  EUV.  Kiryukhin  /9/. 

As  to  the  structure  of  clouds,  it  is  known  that  the  most  symmetrical 
distributions  are  observed  for  St,  and  the  most  skewed  ones  -  for  Ns. 
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The  other  cloud  forms  occupy  an  intermediate  position. 

V.&.  Zaytsev  /l / ,  for  cauliflower  clouds,  and  A«M.  Borovikov  /l/, 
for  the  remaining  forms,  have  established  the  fact  that  the  bases  of 
clouds  with  noticeable  vertical  development  are  characterized  by  an 

almost  symmetrical  distribution  curve.  For 


the  central  and  top  parts  of  clouds,  however 

,nl  &tL*Cpmemf  a  distribution  manifesting  a  moderate  or 

noticeable  skew  in  the  direction  of  large 

drops  is  to  be  registered.  This  phenomenon 

p  i  s  so  clearly  displayed  that  it  is  possible 

Fig .11 .  Particle  distiri-  to  observe  more  symmetrical  distributions 
bution  functions: 

I  -  radiation  fog;  II  -  (1,  1*08,  1.12,  1.15,  1*23)  at  the  lower 

advection  fog. 

boundary  of  Ns  than  at  the  upper  boundary  of  St  (1,  1.16,  1.24,  1-32, 
1*57)*  The  parenthesized  me  an  sizes  result  from  the  observations 
carried  out  by  A*M.  Borovikov.  The  sample  volumes  are:  1169  and  1029 
drops  for  Ns  and  St,  respectively. 

With  regard  to  the  shape  of  the  experimental  distribution  curves, 
we  were  unable  to  establish  a  difference  between  the  clouds  and  fogs 
on  the  whole.  Nevertheless,  a  conclusion  suggests  itself  about  the 
structural  similarity  between  the  lower  boundary  zone  of  clouds  with 
a  considerable  vertical  development  and  radiation  fogs,  on  the  one  hand 
and  between  the  middle  and  upper  sections  of  the  same  type  of  clouds 


and  advection  fogs,  on  the  other  hand. 

The  number  of  observed  cases  is,  however,  insufficient  to  consider 
the  question  of  the  mentioned  structural  differences  as  satisfactorily 
answered.  It  is  hoped  that  further  accumulation  of  empiric  data  and 
the  application  of  the  method  of  moments  will  enable  future  investigat¬ 
ors  to  estimate  the  structural  characteristics  of  fogs  and  clouds,  not 
only  qualitatively,  but  quantitatively  as  well,  depending  on  the  condit 
ions  of  their  formation. 


Conclusion  s« 

The  experimental  problem  of  size  determination  of  patural  fog  and 
cloud  particles  by  sampling  method  was  dkscussed  in  this  paper  from  the 
positions  of  statistical  probability.  The  principal  results  and  conclus¬ 
ions  of  this  research  consist  in  the  following: 

1*  It  has  been  established  that  the  American  method  of  experiment¬ 
al  curve  construction  developed  in  Houghton’s  works,  in  contrast  with 
the  method  adopted  in  the  USSR  and  in  Europe,  is  statistically  inconsist¬ 
ent. 

2.  It  has  been  shown  that  the  evaluation  of  the  experimental 

distribution  peakedness  as  carried  out  by  Houghton  is  ’wrong. 

experimental 

3*  A  new  characteristic  —  the  height  of/curve  --  which  was  utilized 
in  this  paper  in  the  solution  of  a  number  of  concrete  problems,  has  been 
introduced.  As  a  result  of  a  comparison  with  the  experimental  data  it 
has  been  shown  that  a  value  approaching  unity  may  be  used  as  the  most 
probable  value  of  the  altitude  of  an  experimental  curve. 

4.  a  new  class  of  single-peak  curves,  whose  shape  is  a  function 
of  two  parameters,  was  proposed.  For  this  class  of  curves  manifesting 
a  peakedness  close  to  unity  a  relationship  Was  demonstrated  between  the 
height  of  the  curve  and  its  width. 

5*  The  influence  of  the  finitne  ss  of  the  interval  width  was  investi¬ 
gated,  as  a  result  of  which  an  effective  selection  both  of  the  interval 
widths  and  the  number  of  intervals  was  found  to  be  possible  for  natural 
fog  and  cloud  experimental  distributions. 

6.  The  effect  of  the  finitne ss  of  the  sample  volume  was  investigated. 
It  was  shown  thereby  that  the  frequencies  nj_.  If  the  intervals  are  small 
enough,  may  be  appoximately  considered  as  independent,  and  &  polydisperse 
fog  or  cloud  may  be  assumed  as  consisting  of  a  finite  number  of  monodis- 
perse  components  'jr 

C The  possibility  of  estimating  the  accuracy  at  a  given  point  of 

£>? 


of  an  experimental  curve  with  due  consideration  of  measurement  reliabili¬ 


ty  was  substantiated  for  the  case  when  the  width  of  the  intervals  is 
smaller  than  the  most  probable  size. 

7 •  An  analysis  of  experimental  distributions  has  shown  that  an 
incorrect  statistical  approach  to  the  problem  of  particle  size  determinat¬ 
ion  for  natural  fogs  and  clouds  often  leads  to  the  publication  of  data 
which  turn  out  to  be  inadequate  in  view  of  the  extremely  small  number 
of  measured  drops. 

8.  Xa  investigating  the  amplitudes  of  multimodal 

)t  N  THAT 

experimental  di stri button}  /  K&'bler  has  wrongly  untilized  the  errors 

with  low  reliability  and  has  infringed  upon  the  principle  of  the  random¬ 
ness  of  samples  which  caused  him  to  formulate  an  erroneous  hypothesis. 

9*  A  sufficiency  test  w^s  established  permitting  rational  select¬ 
ions  to  be  made  of  the  sample  volume  for  a  newly  plotted  empirical  curvey 
depending  on  the  statistical  error  in  the  interval  of  the  greatest 
particle  concentration,  the  ’width  of  the  interval,  reliability  of 
measurement,  and  the  most  probable  size  and  height  of  curve. 

10.  Proposed  was  a  method  of  experimental  curve  multimodality  analys¬ 
is  with  the  view  to  establishing  the  character  of  the  curve  (method  of 
the  error  envelopes),  which  is  a  further  development  of  Prof.  A»K. 

Mltropol ' skiy ' s  ideas. 

11.  The  problem  concerning  the  multimodality  of  experimental  dis¬ 
tributions  in  connection  with  the  character  of  their  statistical  filling 
was  investigated.  The  process  of  complete  shaping  of  the  experimental 
curve  was  examined  and  found  to  be  consisting  of  three  stages  substant¬ 
ially  differing  by  character  of  the  sample's  statistical  repre  sentative- 
ne  ss. 

12.  It  was  demonstrated  that  the  multimodal  character  of  the  over¬ 
whelming  majority  of  experimental  distributions  involving  the  particles 
of  natural  fogs  and  clouds  is  not  authentic.  The  reason  for  the  experi- 
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mentally  observed  accidental  mul timodality  is  to  be  sought  in  the  inade¬ 
quacy  of  the  sample  volume^ 

13*  In  accordance  with  the  basic  principles  of  mathematical  statist¬ 
ics  a  new  physically  substantiated  system  of  size  distribution  means  has 
been  developed. 

14.  The  mean  sizes  were  for  the  first  time  used  in  finding  a  curve 

an 

,.rhich  in  the  best  possible  manner  smooths  inve  stigated  experimental 

distribution  for  the  case  of  natural  fogs  (method  of  mean  size  ratios, 
or  method  of  moments). 

15*  The  method  of  19  e  an  size  ratios  was  applied  to  establish  the 

radiation 

structural  distinctions  between  distributions  for  fogs  of 
and  advectiofl  origin.  A  comparison  with  vertically  developed  clouds 
has  shown  that  radiation  fogs  are  structurally  similar  to  the 

lower  boundary  zone  of  clouds,  whereas  advection  fogs  are  analogous 

to  the  central  and  upper  portions  of  clouds. 

16.  Some  results  described  in  this  paper  were  compared  with  the 
experimental  data  published  in  the  studies  of  other  authors.  A  compre¬ 
hensive  statistical  survey,  which  comprizes  the  experimental  part 

of  this  study,  was  carried  out  in  order  to  check  the  principal  results. 
The  method  of  serial  samples  was  selected  as  the  working  method  and  was 
found  to  have  entirely  proved  its  value.  Comparison  with  experimental 
data  has  shown  satisfactory  agreement. 
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RESULTS  OP  AN  ANALYSIS  OP  A  WIND  TUNNEL  USED  IN 
A  MAIN  GEOPHYSICAL  OBSERVATORY 
By 

P.  V.  D'yachenko  and  A.  I.  Kameneva 

The  data  concerning:  the  GGO  wind  tunnel,  which  is  used 
as  a  calibrating:  device  in  the  hydrometeorological  service 
of  the  USSR,  are  explained  in  detail  in  this  article. 

The  GGO  wind  tunnel  occupies  a  special  position  in  the 
hydrometeorological  service.  It  has  been  used,  and  is  being 
used,  not  only  for  the  ordinary  calibration  of  instruments, 
measurements  of  wind  velocity,  but  also  for  various  types  of 
experimental  investigations  and  scientific-method  develop¬ 
ments. 

The  tunnel  was  constructed  in  IC32  under  the  direction 
1 

of  M.  A.  Deraent'ev  .  It  was  first  analyzed  by  A .  I.  Chuvashev 
and  Z.  F.  Val'ta.  Subsequent  analyses  in  1935>  1937.  and  19^1 
were  performed  by  K.  A.  Auskulat,  M.  A.  Agapitov,  A.  I.  Kameneva, 
and  M.  M.  Klestova.  These  periodic  analyses  were  done  in  order 
to  verify  the  invarialibility  of  its  characteristics.  In  this 
instance,  low  velocities  were  not  studied.  Low  velocities  were 
included  in  the  research  program  for  1952. 

The  results  of  the  latter  analysis  of  the  GGO  wind  tunnel, 
which  was  conducted  by  the  authors  of  this  article,  are  briefly 
described  below. 

The  results  of  the  above-mentioned  previous  investigations 
are  also  cited  below  in  order  to  provide  a  comparison  of  the 
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parameters  of  the  tunnel.  These  results  were  taken  from 
manuscripts  and  archive  materials. 

A  Brief  Description  of  the  Wind  Tunnel 

The  GGO  wind  tunnef,  a  schematic  drawing;  of  which  is  shown 
in  Figure  1,  is  of  the  open  tunnel  type.  It  has  an  open 
working  part,  which  is  located  in  a  closed  chamber  1,  the 
dimensions  of  which  are  4.5  x  4.0  x  2 . o  m.  The  tunnel  is 
mounted  in  a  larger  casing,  whose  dimensions  are  17  x  7  x  6  m. 
Its  basic  component  parts  are  the  tunnel  collector  or  nozzle  2, 
which  is  the  suction  part,  located  before  the  working  part;  the 
working  part  3,  where  the  verifying  and  calibrating  instruments 
are  mounted;  the  duffusor  cone  4,  which  is  the  intermediate 
part  between  the  working  part  and  the  air  screw,  and  with  the 
wide  part  facing  the  air  screw;  a  six-bladed  air  screw  5, 
producing  a  determined  variation  in  the  pressure  necessary 
for  the  maintenance  of  the  velocity  in  the  tunnel;  a  return 
channel  or  tube  6  through  the  casing  in  which  the  tunnel'  is 
mounted;  and  an  electric  drive  (converter,  motor,  and  switch 
panel)  safeguarding  the  spin  of  the  air  screw. 

The  collector,  determining  the  form  of  the  working  cross 
section  of  the  tunnel,  has  an  octahedral  cross  section.  The 
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diameter  of  the  working  cross  section  of  the  tunnel,  measured 
along  the  circumference  inscribed  in  the  octahedron,  is  equal 
to  1  m. 

The  operation  of  the  wind  tunnel  consists  of  the  following. 

The  air  screw,  when  set  in  motion,  sucks  the  air  from  the 
casing  in  which  the  wind  tunnel  is  mounted.  The  air  stream 
passes  through  the  collector,  containing  the  rectifying  grid;  enter 
the  working  part  at  first,  then  it  goes  into  the  diffusor  cone, 
from  where  it  is  deflected  by  the  air  screw  back  into  the  same 
casing  which  serves,  in  this  fashion,  as  a  channel  for  the 
return  of  the  air  back  into  the  collector. 

The  six-bladed  air  screw,  which  sets  the  air  in  motion, 
is  powered  by  a  direct-current  6.5  kilowatt  (220  volt)  electric 
motor.  This  electric  motor  is  fed  through  a  converter  from  a 
220  volt,  alternating  current  supply  system.  The  electric 
supply  is  regulated  through  a  special  switch  panel.  The  complete 
assemby  makes  it  possible  for  the  electric  motor  to  operate  in 
ia*  very  low  revolutions,  being  decreased  2.5  times  at  the  rotary 
drive  of  the  screw.  Owing  to  this,  it  is  possible  to  obtain 
very  low  velocities  in  the  tunnel. 

The  ability  to  obtain  small  air  stream  velocities,  as  well 
as  even  changes  of  velocity  in  the  entire  range,  is  also  a 
basic  feature  of  the  GGO  wind  tunnel. 

The  ease  of  mounting  the  instruments  in  the  air  stream,  the 


ease  of  looking  out  after  them,  and  the  remote  control  through 
the  air  screw  motor,  which  makes  it  possible  to  truthfully 
determine  any  velocity  found  in  the  closed  chamber  of  the 
tunnel,  should  be  added  to  the  positive  operation  qualities 
of  the  wind  tunnel. 

The  maximum  velocity  of  the  air  stream  produced  in  the 
tunnel  consists  of  20  m/sec. 

The  basic  requirements  for  an  artificial  air  stream  in 
the  tunnel,  as  concerns  accuracy  and  reliability  of  experimental 
results,  are  a  uniform  distribution  of  the  velocities  in  the 
working  field  of  the  tunnel  with  respect  to  different  cross 
sections,  and  a  preservation  of  the  uniformity  of  the  air  stream, 
both  in  quantity  as  well  as  directional  velocity. 

The  pressure  and  velocity  are  distributed  uniformly  with 
respect  to  magnitude  in  first-rate  wind  tunnels  with  any  cross 
section,  and  should  not  deviate  from  the  mean  magnitude  by  more 
than  +  1  With  respect  to  direction,  the  air  stream  should  be 

parallel  to  the  axis  of  the  tunnel.  The  maximum  deviation  of 
the  angles  of  taper  in  the  vertical  plane  should  not  be  more  than 
0.25°^2^.  .A  change  of  these  parameters  should  bring  about  a 
zero  value. 

It  is  necessary  to  note,  however,  that  an  evaluation  of  the 
nature  of  the  stream  produced  in  the  wind  tunnel  will  be  dependent 


upon  the  accuracy  of  the  instruments?  used  for  measuring  . the 
velocity.  In  our  analysis,  which  was  devoted  to  the  measurement 
of  velocity,  we  used  the  TsAGI  micromanometer  in  pair,  coupled 
with  an  air  tube  and  an  AIl-TTnion  State  Standard  63?6-52>  fyPe  ® 
vane  anemometer. 

Analysis  o^  the  Field  of  Velocity  in  Average  Velocities 

The  study  of  the  field  of  velocity  was  conducted  with  steady, 
average  (5»  10,  15  m/sec),  and  low  (0.4,  0.6,  1.2,  and  2.0  m/sec) 
velocities . 

The  working  field  of  the  channel  in  average  velocities  was 
studied  by  ordinary  methods  along  the  axis  of  both  air  tubes. 

One  of  them  (cone  shaped  calibrating)  was  mounted  securely  in 
the  center  of  the  working  field,  and  the  second  (spherical, 
Prandtl-type) ,  located  in  a  simplified  coordinator,  could  be 
moved  to  any  point  in  the  selected  cross  section  of  the  stream. 

The  readings  of  the  velocity  pressures  in  every  individual 
point  of  the  field  were  taken  with  the  above-mentioned  TsAGI 
micromanometer,  which  was  connected  to  the  movable  tube  at  an 
angle  of  inclination  of  about  8°  (sin  a  =  0.125).  The  readings 
of  the  velocity  pressures  in  the  center  of  the  working  field 
were  taken  simultaneously  by  the  second  calibrating  micromanometer 
of  the  same  type,  which  was  connected  to  the  stationary  calibrating 
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Figure  2.  Deviations  of  the  field's  coefficients  from  the 
center  (%),  velocity  5  m/sec. 


Figure  3.  Deviations  of  the  field's  coefficients  from  the  center  (%) 
Velocity:  a  =  10  m/sec;  b  =  15  m/sec. 
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tube 


During  the  processing,  the  readings  of  the  movable  tube 
were  reduced  to  the  readings  of  the  calibrating  tube,  located 
in  the  center,  with  the  help  of  a  transfer  coefficient,  which 
was  obtained,  first  of  all,  in  the  results  of  the  movable  tube 
with  respect  to  the  calibrating  tube,  and  of  verification  of  the 
micromanometers  which  were  used  in  pair  with  the  tubes. 

The  velocity  distribution  was  analyzed  in  one  cross  section, 
going  through  the  center  of  the  working  field  of  the  tunnel,  and 
along  several  contour  lines  spaced  10  cm  from  each  other. 

It  is  necessary  to  note  that,  during  the  measuring  of  the 
field  along  the  vertical  diameter  below  the  center,  an  adverse 
effect  of  the  stationary  tube  holder  upon  the  readings  of  the 
movable  tube  was  noticed.  Therefore,  instead  of  the  stationary 
calibrating  tube  being  used  for  the  measuring  of  the  velocity 
in  the  center  of  the  field,  the  orifice  in  the  collector, 
graduated  according  to  the  calibrating  tube,  was  used.  This 
method  of  measuring  the  velocity  is  used  with  the  checking  of 
meteorological  instruments,  and  is  mentioned  in  the  pertinent 
handbook^ . 

The  uniformity  of  the  velocity  fields  was  characterized  by 
the  values  of  the  coefficients  obtained  as  the  ratio  Vf>  ,  where 
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v  is  the  velocity  of  the  air  stream  at  individual  points  of 
the  field  with  respect  to  the  movable  tube,  and  where  vu  is 
the  simultaneously  measured  velocity  in  the  center  of  the 
field  with  respect  to  the  calibrating  tube. 

It  is  possible  to  visualize  the  accuracy  of  the  analysis 
of  the  tunnel's  field  with  given  velocities  by  the  root  mean 
square  deviation  a  of  a  series  of  coefficients  obtained  in  the 
center  of  the  working  field  of  the  tunnel  during  the  calibration 
of  the  movable  tube  by  the  calibrating  tube,  as  is  shown  in 
Table  I. 

It  can  be  seen  from  Table  I  that  the  root  mean  square 
deviation  of  the  series  of  coefficients,  characterizing  the 
accuracy  of  the  coefficients  obtained  through  analysis  of  the 
field  in  given  velocities,  consists  of  from  0.6  to  0.9 fo. 

The  results  of  the  analysis  of  the  velocity  field  in  average 
velocities  are  shown  in  Table  II  (see  1952),  and  are  illustrated 
by  the  diagrams  in  Figures  2  and  3*  These  diagrams  give  the 
deviations  of  the  coefficients  for  individual  points  of  the 
field  from  the  coefficient  at  the  center,  which  was  taken  as 


unity,  in  percentages 


■  fc-') 


100^. 


It  follows  from  Table  2  and  the  diagrams,  that  the 


coefficients  in  the  working  part  of  the  tunnel's  field  (in  a 
girth  of  40  x  80  cm)  with  given  velocities  are  changed 
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Table  2 


A  -  value  of  the  coefficients  of  the  field  along  the  horizontal  diameter 
B  -  Spacing  of  the  points  in  the  stream  from  the  axis  of  the  tunnel,  m. 
C  -  counterclockwise  along  the  streams 
D  -  Month  and  Year 
E  -  Velocity,  m/sec 

F  -  Spacing  of  the  points  in  the  stream  from  the  axis  of  the  tunnel,  m. 
G  -  Month  and  Year 
H  -  Center 

I  -  Clockwise  along  the  stream 
J  -  Velocity,  m/se.c 

K  -  Value  of  the  coefficients  of  the  field  along  the  vertical  diameter 
L  -  Spacing  of  the  points  in  the  stream  from  the  axis  of  the  tunnel,  m. 
M  -  Upwards  from  the  axis  of  the  tunnel 
N  -  Month  and  Year 
O  -  Velocity,  m/sec 

P  -  Spacing  of  the  points  in  the  stream  from  the  axis  of  the  tunnel,  m. 
Q  -  Center 
R  -  Month  and  Year 

S  -  Downwards  from  the  axis  of  the  tunnel 
T  -  velocity,  m/sec 


insignificantly,  within  limits  of  V%.  Maximum  deviations 
(3  to  4  rfo)  of  the  coefficients  from  the  center  are  observed  only 
at  the  extreme  points  of  the  field. 

During  the  determination  of  the  coefficients  of  the  field 
with  the  help  of  the  graduated  orifice  in  the  tunnel’s  collector, 
the  uniformity  of  the  velocities  in  the  working  field  is  maintained 
within  the  very  same  limits.  This  makes  it  possible  to  analyze 
the  tunnel’s  field  not  by  two  air  tubes,  which  causes  a  shading 
at  certain  points  of  the  field,  but  by  one  tube,  where  this 
presents  the  possibility  of  using  the  calibrated  orifice  in  the 
collector  as  a  control  point  in  the  center  of  the  working  field. 

The  coefficients  of  the  field,  which  were  obtained  during 
its  analysis  in  1937  and  1941,  are  shown  in  Table  II  as 
confirmation  of  the  characteristic  of  the  working  field,  obtained 
in  1952. 

A  relative  comparison  of  the  coefficients,  obtained  in 
different  years,  shows  that  these  coefficients  diverge  among 
themselves  by  1  to  2$.  Divergences  in  3 $  are  single  cases. 

This  makes  it  possible  to  come  to  the  conclusion  that,  for  all 
practical  purposes,  the  distribution  of  the  velocities  in  the 
cross  section  of  the  channel,  passing  through  the  center  of  the 
working  field,  can  be  considered  as  constant  in  all  examined 
velocities  within  the  limits  of  accuracy  of  the  experiment. 


On  the  basis  of  such  a  conclusion,  which  was  arrived  at 
in  1952,  other  cross  sections  of  the  tunnel  (at  a  distance  of 
0.3  ra  from  the  center  plane  and  heading  in  the  direction  of  the 
diffusor  and  collector)  were  not  analyzed.  They  were  examined 
in  the  preceding  years,  and,  with  respect  to  the  uniformity 
of  the  distribution  of  the  velocities,  they  differed  very  little 
from  the  cross  section  passing  through  the  center  of  the 
tunnel . 

An  .Analysis  of  the  Tunnel's  Field  in  Small  Velocities 
Two  slide-wire  axled  vane  anemometers  from  the  MoscobU 
hydrometric  instrument  factory  were  used  for  an  analysis  of  the 
velocity  field  in  low  velocities.  The  anemometers  made  it 
possible  to  measure  the  mean  velocity  of  a  directed  air  stream. 

The  sensitivity  threshold  of  the  anemometers  consisted  of 

0.2  m/sec,  and  the  upper  limit  of  measurement  consisted  of  5  m/sec. 

A  simplified  coordinator  in  the  form  of  a  rectangular 
frame  with  movable  clamps,  which  was  to  be  located  outside  of 
the  tunnel's  field,  was  needed.  To  fill  this  requirement,  we 
strung  out  a  copper  wire  at  the  necessary  distances.  The  vaned 
anemometer,  with  its  holding  clamp  loosened,  could  easily 
displace  to  various  points  of  the  tunnel's  working  field  along 
this  wire. 
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Because  the  size  of  the  anemometer  prevented  the  mounting 
of  a  second  instrument  as  a  regulating  anemometer  into  the 
center  of  the  working  field,  the  velocity  of  the  stream  in 
the  center  was  measured  by  one  and  the  same  anemometer  -  which 
was  periodically  placed  in  the  center  of  the  working  field  - 
through  2  to  3  measurements  at  individual  points  of  the  field. 

The  indications  of  the  anemometer  were  read  after  holding 
it  for  100  seconds  at  the  examined  point  of  the  field. 

In  other  respects,  the  measurements  were  conducted  along 
the  same  plan  which  was  used  for  the  analysis  of  the  field  in 
average  velocities.  The  coefficients  were  calculated  in  an 
analogous  manner. 

In  view  of  the  fact  that  the  sensitivity  threshold  of  the 
anemometers  was  close  to  0.2  m/sec,  it  was  safe  to  compute  the 
measurements  only  from  0.4  m/sec.  Therefore,  the  field-  was 
measured  at  velocities  of  0.4,  0.7,  1.0,  and  2.0  m/sec. 

In  order  to  control  the  results  obtained  at  0.4  m/sec,  the 
measurements  were  made  twice  by  two  different  anemometers. 

The  representation  concerning  the  accuracy  of  measprement 
by  vaned  anemometers  gives  in  Table  III  the  root  mean  square 
deviation  of  a  series  of  velocities,  which  were  obtained  by 
the  readings  of  two  anemometers  periodically  placed  in  the 
center  of  the  working  field  of  the  tunnel  during  the  process 
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A  -  average  magnitude  of  velocity,  m/sec 
B  -  root  mean  square  deviation  (+) 

C  -  absolute,  m/sec 
D  -  relative,  % 

E  -  Number  of  readings 


Figure  4.  Deviations  of  the  field's  coefficients 
from  the  center  (%),  velocity  0,4  m/sec 
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of  analysis. 

In  conformity  with  Table  III,  the  root  mean  square 
deviation  of  the  measurement  of  velocity  by  slide-wire  axled 
anemometers  does  not  exceed  +  0.01  ra/sec  in  low  velocities. 
Therefore,  they  are  fully  adequate  for  comparative  measure¬ 
ments  . 

The  results  of  the  analysis  of  the  field  in  low  velocities 
are  shown  in  Table  IV  and  in  the  diagram  of  Figure  4. 

From  a  comparison  of  the  diagrams  in  Figures  3  and  4,  it 
follows  that  in  low  velocities  of  the  order  of  0.4  /sec,  the 
field  of  velocity  in  the  center  part  (cross  section  40  x  80  cm) 
is  less  uniform  than  at  average  velocities,  particularly  at  the 
edges.  There  is  a  significant  decrease  of  velocity  in  the  upper 
part  of  the  cross  section.  The  deviations  of  the  coefficients 
at  individual  points  of  the  field  away  from  the  center  amount  to 
from  52  to  60^.  Inversely,  an  increase  in  velocities  can  be 
observed  in  the  lower  part  of  the  cross  section. 

In  this  particular  case,  the  coefficients,  as  compared  with 
the  center,  are  greater  by  from  5  to  7^.  A  repeated  examination 
of  the  field  at  a  velocity  of  0.4  m/sec  by  a  second  anemometer 
confirmed  the  results  obtained  the  first  time.  The  field 
gradually  straightens  out  at  velocities  greater  than  0.4  m/sec, 
and,  as  can  be  seen  from  Table  IV,  beginning  with  a  velocity  of 
1.2  m/sec,  the  coefficients  differ  among  themselves  by  no  more 
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than  2%,  which  is  within  the  limits  of  the  accuracy  of 
measurement  of  the  velocity  by  a  vaned  anemometer.  Exceptions 
are  the  extreme  points  of  the  field  up  from  the  axis  of  the 
tunnel,  where  there  is  still  a  decrease  of  velocity  in 
comparison  with  the  center. 

It  can  be  assumed  that  the  effect  of  a  difference  of 
temperatures  in  the  closed  chamber  and  in  the  casing  of  the 
tunnel  manifests  itself  in  low  velocities.  This  difference  of 
the  air  stream  at  low  velocities  produces  an  additional 
velocity  field  in  the  working  volume,  owing  to  the  temperature 
heterogeneities,  which  bring  about  a  different  density  of  the 
air  in  the  working  field  of  the  tunnel. 

Hence  follows  the  practical  conclusion  that  the  (H40  wind 
tunnel  cannot  be  suitable  for  the  checking  of  vaned  anemometers 
at  velocities  below  that  of  0.4  m/sec,  on  account  of  the  fact 
that  the  indicated  difference  of  temperatures  at  these  velocities 
will  cause  a  still  greater  non-uniformity  in  their  distribution. 

Velocity  Fluctuations  of  the  Air  Stream 
Important  characteristics  of  the  nature  of  the  air  stream 
are  the  fluctuations  of  the  velocity,  which  are  always  observed 
to  a  greater  or  lesser  degree  around  a  certain  average  value. 

The  inconstancy  of  the  voltage  in  the  network  and  the 
natural  oscillations  of  the  closed  chamber  under  the  influence 


of  the  stream  can  he  noted  as  causes  of  these  fluctuations. 


The  magnitude  of  the  fluctuation  of  the  air  stream  was 
determined  in  the  process  of  calibrating  two  air  tubes,  placed 
in  the  center  of  the  working  field  of  the  tunnel,  by  the 
fluctuation  of  the  alcohol  column  in  a  TsAGI  micromanometer 
(sin  a  =  O.I25)  at  various  velocities  of  the  stream  -  6,  10, 

12,  aY\d  15  m/sec.  Ten  instantaneous  indications  were  read 
simultaneously  by  both  micromanometers  in  each  of  the  indicated 
velocities.  A  micromanometer  was  connected  to  each  of  the  two 
air  tubes. 

The  absolute  value  of  the  velocity  fluctuation  was  character¬ 
ized  by  the  magnitude 

-  v„ 


max 


'mm 


m/ sec 


relative  to 

vmax  ~  vmin 

v  +  v 
max  mi n 

The  results  are  shown  in  Table  5* 

Due  to  the  fact  that  the  data  cited  in  Table  5  includes  the 
natural  oscillations  and  inertia  of  the  micromanometer,  it  is  not 
sufficiently  accurate,  and  can  be  used  only  for  a  qualitative 
judgement  concerning  the  fluctuations  of  the  air  stream. 

The  velocity  fluctuations  of  the  air  stream,  as  is  evident 
from  the  table,  do  not  exceed  0.1  m/sec  in  average  velocities. 
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Table  5 


(fi)  Ch'OpoCTb,  M/ceK. 

(§)  3HaqeHHe  Koae6annii  (+) 

i 

Mhcjio 

iiaO.ito.jcmm 

Cg) 

MdKCHMd/IbHdfl 

to 

MHHHM3  JlbHdfl 
(?) 

aCcoaioTHoe, 
(g}  M/ceK. 

OTHOCHTeJlbHOe. 
(?)  °/» 

5,92 

5,96 

5,84 

0,06 

1,0 

160 

9,81 

9,87 

9,71 

0,08 

0,8 

160 

11,84 

11,93 

11,77 

0,08 

0,7 

160 

14,96 

15,06 

14,84 

0,11 

0,7 

160 

A  -  velocity,  m/sec 

B  -  magnitude  of  fluctuation  (+  ) 

C  -  maximum 

D  -  minimum 

E  -  absolute,  m/sec 

F-  relative,  % 

G  -  number  of  observations 
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If  we  are  to  take  Into  consideration  the  fact  that,  in  practice, 
it  is  proposed  to  take  several  readings  hy  the  mi croma nometer 
in  order  to  eliminate  random  errors  owing  to  these  fluctuations, 
then  it  is  possible  to  assume  that  the  effect  of  the  fluctuations 
upon  the  velocity  measurement  ip  insignificant. 

Analysis  of  the  Direction  of  the  Stream 

It  was  proposed  to  use  a  special  cylindrical  nozzle, 
which  was  produced  by  the  "Etalon"  factory,  for  an  analysis  of 
the  direction  of  the  air  stream  in  the  working  field  of  the 
tunnel.  However,  this  nozzle,  in  view  of  the  insufficient 
length  of  its  adapter,  did  not  make  it  possible  to  measure  the 
direction  without  disturbing  the  stream  itself.  Therefore, 
the  analysis  was  limited  only  to  the  solution  of  the  problem 
concerning  the  capability  of  applying  the  nozzle  in  low 
velocities,  i.e.,  for  a  practical  determination  of  its  sensitivity. 

This  cylindrical  nozzle,  similar  to  the  one  described  in 

4 

Popov’s  work  ,  has  three  orifices,  of  which  the  two  openings  on 
the  edges  are  symmetrically  located  at  angles  of  45°  froin  the 
third,  or  center,  orifice. 

The  direction  of  the  stream  is  determined  by  the  two  outer 
openings,  which  are  connected  to  the  micromanometer.  The 
operating  principle  of  the  nozzle  is  based  on  the  fact  that  the 
pressures  in  the  outer  openings  are  equal  if  they  are  symmetrical 
with  reference  to  the  stream.  The  center  orifice  of  the  nozzle, 


under  conditions  of  its  calibration  by  the  calibrating  air  tube, 
serves  for  tbe  determination  of  the  stream  velocity  in  that 
plane  where  the  direction  is  also  determined.  For  that  reason, 
one  of  the  outer  orifices,  after  its  orientation  in  the  air 
stream,  is  connected  to  the  micromanometer  in  addition  to  the 
center  orifice.  The  orientation  of  the  nozzles  with  respect 
to  the  air  stream  was  done  as  follows:  They  were  placed  with 
the  orifices  facing  against  the  stream,  and  rotated  about  their 
axis  in  one  or  the  other  direction  until  that  time  when  the 
pressures  in  the  outer  points  were  not  equal.  The  indicated 
position  of  the  nozzle  at  the  control  point  of  the  stream  Was 
taken  as  the  origin  (deflection  angle  of  the  nozzle  a  =  0°) 
and  was  fixed  by  a  reading  along  the  graduated  circle  of  the 
nozzle.  The  deflection  angle  up  to  this  position  also  determines 
the  direction  of  the  stream  in  the  given  points  in  relation  to 
the  direction  at  the  control  point,  which  was  taken  as  the 
initial  point. 

In  order  to  determine  the  sensitivity  of  the  nozzle,  the 
outer  orifices  were  connected  to  the  micromanometer  at  an  angle 
of  inclination  a.  =  3°* 

The  nozzle,  having  been  oriented  in  the  stream  with  respect 
to  direction,  i.e,  having  attained  equal  pressures  in  the  two 
outer  openings,  was  turned  at  different  angles  of  ratotation  in 
both  directions,  and  the  corresponding  indications  of  the 
micromanometer  were  read  in  each  instance.  The  results  of  the 
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Table  6 


CKOpOCTb 
B03jyiuHoro 
noTOKa,  M/ceK. 

CpeAHee  n3MCHeHiie  nokaaaHnii  MiiKpoxaHOkieTpa 
e)  upii  noBopoTC  wacaaK-a  Ha'l° 

Mnc.no 

iiaS.iio^eiiHH 

® 

B  AeJICHKRX 
(Q  MHKpOM3HOMCTpa 

B  MHJIJIIlMCTpaX 
vi/BOilflHoro  CTo^6a 

0,9 

0,1 

0,004 

■ 

12 

'  2.1 

0,4 

0,016 

6 

5,0 

3,3 

0,146 

6 

10,0 

12,6 

0,559 

4 

A  -  velocity  of  air  stream,  m/sec 

B  -  average  change  of  indication  of  micromanometer  with  a 
rotation  of  the  nozzle  by  1° 

C  -  in  micromanometer  divisions 

D  -  in  mm  of  water  column 

E  -  number  of  observations 


determination  of  the  sensitivity  of  the  nozzle  are  listed  in 
Table  VI. 

It  follows  from  Table  VI  that  the  micromanometer  ,  with  a 
rotation  of  the  nozzle  in  1°  in  low  velocities,  changes  itR 
indication  in  tenths  on  the  scale  of  division.  Therefore,  the 
measurement  of  the  direction  at  velocities  below  that  of  2  m/sec 
will  be  accompanied  by  larger  errors,  and  it  is  not  appropriate 
to  use  the  described  method. 

By  generalizing  the  results  of  the  analysis  of  the  GGO 
wind  tunnel  in  relation  to  the  field  of  velocity  and  its 
operating  characteristics,  we  may  come  to  the  conclusion  that, 
regardless  of  its  obsolete  construction,  it  still  belongs  to 
that  group  of  first-rate  wind  tunnels  where  it  is  possible,  with 
the  presence  of  appropriate  measuring  apparatus,  to  conduct 
accurate  experiments  and  to  conduct  a  checking  of  instruments 
in  velocity  limits  from  0.5  to  20  m/sec.  An  intrinsic  defect 
of  the  tunnel  is  the  low  upper  limit  of  the  produced  velocities. 
This  defect  became  particularly  felt  after  the  creation  of  new 
instruments  during  the  past  few  years  which  can  measure 
velocities  up  to  40  m/sec  and  above.  In  order  to  expand  the 
capability  of  this  wind  tunnel,  it  is  desirable  to  increase 
the  upper  velocity  limit,  or  to  construct  a  new  tunnel. 
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TESTING  OF  HAND  ANEMOMETERS  AT  THE  PLANT 

This  article  discusses  shortly  the  results  of  work  which  made  it  possible 
to  produce  hand  cup  anemometers  GOST  6376-52,  type  A  that  are  provided  with  a 
standard  curve  at  the  plant. 

Until  1958,  hand  cup  anemometers  GOST  6376-52,  type  A  manufactured  at 
the  Moscow  Hydrometeorological  Instrument  Plant,  before  being  shipped  away, 
used  to  be  tested  in  the  wind  tunnel  of  the  plant.  Four  anemometers  were 
mounted  in  the  cross-section  of  the  tunnel.  As  a  result  of  the  test  for  each 
anemometer,  a  calibration  curve  was  plotted  which  made  it  possible  to  express 
the  instrument's  reading  as  an  actual  air  stream  velocity.  Each  anemometer 
to  be  shipped  out  was  provided  with  such  a  calibration  curve. 

(Information  on  the  anemometer  as  well  as  details  on  its  conventional 
testing  are  set  forth  in  the  Specifications  /l /  now  in  effect). 

In  view  of  the  standardized  manufacturing  of  hand  anemometers  and  of  the 
small  difference  between  individual  calibration  curves,  the  Moscow  Hydrometeoro¬ 
logical  Instrument  Plant  proposed  that  anemometers  be  produced  that  are  provided 
with  a  single  standard  calibration  curve. 

In  order  to  implement  this  proposal,  this  writer  and  plant  workers  con¬ 
ducted  a  special  joint  investigation  in  the  plant  control  shop.  It  was  revealed, 
first  of  all,  that  the  most  accurate  standard  calibration  curve  is  obtainable  by 
testin^fnemometers  in  an  UPAR  —  an  installation  for  the  testing  of  hand  anemo¬ 
meters  —  which  is  also  described  in  Specifications  /l / .  In  the  UPAR  the  test 
error  is  less  than  in  the  other  aerodynamic  installations  of  the  plant.  This 
is  because  in  the  UPAR,  anemometers  are  tested  singly  and  not  in  groups,  the 
UPAR  is  little  affected  by  changing  external  conditions  (temperature  and 
pressure) ,  and  a  mechanical  adjustment  ensures  uniform  mounting  of  the  hand 


anemometers  themselves. 

These  circumstances  are  dealt  with  in  detail  in  the  paper  / 2/ . 

To  derive  a  standard  calibration  curve  129  calibration  curves  for  stand¬ 
ard  anemometers  manufactured  before  1958  were  used.  All  the  curves  were  obtain¬ 
ed  from  UPAR  tests  of  these  instruments. 

Analysis  of  these  curves  showed  that  a  single  calibration  curve  will  not 
suffice  to  produce  accurate  instruments,  since  test  results  of  individual  instru¬ 
ments  may  deviate  substantially  from  this  curve.  Therefore,  all  the  instruments 
were  divided  on  a  sensitivity  basis  into  three  groups  according  to  the  character 
of  their  curves.  In  the  first  group  were  placed  9  instruments  whose  highest 
readings  (given  20  divisions  per  second)  were  18.1  to  18.5  m/sec.;  in  the  second 
group,  70  instruments  with  the  maximum  readings  18.6  to  19.0  m/sec.;  and  in  the 
third  group  were  placed  50  instruments  with  a  maximum  velocity  reading  of  19.1 
to  19*5  m/sec. 

Standard  curves  for  each  group  were  obtained  by  taking  the  arithmetic 
mean  of  the  values  read  from  the  calibration  curve  of  each  anemometer  for  the 
actual  air  stream  velocity;  these  values  correspond  to  readings  on  the  scale  at 
every  other  division  per  second  (l  to  20  divisions  per  second).  Thus  were 
obtained  three  standard  curves.  The  mean  maximum  velocity  proved  to  be  18.4 
m/sec  for  the  first  standard  curve,  given  anemometer  readings  of  20  divisions 
per  second,  18.8  m/sec  for  the  second,  and  19.2  m/sec  for  the  third. 

To  characterize  a  number  of  calibration  curves  which  had  been  used  for  the 
plotting  of  individual  standard  cruves,  the  root-mean-square  error  £T  was 
calculated  for  the  values  read  from  the  standard  curve  for  the  actual  air  stream 
velocities  which  correspond  to  readings 

on  the  anemometer  scale  of  1,  10,  and  18  divisions  per  second.  Talble  1  gives 
the  results  obtained. 
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9 

2 
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0,05 

3,7 
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70 

3 

1,36 

0,05 

! 

3,7 

9,84 

0,08 

0,8 

17,37 

0,14. 

0,8 
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As  Table  1  shows,  the  root-mean-square  error  of  a  number  of  calibration 
curves  for  individual  anemometer  groups  does  not  exceed  +  0.14  m/sec  or  +  0.8 
of  the  velocity  in  the  vicinity  of  maximum  velocities. 

Table  2  shows  the  distribution  of  deviations  of  individual  calibration  curves 
from  the  mean  (standard)  curve  which  is  characterized  by  the  distribution  of 
deviations  of  individual  actual  velocity  values  as  read  from  individual  calibra¬ 
tion  curves,  from  the  arithmetic  means  for  each  group  of  instruments. 

. 

Table  2 


Distribution  of  deviations  of  individual  calibration  curves  from  the  mean 
(standard)  curve  in  percents  of  the  whole  number 
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Total 

number 

•<  lo 

<2a 

<3a 

1 

1,31 

78 

100 

100 

9 

9,40 

66 

100 

100 

9 

16,59 

78 

89  • 

100 

9 

2 

1,35 

76 

96 

100 

70' 

9,65 

73 

97 

100 

70 

17,02 

68 

99 

100 

70 

3 

1,36 

78 

92 

100 

50 

9,84 

70 

96 

100 

50 

17,37 

68 

94 

100 

50 

/a  C, 


3 


From  the  data  of  Table  2,  the  overwhelming  majority  of  deviations  of 
individual  curves  from  the  standard  curve  fall  within  the  range  of  twice  the 
root-mean-square  error  and  amounts  to  89  to  99$  for  all  the  anemometer  groups 
at  velocities  close  to  the  upper  limit. 

In  this  context,  anemometers  whose  individual  test  data  will  deviate  from 
the  standard  calibration  curve  by  more  than  2  &  can  be  considered  not  to  con¬ 
form  to  our  series,  and  cannot  be  provided  with  the  standard  curve  that  we 
obtained. 

For  any  anemometer  reading,  the  numerical  expression  for  twice  the  root- 
raean-square  error  (we  have  in  mind  the  highest  error  obtained  for  a  group  of  3 
anemometers)  will  be  approximately  (0.1  +  0.01  v)  m/ sec,  where  v  is  the  actual 
air  stream  velocity. 

The  quantity  (0.1  +  0.01  v)  m/sec  is  taken  as  a  permissible  variation  of 
an  individual  calibration  curve  from  the  most  suitable  standard  curve.  This 
permissible  variation  is  less  than  the  error  in  a  conventional  wind  tunnel  group 
test  of  anemometers. 

Evidently,  when  we  arrange  individual  calibration  curves  between  standard 
curves,  and  then  substitute  the  most  suitable  standard  curve  for  each  individual 
curve,  we  are  introducing  an  error  not  exceeding  0.2  m/sec. 

It  should  be  noted  that  by  increasing  the  number  of  standard  curves  to 
five  or  more  in  lieu  of  the  recommended  three,  it  would  seem  at  first  glance 
that  we  should  accomplish  better  conformity  of  the  standard  curve  to  the 
individual  characteristic  of  an  anemometer.  Such  a  gain,  however,  would  be  but 
an  apparent  one.  Actually  the  bringing  closer  together  of  velocity  calibration 
curves  through  an  increase  in  their  number,  is  limited  by. the  error  magnitude 
of  each  individual  anemometer.  If  this  error  exceeds  deviations  in  velocity 
between  selected  standard  curves,  there  will  be  no  further  gain  in  accuracy 
and  there  is  no  point  in  making  use  of  such  standard  curves  that  are 

closer  together.  On  the  other  hand,  it  follows  from  the  data  considered  above 
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that  a  single  standard  curve  is  not  enough. 

In  order  to  test  whether  the  obtained  standard  curves  agree  with  individual 
test  data,  the  results  of  UPAR  control  tests  of  412  hand  anemometers,  taken 
from  the  plant  files  were  used.  From  readings  taken  on  these  instruments  at 
the  highest  velocity  at  which  maximum  deviations  from  standard  curves  are  most 
likely,  the  most  suitable  calibration  curve  was  selected  for  each  instrument. 

The  deviations  proved  to  be  as  follows:  115  anemometers  (2?. 9$)  fell  along  the 
selected  standard  curve,  216  (52.4$)  showed  deviations  of  0.1  m/sec,  63  (15 • 3$) 
of  0.2  m/sec,  11  (2.7$)  of  0.3  m/sec,  and  7  (1.7$)  of  more  than  0.3  m/sec. 

Thus,  the  overwhelming  majority  of  instruments  could  be  provided  with  a 
standard  curve. 

Based  on  this  work,  a  Procedure  was  written.  In  conformity  with  this 
Procedure,  hand  anemometers  GOST  6376-52,  type  A  are  provided  with  a  standard 
curve  upon  shipment. 

All  instruments  produced  are  tested.  The  test  includes  an  external  exa¬ 
mination,  a  determination  of  the  sensitivity  threshold,  and  a  recording  of  two 
anemometer  readings  in  the  UPAR.  Headings  are  checked  at  two  velocities:  at 
any  speed  within  the  range  1  to  5  m/sec  and  at  17  to  20  m/sec. 

A  anemometer  which  has  met  GOST  specifications  as  to  external  examination 
and  testing,  is  provided  with  the  most  suitable  standard  curve.  Such  a  curve 
is  selected  from  test  data  of  anemometer  readings  by  comparison5 ?hese  data 
with  each  of  three  standard  curves.  Readings  should  coincide  with  the  selected 
curve  or  deviate  from  it  by  no  more  than  (0.1  +  o.l  v)  m/sec. 

Instruments  which  meet  GOST  specifications  but  deviate  from  standard  curves 
by  a  value  greater  than  indicated  above,  are  provided  with  individual  calibra¬ 
tion  curves  obtained  by  the  full  testing  of  these  instruments  in  conformity  with 
the  current  Specifications  for  the  Testing  of  Meteorological  Instruments  /l/. 

The  above-described  production  method  for  hand  anemometers  has  been  in  effect 

at  the  Moscow  Hydrometeorological  Instrument  Plant  since  1958.  The  method  makes 
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it  possible  to  effect  substantial  savings. 


and  to  produce  instruments  that  are  more  comparable  to  each  other,  since  the 
Control  Bureau  of  UGMS  (Local  Administration  of  the  Hydrometeorological  Service) 
is  also  testing  anemometers  in  UPAR's  -which  have  been  calibrated  in  the  GGO 
(Main  Geophysical  Observatory)  wind  tunnel. 

Analysis  and  generalization  of  the  results  of  this  method  at  the  plant 
will  make  possible  the  future  control  of  output  uniformity  and  quality  from 
the  distribution  of  anemometers  among  the  standard  curves,  as  well  as  the  use 
of  this  method  for  testing  other  types  of  anemometers. 
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ERRORS  OF  THERMOELECTRIC  BALANCOMETER  TESTING 

This  article  contains  the  results  of  experimental  investigations 
on  errors  with  which  the  parameters  of  balancometers  used  in  hydro¬ 
meteorological  network  stations  are  determined  in  testing. 

The  investigation  of  error  in  testing  thermoelectric  balanco¬ 
meters,  just  as  much  similar  work,  being  carried  on  by  the  Central 
Testing  Laboratory,  GGO ,  is  intended  to  evaluate  errors  with  which 
the  parameters  of  balancometers  used  in  hydrometeorological  network 
stations  are  determined  in  testing.  Accidental  errors  of  the  con¬ 
version  factor,  balancometer  sensitivity  and  wind  correction  factors 
were  evaluated.  In  addition,  balancometer  sensitivity  dependance 
on  the  intensity  and  direction  of  incident  radiation,  on  the  tempera¬ 
ture  of  the  medium  and  on  the  direction  of  the  air  flow  striking  it 
was  investigated. 

Ten-battery  balancometers,  prepared  in  the  GGO  Experimental  Work¬ 
shop  in  1951,  1956  and  1957,  were  studied.  The  balancometers  satisfy 
the  requirements  imposed  on  them  in  testing  with  the  exception  of 
balancometer  No.  31A2;  the  sensitivities  of  its  sides  differ  by  10% 
(permitted  tolerance  5%). 

The  accidental  errors  were  evaluated  as  the  mean  quadrantal  error 
of  a  number  of  values  obtained  in  repeated  balancometer  tests.  Test¬ 
ing  was  carried  out  in  conformity  with  the  Manual  on  testing  (l). 

By  one  testing  the  entire  cycle  of  measurements  carried  out  according 
to  the  Manual  i s  understood. 

ACCIDENTAL  ERROR  OF  THE  CONVERSION  FACTOR  AND  SENSITIVITY 

In  order  to  evaluate  accidental  error  in  conversion  factor  and 
sensitivity  determination^ 6  balancometers  were  tested.  3  galvano¬ 
meters,  prepared  by  the  GGO  Experimental  Workshop  and  the  Leningrad 


Hydrometeorological  aevice  Factory,  were  used  with  the  balancometers. 
A  thermoelectric  actinometer  served  as  a  model  in  testing.  The  tests 
were  repeated  ten  times. 

The  conversion  factors  of  balancometers  coupled  with  gal vanometers 
were  determined  by  all  the  methods  recommended  in  the  Manual  on 
testing,  namely:  (l)  under  laboratory  conditions  with  radiation  ob¬ 
tained  from  a  lamp,  (2)  under  natural  conditions  with  an  installation 
in  which  the  detection  surfaces  of  the  balancometer  are  placed  per¬ 
pendicular  to  the  incident  radiation  and  protected  from  the  action 
of  scattered  radiation  and  wind,  and  (3)  under  natural  conditions 
without  an  installation  with  the  detection  surfaces  of  the  balanco¬ 
meter  placed  horizontally  just  as  when  observations  are  made. 


The  sensitivity  of  balancometers  (without  coupled  galvanometer) 
was  determined  only  under  laboratory  conditions. 

The  conversion  factors  of  the  devices  were  determined  under  labo¬ 
ratory  conditions  at  the  installation  TS.L.P.  in  November,  1957. 

The  mean  conversion  factor  values  obtained  for  each  device  couple 
from  ten  testings  and  the  mean  quadrantal  errors  of  a  number  of 
measurements  for  the  sides  of  balancometer  No.  1  are  listed  in 
Table  1.  The  mean  quadrantal  errors  of  the  conversion  factors 
calculated  for  each  of  the  sides  and  for  the  mean  values  are  given 
in  Table  2. 

The  results  given  in  Table  2  indicate  that  the  mean  quadrantal 
error  of  the  conversion  factor  does  not  exceed  2.2%.  The  mean 
quadrantal  error  of  the  conversion  factor  (the  average  for  two  sides) 
calculated  for  all  balancometers  tested  (repeated  60  times)  is  1.2%. 


/// 


Table  1 


Mean  conversion  factor  values  a  and  mean  quadrantal  error  o  (  %  ) 
obtained  in  testing  under  laboratory  conditions 


Instrument 

3413 

3429 

3142 

3324 

3136 

3169 

No. 

37 

37 

1178 

1178 

781 

781 

a 

0,0140 

0,0138 

0,0162 

0,0150 

0,0184 

0,0166 

a 

1.8 

1.2 

0,9 

1,2 

0,6 

1.2 

Table  2 

Mean  quadrantal  errors  a  of  the  conversion  factors  of  the  balancometers  % 


Instrument 

No. 

3413 

3429 

3142 

3324 

3136 

3169 

acp 

37 

1178 

781 

781 

Side  1 

1,8 

1,2 

0,9 

1,2 

0,6 

1.2  . 

1.2 

Side  2 

2,2 

1,1 

1,5 

1.3 

0.8 

1.0 

1.4 

average 

2,2 

1,0 

1.1 

1,1 

0,6 

0,9 

1,2 

Table  3 


Mean  quadrantal  error  o  of  balancometer  sensitivity  % 


Instrument 

No. 

3136 

3159 

3114 

3142 

3137 

3169 

Side  1 

1,0 

1,3 

1,0 

1,6  ' 

0,8 

1,4 

Side  2 

1,0 

1,0 

1.5 

2,1 

1,2 

1.7 

average 

0,8 

0.9 

1.1 

1,8 

1.4 

The  distribution  of  remaining  errors  is  as  follows:  errors,  not 
exceeding  O'  constitute  70%,  those  not  exceeding  2CT ,  93%  and  those 
not  exceeding  3  O'  ,  100 %  of  all  the  errors  involved. 

For  evaluating  errors  in  determination  of  balancometer  sensitivity, 
the  results  of  balancometer  tests  repeated  ten  times  at  the  instal¬ 
lation  TS.B.P.  in  1951  were  used. 

The  mean  quadrantal  errors  of  sensitivity  calculated  for  each  of 
the  sides  and  for  the  mean  sensitivity  values  are  given  in  Table  3. 


It  is  evident  from  the  data  listed  in  Table  3  that  the  mean  quad¬ 
rantal  error  of  sensitivity  does  not  exceed  2.1%.  The  mean  quadran¬ 
tal  error  calculated  for  sensitivity  values  (the  average  for  two 
sides)  on  all  tested  balancometers  is  1.2%. 

The  distribution  of  remaining  errors  is  as  follows:  errors  not 
exceeding  CT  constitute  7$%,  those  not  exceeding  aler,  95%  and  those 
not  exceeding  3  tr ,  100%  of  all  the  errors  involved. 

Thus,  accidental  errors  in  determining  conversion  factors  of 
balancorneter-galvanometer  couples  and  in  determining  balancometer 
sensitivity  have  one  and  the  same  order. 

Balancometer  conversion  factors  were  determined  for  natural  con¬ 
ditions  during  the  expedition  on  Lake  Sevan  in  September,  1957.  The 
testing  of  the  devices  both  with  and  without  installations  was  re¬ 
peated  ten  times.  Because  of  the  large  amount  of  work  involved  in 
testing  under  natural  conditions  the  conversion  factors  were  deter¬ 
mined  for  only  one  of  the  sides  (No.  1)  of  each  balancometer. 

The  mean  conversion  factor  values  obtained  as  a  result  of  10 
tests  in  natural  conditions  in  an  installation  and  the  mean  quadran- 
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tal  errors  of  a  number  of  them  are  listed  in  Table  4. 

From  the  data  given  in  Table  4  it  is  evident  that  'the  mean  quad- 
rantal  errors  in  the  determination  of  conversion  factors  under  natur¬ 
al  conditions  in  an  installation  do  not  exceed  2.3 %.  The  mean  quad- 
rantal  error  calculated  from  the  results  of  tests  on  all  balancometers 
(repeated  60  times)  is  1.6%.  The  distribution  of  remaining  errors 
is  as  follows:  errors,  not  exceeding  o"  ,  constitute  8 2 %,  those  not 
exceeding  2 & >  95 %  and  those  not  exceeding  3 or,  100%  of .  all  the 
errors  involved. 

In  determining  balancometer  conversion  factors  under  natural  con¬ 
ditions  without  an  installation  the  test  was  conducted  with  the  sun 
at  an  altitude  of  not  less  than  40°.  The  balancometers  were  placed 
so  that  their  detection  surfaces  were  horizontal  and  the  handle  would 
be  directed  toward  the  side  away  from  the  sun.  The  wind  velocity 
during  the  test  varied  from  0  to  3  meters  per  second.  The  wind  was 
intermittent . 


From  the  data  given  in  Table  5  it  is  evident  that  the  .mean  quad- 
rantal  error  does  not  exceed  2.3%.  The  mean  auadrantal  error  cal¬ 
culated  from  the  results  of  all  balancometer  tests  is  1.7%.  The 
distribution  of  remaining  error  is  as  follows:  errors,  not  exceeding 
Of,  constitute  80%,  those  not  exceeding  2  <F ,  98%  and  those  not  ex¬ 
ceeding  3 cr  >  100%  of  all  the  errors  involved. 

The  mean  conversion  factor  for  side  No.  1  and  the  mean 

ouadrantal  errors  obtained  in  the  determination  of  conversion  factors 
by  different  methods  are  given  in  Table  6. 


•  Table  4 

Mean  conversion  factor  a  values  and  mean  quadrantal  errors  a  (  %  )  obtained 
in  testing  under  natural  conditions  in  an  installation 


Instrument 

3413 

3429 

'3142 

3324 

3136 

3169 

No. 

37 

37 

1178 

1178 

781 

781 

a 

0,0140 

0,0137 

0,0164 

0,0149 

0,0181 

0,0163 

<3 

.  1,1 

0,9 

2,3 

1,0 

1.6 

2.0 

Table  5 

Mean  conversion  factor  a  values  and  mean  quadrantal  errors  a  (%)  obtained  in 
testing  under  natural  conditions  without  an  installation 


Instrument 

3413 

3429  . 

3142 

3324 

3136 

3169 

No. 

37 

37 

1178 

1178 

781 

781 

a 

0,0141 

0,0137 

0,0163 

0,0149 

0,0182 

0,0162 

a  # 

1,6 

1,8 

•  1,3 

1,3 

.1,5 

2.3 

Table  6 

Balancometer  conversion  factors  obtained 'by  various  methods  and  the  mean 

quadrantal  errors  a,  %  _ 


_ _  Instrument  No. 

Method  3413 


In  the  laboratory,  th 
Under  field  f 
conditions  1 


in  the  in-  6>6i46 
stallation 
without  the  0,0141 
installation 


3429 

3142 

3324 

3136 

3169 

37 

1178 

1178 

781 

781 

0,0138 

0,0137 

0,0162 

0,0164 

0,0150 

0,0149 

0,0184 

0,0181 

0,0166 

0,0163 

0,0137 

0,0163 

0,0149 

0,0182 

1 0,0162 

//S' 


The  compounding  of  mean  conversion  factor  values  determined  by 
different  methods  (Table  6)  indicates  that  divergence  between  values 
determined  by  different  methods  is  within  the  limits  of  accidental 
error  and  does  not  have  a  definite  sign.  This  confirms  the  possibil¬ 
ity  of  testing  devices  by  any  one  of  the  demonstrated  methods.  The 
smallest  accidental  error  value  was  obtained  while  testing  under  lab¬ 
oratory  conditions  as  might  be  expected  since  radiation  currents 
under  laboratory  conditions  are  more  constant.  It  must  be  remembered 
that  in  testing  in  the  laboratory  10  pair  of  comparative  readings 
were  taken,  while  under  natural  conditions  20  pair  were  taken.  In 
testing  under  natural  conditions  without  an  installation  one  reading 
was  understood  as  being  the  average  of  three  readings  taken  one  im¬ 
mediately  after  the  other  as  is  the  custom  in  observation. 

For  evaluation  of  accidental  error  in  the  individual  balancometer 
readings  for  each  balancometer,  the  mean  quadrantal  error  of  the  con¬ 
version  factor  determined  by  one  pair  of  balancometer  readings  (in 
sun,  in  shade)  is  calculated.  The  values  obtained  are  given  in 
Table  7. 

The  distribution  of  remaining  individual  reading  errors  in  con¬ 
version  factor  determination  which  were  calculated  from  test  results 
on  all  balancometers  is  given  in  Table 

The  compounding  of  magnitudes  of  errors,  given  in  Table  7, 
indicates  that  individual  reading  error  in  conversion  factor  determi¬ 
nation  under  laboratory  conditions  is  considerably  less  than  indivi¬ 
dual  reading  error  under  natural  conditions.  Thus,  it  proves  the 
necessity  of  taking  a  large  number  of  readings  in  testing  under  natu¬ 
ral  conditions  if  an  installation  is  not  used  for  testing. 


Table  7 

Mean  quadrantal  errors  of  the  individual  reading  in  determining  the  conversion 

factor,  % 


Method 

3413 

'37' 

3429 

"37 

EH 

°cp 

In  the  laboratory  ^  the  in_ 

1  ‘  l.T 

1,4 

1,0  . 

1.4 

0,9 

1,3 

1 

1,3 

Under  field  /  stallation 

2,5 

2,4 

2,9 

2.4 

2,9 

2,9 

2,7 

conditions  l  without  the  3,2  • 

installation 

3,4  ' 

3,8 

3.2 

3,3 

4,3 

3,5 

Table  8 

Distribution  of  remaining  errors  of  the  individual  reading  in  determining  the 

the  conversion  factor,  % 


No.  of 
cases 

a 

2s 

3s 

78 

97 

100 

600 

>n  SI 

98 

100 

1200 

81 

97 

99 

1208 

Method 


In  the  laboratory 

Under  field 
conditions 


\  without  the  in¬ 
stallation 


The  distribution  of  remaining  errors  given  in  Table  8  indicates 
that  practically  all  deviations  do  not  exceed  double  the  mean • ouadran- 
tal  error  which  in  testing  under  laboratory  conditions  is  2.6$,  in 
testing  under  natural  conditions  in  an  installation,  5.4$  and  in 
testing  under  natural  conditions  without  an  installation,  7$.  Thus, 
in  determining  balancometer  conversion  factors,  values  differing 
from  the  mean  value  by  more  than  26  should  be  disregarded. 

BALANCOMETER  SENSITIVITY  DEPENDANCE  ON  RADIATION  INTENSITY 

In  studying  balancometer  sensitivity  dependance  on  radiation 
intensity  for  4  balancometers ,  the  sensitivity  was  determined  at 
radiation  intensities  of  1.5,  0.75  and  0.15  cal/cm^  per  minute. 

The  source  of  radiation  was  a  1000  kilowatt  projector  lamp  and  the 
flow  of  radiation  was  condensed  by  means  of  two  lenses.  The  inten¬ 
sity  of  the  radiation  was  varied  by  means  of  a  rotating  disc  with 
sectors  removed.  This  disc  was  proposed  by  B.  A.  Eisenstadt  (2). 
Balancometer  sensitivity  determination  was  carried  out  by  the  met¬ 
hod  ordinarily  used  in  testing  with  one  exception:  that  illumination 
of  the  devices  was  effected  by  means  of  a  nonparallel  beam.  In  the 
installation  used  in  TS.L.P.  for  device  testing,  the  intensity  of 
parallel  beams  of  radiation  obtained  from  1000  watt  lamps  is  0.6- 
0.7  cal/cm^  per  minute.  In  order  to  obtain  larger  radiation  inten¬ 
sity  values  pencils  of  light  were  created.  This  led  to  the  fact 
that  surface  illumination  of  the  balancometer  was  uneven;  the  inten¬ 
sity  diminished  toward  the  edges.  Because  of  this,  the  mean  intensi¬ 
ty  values  of  the  radiation  striking  the  surface  of  the  balancometer 
and  actinometer  are  not  identical  and  the  absolute  balancometer 
sensitivity  values  do  not  precisely  correspond  to  the  true  values. 
However,  since  illumination  distribution  in  a  given  experiment 

remains  constant,  it  is  possible  to  compound  the  sensitivity  values 


obtained  with  different  radiation  intensity  values.  With  each 
radiation  intensity  value,  the  sensitivity  of  the  balancometer  was 
determined  three  times  during  which  time  the  mounting  of  the  balanco¬ 
meter  and  actinometer  on  the  rotary  table  of  the  installation  remained 
unchanged.  The  balancometer  sensitivity  values  in  millivolts  at 
1  cal/cm^  per  minute  obtained  with  radiation  intensities  of  1.5,  0.75, 
and  0.15  eal/cm^  per  minute  are  given  in  Table  9.. 

It  is  evident  from  the  data  given  in  Table  9  that  divergences  in 
balancometer  sensitivity  values  obtained  with  different  radiation 
values  are  found  to  be  within  the  limits  of  magnitude-  of  accidental 
errors  in  measuring.  Therefore,  it  can  be  assumed  that  within  radi¬ 
ation  intensity  limits  from  0.15  to  1.5  cal/cm^  per  minute,  balanco¬ 
meter  sensitivity  does  not  depend  on  radiation  intensity. 

Tn  addition,  the  sensitivities  of  three  balancometers  (No.  3324, 
3142  and  3160)  at  radiation  intensities  of  0.3  and  3  cal/cm2  per 
minute  wore  compared.  Tt  appeared  from  this  that  with  an  increase 
of  radiation  intensity,  balancometer  sensitivity  increases.  If 
balancometer  sensitivity  at  a  radiation  intensity  of  0.3  cal/cm^  per 
minute  is  taken  to  be  100%,  the  sensitivity  of  balancometer  No.  3324 
at  a  radiation  of  3  cal/cm^  per  minute  increased  by  5^,  of  balanco¬ 
meter  No.  3142,  2%  and  that  of  balancometer  No.  3169,  6%.  The  re¬ 
sults  obtained  contradict  the  theoretical  calculations  of  K .  Y.  Kon¬ 
dratiev  (3)  according  to  which  balancometer  sensitivity  diminishes 
with  radiation  balance  increase. 


Table  9 


o 

Balaneometer  sensitivity  in  millivolts  at  1  cal/ cm  per  minute  with  different 

values  of  radiation  intensity _ 


Instrument 

No. 

,  Radiation  intensity,  cal/  cm  -min 

1.5 

0,75 

0,15 

3324 

7,04 

6,94 

6,93 

!  .  1 

7,06 

6,95 

6,93 

7,04 

6,96 

6,8S 

average  7,05 

6,95 

6,91 

3142 

6,81 

6,78 

6,92 

V 

6,85 

6,81 

7,02 

6,86 

6,84 

6,88 

average  6,84 

6,81 

6,94 

3169 

9,32 

9,16 

9,14 

9,32 

9,21 

9,14 

9,30 

9,20 

9,15 

average  9,31 

9,19 

9,14 

3429 

6,70 

6,67 

6,71 

6,69 

6,64 

6,68 

6,68 

6,58 

6,72 

average  6,69  1 

6,63 

6,70 

BALANCOMETER  SENSITIVITY  DSPENDANCE  ON  DIRECTION  OF  INCIDENTAL 

RADIATION 

Balaneometer  sensitivity  dependance  on  the  direction  of  incidental 

radiation  was  investigated  under  laboratory  conditions. 

„  „  „  iNDicAT/aJVS 

the  correction  factors  for  balaneometer 

at  various  angles  of  incidence  were  determined  for  6  balanco- 
meters.  In  this  case  the  balancometers  were  installed  so  that  at  an 
angle  of  incidence  of  90°  (at  an  elevation  of  the  sun  of  0°),  the 
handle  of  the  balaneometer  was  directed  tbward  the  side  opposite  to 
the  source  of  light.  The  correction  factors  were  determined  twice. 
Tn  addition,  for  one  of  the  balancometers  (No.  3169 )j  the  correction 
factors  for  dependance  on  direction  of  incidental  radiation  were  de¬ 
termined  also  by  rotating  the  balaneometer  over  the  azimuth  by  90°. 


The  correction  factor  values  obtained  for  dependance  on  the  elevation 
of  the  sun  are  given  in  Table  10.  The  correction  factors  for  balan¬ 
cometer  No.  3169,  determined  for  two  azimuths  differing  by  90°,  are 
listed  in  the  last  two  lines  of  Table  10. 

From  the  data  given  in  Table  10  it  is  evident  that  balancometer 
sensitivity  depends  on  the  direction  of  incidental  radiation;  here 
the  correction  factor  for  exposures  of  the  balancometer  is  always 
greater  than  unity.  Sensitivity  dependance  on  the  direction  of  in¬ 
cident  radiation  is  given  for  elevations  of  the  sun  of  50°  or  less, 

therefore,  for  all  practical  purposes,  of  the  balancometer 

are  always  somewhat  '  The  data  obtained  confirm  the  inexpedi¬ 

ency  of  measurement  with  an  unshaded  balancometer,  especially  at  small 
elevations  of  the  sun. 


Table  10 

Correction  factors  for  indications  of  the  balancometer  depending  on  the 


elevation  of  the  sun 


V 

Instrument 

No. 

Sun's  altitude,  degrees 

2,5 

5 

10 

20 

30 

40 

50 

•  3413 

1,86 

1,17 

1,15 

1,09 

1,06 

1,05 

1,03 

3429 

1,20 

1,18 

1,11 

1,06 

1,02 

1,02 

1,02 

3142 

1,19 

1,16 

1,08 

1,04 

1,03 

1,02 

1 ,02 

3324 

1,26 

1,20 

1,17 

1,08 

1,06 

1,04 

1,02 

3136 

1,73 

1,28 

1,12 

1,05 

1,03 

1,02 

1 ,01 

3169 

1,45 

1,28 

1,18 

1,08 

1,04 

1,02 

1,01 

3169 

1,50 

1,41 

•  1 '  17 

1,06 

1,02 

1,01 

1,00 

70 


Remarks 


1,00 

1,00 

1,01 

1,00 

1,00 

1,00 

1,00 


Rotated  along 


the  azimuth 


by  90° 


ERRORS  IN  WIND  CORRECTION  FACTOR  DETERMINATION 

For  the  evaluation  of  accidental  errors  of  wind  correction  factors, 
Ehe  factors  were  determined  in  a  wind  tunnel  ten  times  for  7  balanco- 
meters.  In  this  case  the  correction  factors  were  determined  for  each 
of  the  sides  for  two  of  the  balancometers ,  for  the  others,  over  one 
side  as  required  by  the  Manual. 

The  mean  wind  correction  factor  values  at  wind  velocities  of  1, 

5,  10  and  15  meters  per  second  are  given  in  Table  11.  In  Table  12 
the  mean  quadrantal  error  of  correction  factor  determination  at  wind 
speeds  of  5,  10  and  15  meter/second,  expressed  in  percent  of  correc¬ 
tion  factor  size  are  listed. 


Table  11 

Mean  values  of  wind  correction  factors 


Instrument 

No. 

Velocity 

7,  m/sec 

1 

5 

10 

15 

3413 

1,02 

1 ,09 

1,17 

1,24 

3429  (CTOpOHa  1) 

1,02 

1,06 

1,11 

1,16 

3429  (  .  2) 

1,02 

1,08 

1,14 

.1,21 

3142  (  .  1) 

1 ,04 

1 ,16 

1,28 

1,38 

3142  (  .  2) 

1 ,04  ' 

1,16 

1,29 

1,39 

3324 

1,02 

1,08 

1,15 

1,22 

3136 

1 ,02 

1,10 

1,19  ' 

1,27 

3169 

1,02 

1,11 

1,20 

1,27 

3204 

1,02 

1  ,07 

1,12 

1,17 

The  mean  quadrantal  errors,  calculated  from  all  balancometer  test 
results,  are:  0.01  (1.3%)  for  a  wind  speed  of  5  meters  per  second, 
0.02  (1.7%)  for  a  wind  speed  of  10  meters  per  second  and  0,2  (1.7%) 


for  a  wind  speed  of  15  meters  per  second. 


A?  2. 


The  distribution  of  remaining  errors  is  as  follows:  for  a  wind 
speed  of  5  meters  per  second,  the  errors  not  exceeding  c  constitute 
73%,  those  not  exceeding  2  CT  ,  100%;  for  a  wind  speed  of  10  meters 
per  second,  the  errors  not  exceeding  or  constitute  79%  and  those  not 
exceeding  2  O',  100%;  for  a  wind  speed  of  15  meters  per  second,  errors 


Table  12 

Mean  quadrantal  errors  of  wind  correction  factors  of  balancometers  (  %  )  at 
at  wind  speeds  of  5,  10  and  15  meters  per  second 


Instrument 

No. 

velocity,  m/sec 

5 

10 

15 

3413 

!  '  1,3 

1,8 

1.9 

3429 

1.3 

2,0 

1,8 

3142 

1,3 

1.7 

1,8 

3324 

.0,7 

1.5 

1,6 

3136 

1,8 

2.1 

2,0 

3169 

1,3 

1,5 

1.4 

3204 

1,2 

1,5 

1.5 

not  exceeding  <r  ,  constitute  80%,  those  not  exceeding  2  0*  ,  9 8%  and 
those  not  exceeding  3 O'  ,  100%. 

The  correction  factors  were  determined  for  each  of  the  sides  for 
balancometers  No.  3429  and  3142.  According  to  the  data  listed  in 
Table  11  it  is  evident  that  the  correction  factors  for  both  sides  of 
balancometer  No.  3142  practically  coincide  while  for  balancometer 
No.  3429  they  differ  by  0.05  (4%)  at  a  wind  velocity  of  15  meters  per 
second.  It  must  be  noted  that  the  sensitivitigg  of  the  sides  of 
balancometer  No,  3142  differ  by  10%  while  the  sensitivities  of  the 


sides  of  balancometer  No.  3429  differ  by  no  more  than  1.5%. 

The  results  outlined  above  were  obtained  for  wind  tunnel  condi¬ 
tions  in  which  the  balancometer  was  installed  (by  eye)  so  that  the 
detecting  surfaces  were  parallel  to  the  airflow;  at  the  same  time 
the  detecting  surfaces  were  located  vertically  and  illuminated  from 
the  side. 

IN£>ICA7 

As  is  known,  balancometer  under  natural  wind  conditions 

are  very  unstable  and  as  a  result  it  is  possible  to  determine  correc¬ 
tion  factors  under  natural  conditions  only  from  a  large  number  ef 
readings.  The  compounding  of  wind  correction  factors  in  determining 
them  under  natural  conditions  and  in  the  wind  tunnel  (4)  indicated 
that  the  mean  correction  factor  values  obtained  from  a  large  number 
of  determinations  made  under  natural  conditions  correspond  well  with 

the  values  obtained  in  wind  tunnel  testing.  One  of  the  causes  of 
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balancometer  inconstancy  at  the  same  wind  velocity  may  be 

a  difference  in  the  degree  of  air  turbulence  and  the  fact  that  the 
air  flow  is  not  always  directly  parallelly  to  the  detecting  surfaces 
of  the  balancometer^  but  strikes  them  at  different  angles. 

In  order  to  evaluate  the  influence  of  the  direction  of  the  air 
flow  striking  the  balancometer,  wind  correction  factors  were  deter¬ 
mined  for  three  balancometers  with  different  orientations  of  the  bal¬ 
ancometer  in  the  air  current  of  the  wind  tunnel.  The  detection  sur¬ 
faces  of  the  balancometers  were  situated  in  the  airflow  at  different 
angles  of  attack  so  that  the  current  of  air  was  directed  first 

on  the  dark  and  then  on  the  illuminated  detection  surface.  The  cor¬ 
rection  factors  were  determined  at  angles  of  incidence  of  2 ,  5,  10 
and  30°.  In  these  experiments  the  balancometers  No.  3413  and  3324 
were  arranged  so  that  their  detection  surfaces  were  always  situated 


in  vertical  planes  and  illuminated  from  the  side.  However,  balan- 
cometer  No.  3572  was  first  situated  horizontally  and  then  tilted 
depending  on  the  angle  of  incidence  chosen.  It  was  illuminated 
from  above.  In  addition,  the  correction  factors  for  balancometer 
No.  3572  at  0°  angle  of  incidence  were  determined  with  the  detec¬ 
tion  surfaces  in  a  vertical  position  (as  is  done  in  ordinary  testing). 
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Wind  correction  factors  for  balancometer  at  different 

angles  of  incidence  and  with  air  current  velocities  of  2,  5,  10  and 
15  meters  per  second  are  given  in  Table  13.  On  the  last  line  of 
Table  13  the  correction  factors 

for  balancometer  No.  3572  when  the  detection  surfaces  are  placed  ver¬ 
tically  are  given. 

Graphs  of  wind  correction  factors  for  balancometer  No.  3324  at 
different  orientations  in  the  air  flow  are  given  in  Figure  1.  For 
better  visualization,  the  curves  obtained  with  different  angles  of 
incidence  are  drawn  on  the  general  coordinates  of  the  axis.  The 
spped  of  the  air  flow  in  raeters  per  second  is  drawn  along  the  ab¬ 
scissa  axis  and  the  correction  factors  along  the  ordinate  axis.  The 
curves  for  cases  whene  the  air  current  was  directed  on  the  illuminated 
side  of  the  balancometer  are  represented  by  broken  lines.  The  curves 
for  cases  where  the  air  flow  was  directed  upon  the  nonilluminated 
side  of  the  balancometer  are  represented  by  broken  lines  with  dots. 

The  curve  for  the  case  where  the  balancometer  is  located  lengthwise 
with  the  flow  is  represented  by  a  solid  line  (from  the  average  value 
of  ten  testings) 

The  results  given  in  Table  13  and  in  the  graphs  (Fig.  l)  indicate 
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Fig.  1.  The  wind  correction  factors  for  balancometer  No.  3324  at  different  orientations 

in  the  flow  of  air 


Wind  correction  factors  of  balancometers 


Table  13 


4-> 

a 

<D 

<D 

O 

Stream  velocity  (m/ sec),  directed  to  the  side 

m  © 

illuminated 

dark 

h 

T3  ©  £ 

•g  MM  9 

5 

10 

15 

2 

5 

10 

15 

&  . 

a  a-s 

3413 

■  0 

1,04 

1,09 

1.17 

1,24 

2 

1,04 

1,09 

1,16 

1,21 

1,03 

1,07 

1,12 

1,16 

5 

1,04 

1,09 

1,16 

1,22 

1,02 

1,02 

1,02 

1.02 

10 

1,06 

1.14 

1,26 

1,36 

1,02 

1,05 

1,08 

1,09 

30 

1,06 

1,13 

1,23 

1,30 

1,05 

1,10 

1,12 

1.13 

3324 

0 

1,03 

1,08 

1,15 

1,22 

5 

1,05 

1,12 

1,20 

1,26 

1,02 

1,04 

1,05 

1,06 

10 

1,06 

1,14 

1,26 

1,37 

1,02 

1,03 

1 ,04 

1,04 

30 

1,04 

1,10 

1,19 

1,28 

1,05 

1,08 

1,11 

1,13 

3572 

0 

1,03 

1,07 

1,14 

1,21 

' 

5 

1,06 

1,13 

1,22 

1,26 

1,02 

1,04 

1,05 

1,07 

10 

1,08 

1,16 

1,28 

1,36 

1,03 

1,06 

1,09 

l.H 

28 

1,04 

1,08 

1,14 

1.19 

1,04 

1,07 

1.12 

1,15 

0 

1,04 

1,09 

1,17 

1,23 

>• 
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that  wind  correction  factors  depend  essentially  on  the  orientation 
of  the  balancometer  in  the  air  current. 


The  correction  factors  are  decreased  if  the  flow  of  air  is  directed 
on  the  illuminated  side  of  the  balancometer  and  are  increased  if  it 
is  directed  upon  the  nonilluminat ed  side.  This  indicates  that  the 
side  opposite  to  the  one  upon  which  the  flow  of  air  is  directed, 
i.e.,  the  side  upon  which  the  turbulent  air  flow  strikes  is  cooled 
to  a  greater  degree.  The  maximum  correction  factor  values  were  ob¬ 
served  at  angles  of  incidence  of  10°,  the  minimum  values  at  angles 
of  incidence  of  5  and  10°.  At  an  angle  of  incidence  of  10°  depen¬ 
ding  on  whether  the  air  flow  is  directed  on  the  illuminated  or  the 
non-illuminated  side,  the  correction  factors  differ  by  8  to  10$  at 
an  air  current  velocity  of  5  meters  per  second  and  by  15  to  19$  at 

an  air  current  velocity  of  10  meters  per  second.  Thus,  the  large 
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variations  in  balancometer  observed  under  natural  condit¬ 

ions  may  be  partially  attributed  to  the  fact  that  the  balancometer 

is  subjected  to  the  action  of  randomly  directed  air  currents.  How- 
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ever,  balancometer  variation  may  not  be  caused  entirely  by 

variance  of  the  airflow  striking  the  balancometer  because,  in  testing 
with  an  installation  where  the  balancometer  is  protected  from  wind 
influence,  the  mean  auadrantal  error  of  the  individual  reading 
(Table  7)  also  has  a  large  value. 

Correction  factors  obtained  at  angles  of  incidence  of  2°  differ 
only  blightly  from  those  obtained  in  an  installation  with  the  balan¬ 
cometer  lengthwise  to  the  current.  It  is  possible  by  this  to  explain 
the  comparatively  small  differences  in  correction  factor  values  obtained 
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in  repeated  testing  of  balancometers  installed  along  the  current  by 
sight . 

The  results  given  in  Table  13  for  balancometer  No.  3572  for  an 
angle  of  incidence  of  0°  indicate  that  the  correction  factors  ob¬ 
tained  with  horizontal  and  vertical  positioning  of  the  detection 
surfaces  practically  coincide. 

BALANCOMETER  SENSITIVITY  DEPENDANCE  ON  TEMPERATURE 

Balancometer  sensitivity  dependance  on  the  temperature  of  the 
medium  has  already  been  studied  (5)  for  three  balancometers  in  which 
the  space  between  the  detection  plates  and  thermobattery  was  filled 
with  rosin.  It  was  established  that  balancometer  sensitivity  increa¬ 
ses  with  a  rise  in  temperature.  At  the  present  t ime shellac  instead 
of  rosin  is  used  in  the  manufacture  of  balancometers.  It  was  useful 
to  study  these  balancometers  bv  the  relationship  of  temperature  coef¬ 
ficient  of  their  sensitivities. 

The  temperature  coefficients  of  balancometer  sensitivities  were 
determined  within  a  temperature  range  of  3  to  4^°  according  to  the 
methods  described  in  the  treatise  (5).  Here  two  balancometers  filled 
with  rosin  (Nos.  3136  and  3169). for  which  the  temperature  coefficients 
had  been  previously  determined,  and  two  balancometers  filled  v/ith 
shellac  (Nos.  3324  and  3429)  were  taken.  The  results  of  determina¬ 
tions  on  the  temperature  coefficients  of  sensitivity  are  given  in 
Table  14. 

The  data  given  in  Table  14  confirms  the  temperature  coefficient 
values  obtained  previously  for  balancometers  filled  with  rosin.  For 


Table  14 

Results  of  determining  the  temperature  coefficient  of  balancometer  sensitivity 


Instrument  No. 

Date 

3136 

3169 

3324 

3429 

25  II 

26  III 

6  III 

18  III 

0.0C09 

0,0009 

0,0013 

0,0014 

0,0014 

0,0008 

0,0006 

0,0015 

0,0012 

-0,0001 

-0,0002 

0,0000 

+0,0003 

Average 

0,0010 

0,0014 

0,0010 

0,0000 

balancometers  filled  with  shellac,  the  temperature  coefficient 
obtained  for  the  sensitivity  of  balancometer  No.  3324  is  of  the 
sane  order  as  that  for  balancoraeters  filled  with  rosin.  Balanco¬ 
meter  No.  3429,  on  the  other  hand,  does  not  have  a  temperature 
coefficient . 

CONCLUSIONS 

1.  Balancometer  conversion  factors  obtained  by  repeated  testings 
under  laboratory  and  natural  conditions  (by  two  methods)  coincide 
very  well.  This  confirms  the  possibility  of  conversion  factor  deter¬ 
mination  by  any  one  of  the  three  methods  recommended  by  the  Manual 
for  testing  meteorological  devices. 

2.  It  is  suitable  to  test  balancomet ers  under  laboratory  condit¬ 
ions  because  here  accidental  errors  of  the  conversion  factor  have  a 
smaller  value.  This  method  also  requires  less  work. 

The  mean  auadrantal  error  in  conversion  factor  determination  is 
1.2 %  in  testing  under  laboratory  conditions,  1.6$.  in  testing  under 
natural  conditions  with  an  installation  and  1.7#  in  testing  under 
natural  conditions  without  an  installation. 

3.  The  mean  quadrantal  error  of  each  individual  reading  in  bal¬ 
ancometer  testing  is  1.3#  in  testing  under  laboratory  conditions, 

2.7%  under  natural  conditions  with  an  installation  and  3.5#  under 
natural  conditions  without  an  installation. 

The  large  accidental  error  -values  of  individual  readings  in 
testing  balancometers  under  natural  conditions  confirm  the  necessity 
of  performing  a  large  number  of  comparative  readings  (not  less  than 
20  pair) . 

4.  Within  the  limits  of  radiation  balance  magnitudes,  encountered 
under  natural  conditions,  balancometer  sensitivity  does  not  depend 

on  the  magnitude  of  measured  radiation  balance. 
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5.  Balancometer  sensitivity  depends  on  the  direction  of  the  in¬ 
cident  radiation:  with  increasing  angle  of  incidence  balancometer 
sensitivity  decreases.  In  connection  with  this,  the  value  of  the 
radiation  balance,  measured  by  a  balancometer,  is  always  lowered. 

In  measuring  with  unshaded  balancomet ers  it  is  necessary  to  intro¬ 
duce  correction  factors  individual  for  each  balancometer. 

6.  The  mean  quadrantal  error  in  determining  wind  correction  fac¬ 
tors  in  a  wind  tunnel  is  1.3 %  at  an  air  flow  velocity  of  5  meters  per 

second  and  1.7%  at  an  air  flow  velocity  of  10  and  15  meters  per  second. 

GAT  ta/tJ 

Balancometer  in  wind  depend  essentiallv  on  the  direc¬ 

tion  of  the  air  currents  striking  the  balancometer.  Therefore,  the 

correction  factors  obtained  in  the  wind  tunnel  are  applicable  only 
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for  correcting  averaged  balancometer 
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